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ABSTRACT 


The  objective  of  this  program  was  to  prepare  a  clutter  research  plan 
for  Army  air  defense  weapons  that  encompasses  the  anticipated  methods  of 
employment  of  these  weapons,  as  well  as  the  problems  associated  with  the 
introduction  of  new  radar  technology  in  these  systems.  Specifically,  the 
investigations  included  studies  of  the  physics  of  electromagnetic  scattering 
from  rough  surfaces,  a  definition  of  clutter  problems  as  related  to  innovations 
in  radar  technology  such  as  electronic  beam  steering,  signal  processing  and 
design,  computer  technology,  study  of  clutter  models  and  simulation  techniques, 
and  a  study  of  clutter  effects  in  a  radar  defense  complex. 

Probability  density  functions,  correlation  functions  and  doppler 
spectra  of  clutter  signals  are  derived  for  the  configuration  in  which  only  a 
few  major  scatterers  are  present.  Results  are  presented  for  both  coherent  and 
noncoherent  detection. 

Methods  of  simulating  clutter,  using  both  analog  and  digital 
approaches,  are  described.  Emphasis  is  placed  on  advanced  forms  of  simulators 
such  as  would  be  required  for  evaluation  of  future  multi-dimensional  radar 
systems  rather  than  on  additive  filtered  noise  for  the  simulation  of  clutter. 

The  simulation  techniques  are  based  on  two  suggested  generic  forms  of  clutter 
models,  called  open-loop  and  closed-loop.  Because  of  the  greater  flexibility 
of  closed-loop  simulation,  this  is  presented  in  greater  detail. 

Techniques  for  performing  radar  clutter  measurements  are  described 
and  some  of  the  significant  problems,  related  primarily  to  terrain  and  cultural 
observables,  are  discussed.  Clutter  problems  as  related  to  radar  systems  and 
innovations  in  radar  technology,  such  as  electronically  steered  phase  arrays, 
are  also  described. 

The  program  plan  for  future  clutter  research  incorporates  the  findings 
of  the  supporting  studies  performed  on  this  program.  The  program  plan  identifies 
the  major  areas  of  development  and  research  required  to  perform  a  meaningful 
program  of  clutter  measurements.  It  is  stressed  that:  1)  the  electromagnetic 
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observables  should  be  measured  in  such  a  way  that  correlation  is  maintained 
among  the  measured  parameters  and  2)  that  environmental  descriptors  are  of 
importance  equal  to  that  of  the  radar  sensor  data  and  must  be  derived  and 
associated  with  it. 
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Section  1 

INTRODUCTION  AND  SUMMARY 


1.1  INTRODUCTION 

The  objective  of  this  program  was  to  prepare  a  clutter  research  plan 
for  Army  air  defense  weapons  that  encompasses  the  anticipated  methods  of 
employment  of  present  weapons,  as  well  as  the  problems  associated  with  the 
introduction  of  new  radar  technology  in  these  systems.  Emphasis  was  placed 
on  terrain  clutter;  other  sources  of  clutter  such  as  sea  and  precipitation 
were  not  treated  explicitly. 

This  task  was  approached  by  identifying  several  specific  topics 
which  are  logical  elements  of  a  clutter  research  program  and  then  by  examin¬ 
ing  each  topic  in  turn  in  its  relation  to  the  others.  The  topics  selected 
for  this  purpose  are:  1)  physics  of  electromagnetic  scattering  from  rough 
surfaces,  2)  definition  of  clutter  problems  as  related  to  innovations  in 
radar  technology  such  as  electronic  steering,  signal  processing  and  design, 

3)  computer  technology,  4)  clutter  models  and  simulation  techniques,  5)  clutter 
effects  in  a  radar  defense  complex.  Each  one  is  discussed  in  the  corresponding 
section  of  the  report  and  their  interrelation  to  form  a  well  integrated  clutter 
research  program  plan  is  described  in  Section  10. 

Clutter  has  been  defined  as  a  conglomeration  of  unwanted  echoes. 

From  a  phenomonological  view,  a  more  accurate  and  fundamental  description 
would  be  the  presence  of  unwanted  scatterers  or  collections  of  scatterers 
within  the  system  environment.  It  is  true  that  these  unwanted  scatterers 
produce  unwanted  echoes.  However,  since  observed  echoes  are  strong  functions 
of  both  illumination  and  observation  systsms,  basic  phenomena  cannot  be  de- 
cribed  in  such  terms.  For  example,  observed  echoes  (clutter)  are  a  function 
of  the  resolution  of  illumination  and  observation  systems.  The  overall 
environment  has  not  changed,  just  the  observation  technique.  This  change  in 
emphasis  removes  the  particular  radar  system  from  consideration  and  views  the 
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clutter  problem  only  in  terms  of  the  basic  phenomena  of  scattering.  This 
point  is  of  utmost  importance  when  simulation  of  a  number  of  radar  types  is 
required  because  changes  in  radar  observables  for  an  environment  in  which  the 
scatterers  (and  their  time  variations)  are  defined  can  be  traced  directly  to 
changes  in  radar  system  parameters.  Moreover,  use  of  a  suitably  defined 
scattering  model  permits  evaluation  of  clutter  effects  on  such  diverse  systems 
as  CW  Doppler  (for  personnel  intrusion  detection)  and  missile  semiactive 
homing  systems.  The  former,  of  course,  requires  a  definition  of  scatterer 
motion  and  cross  section  (amplitude  and  phase),  while  the  latter  requires,  in 
addition,  definition  of  target-scatterer  interactions. 

In  general,  effects  of  clutter  can  bo  broken  into  three  broad  cate¬ 
gories.  The  first  is  the  monostatic  return  from  clutter  itself  where  clutter 
radar  cross  section  is  sufficiently  high  to  obscure  the  true  target  return. 

The  second  is  modification  of  radar  system  parameters  by  clutter-radar  inter¬ 
action,  and  the  third  is  modification  of  true  target  return  by  target -clutter 
interaction.  Examples  of  each  are  readily  available.  For  the  first,  we  have 
the  difficulty  of  discerning  a  low-altitude,  low-radar-cross-section  vehicle 
in  terrain  clutt«r,  particularly  with  a  low-resolution  radar.  For  the  second, 
we  have  the  modification  of  antenna  lobe  structure  and  boresight  of  a  ground- 
based  radar  operating  close  to  grazing  incidence  which  can  increase  the 
apparent  target  scintillation,  or  under  conditions  of  local  terrain  masking 
can  interrupt  the  line  of  sight  between  the  radar  and  target.  For  the  third, 
we  have  the  difficulty  of  differentiating  between  a  target  and  its  image  with 
the  antenna  beam  close  to  terrain.  (Usually  the  Doppler  frequency  difference 
or  range  resolution  is  insufficient  for  effective  discrimination.) 

These  observed  clutter  effects  are  functions  of  radar  system  para¬ 
meters  as  well  as  the  electromagnetic  environment,  and  generally  depend  on 
system  resolution,  both  in  range  and  angle  and  on  range/angle  scan  rates.  In 
fact,  it  is  precisely  this  dependence  on  system  parameters  which  limits  the 
usefulness  of  much  of  the  experimental  clutter  data  obtained  in  the  past. 

While  much  data  exist  which  are  valid  for  particular  systems  employed  to  make 
the  measurements,  they  cannot  be  extrapolated  to  predict  the  performance  of 
systems  with  different  capabilities. 
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Another  factor  of  major  importance,  which  has  all  too  often  been 
neglected,  is  an  adequate  description  of  the  environment  in  which  a  set  of 
clutter  measurements  has  been  made.  Data  on  environmental  variables,  such 
as  terrain  cover  and  masking,  is  of  equal  importance  with  electromagnetic 
variables  for  the  interpretation  and  understanding  of  clutter.  The  major 
problem  in  selecting  environmental  descriptors  is  that  a  set  sufficient  to 
categorize  terrain  with  respect  to  clutter  has  not  been  defined.  This  factor 
stresses  the  need  for  care  in  the  design  and  conduct  of  a  clutter  measurements 
program.  The  emphasis  during  such  a  program  should  be  primarily  on  gathering, 
analyzing  and  interpreting  data  to  more  fully  understand  the  basic  nature  of 
the  interactions  involved.  Clutter  data  per  se,  is  of  limited  value  and 
should  not  be  the  sole  end  product  of  a  measurements  program. 

1.2  SUMMARY 

In  the  analytic  description  of  phenomena,  simplified  models  are 
generally  used  to  represent  the  interactions  occurring  in  nature.  After  postu¬ 
lating  a  model,  presumably  on  the  basis  of  reasonable  assumptions,  the  analyst 
can  draw  conclusions  concerning  the  outcome  of  a  process  given  a  set  of 
stimuli.  This  same  process  occurs  with  respect  to  clutter.  Clutter,  however, 
involves  an  extremely  complicated  set  of  processes  and  no  simple  model  will 
suffice  to  represent  it.  One  simplified  model  has  assumed  Gaussian  statistics. 
It  is  known  from  the  Central  Limit  Theorem  in  probability  theory  that  the  sum 
of  a  number  of  random  functions,  satisfying  reasonable  constraints,  will  tend 
to  a  Gaussian  random  process  in  the  limit  of  an  infinite  number  of  functions. 
This  theorem,  when  applied  to  clutter,  has  been  used  by  analysts  to  formulate 
a  model  of  clutter  as  a  Gaussian  random  process  for  which  only  the  variance, 
mean  and  spectrum  is  required  in  order  to  specify  the  process.  Results  over  the 
past  years  from  clutter  measurements,  which  have  been  based  primarily  on  this 
model,  can  hardly  be  called  definitive. 

In  Section  2  of  this  report  the  manner  in  which  clutter  affects  the 
division  of  radar  functions  in  an  air  defense  complex  is  described.  The 
important  case  of  surface-to-air  missile  systems  and  the  role  that  clutter 
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plays  in  their  effectiveness  in  defending  against  low-altitude,  high-speed 
aircraft  is  discussed.  The  radar  problems  anticipated  are  used  to  focus 
initial  clutter  measurements,  described  in  the  program  plan,  to  areas  of  key 
importance . 

In  Section  3,  two  generic  forms  of  clutter  models  are  des  ribed. 

In  the  first  model,  called  an  "open-loop"  clutter  model,  the  cluttei  signal 
is  formed  and  added  to  a  signal  representing  the  target;  both  are  cou  'led 
to  the  radar  processor.  This  model,  which  includes  the  Gaussian  clutt  r 
representation,  requires  knowledge  of  the  dependent  variations  in  the  c.  itter 
signal  due  to  variations  in  the  radar  system.  A  second  generic  clutter  n.  del, 
called  a  "closed-loop"  model,  represents  the  received  clutter  and  target 
signals  as  an  operation  on  the  transmitted  signal.  This  model  is  a  closer 
representation  of  the  actual  scattering  process,  which  occurs  in  the  formation 
of  clutter  signals,  than  the  open-loop  model.  The  data  on  required  scattering 
parameters,  with  their  spatial  and  time  dependencies,  is  limited  because  of 
insufficient  suitable  data.  The  program  plan  given  in  Section  10  describes  a 
program  of  clutter  measurements  intended  to  provide  the  required  data  on 
clutter  characteristics. 

In  Sections  4,  5  and  6,  probability  density  functions,  correlation 
functions  and  doppler  spectra  of  clutter  signals  are  discussed  for  the  case 
in  which  only  a  few  major  scatterers  are  present.  The  convergence  of  the 
analytical  results  to  the  case  of  many  scatterers  is  also  discussed.  These 
analytical  studies,  in  conjunction  with  the  closed-loop  clutter  model,  form 
a  basis  for  a  clutter  measurements  program. 

In  Section  7  methods  of  simulating  clutter  are  described:  both 
analog  and  digital  techniques  are  considered.  Emphasis  is  placed  on  advanced 
forms  of  simulators  such  as  would  be  required  for  evaluation  of  future  multi¬ 
dimensional  radar  systems  rather  than  on  additive  filtered  noise  simulation  of 
clutter.  Both  open-loop  and  closed-loop  simulation  are  discussed,  although, 
because  of  its  greater  flexibility,  closed-loop  simulation  is  discussed  in 
greater  detail. 


4 


In  Section  8,  techniques  for  performing  radar  clutter  measurements 
are  described  and  some  of  the  significant  problems,  related  primarily  to 
terrain  and  cultural  observables,  are  discussed.  Clutter  problems  related  to 
innovations  in  radar  technology,  such  as  electronically  steered  phased  arrays, 
are  described  in  Section  9. 

The  program  plan  for  future  clutter  research  which  incorporates  the 
findings  of  the  previous  sections  is  given  in  Section  10.  Included  in  the 
program  plan  are  the  major  areas  of  development  and  research  required  to  perform 
a  meaningful  program  of  clutter  measurements.  It  is  stressed  that:  1)  the 
electromagnetic  observables  should  be  measured  in  such  a  way  that  correlation 
is  maintained  among  the  measured  parameters  and  2)  environmental  descriptors 
are  of  importance  equal  to  that  of  the  radar  sensor  data  and  must  be  associated 
with  it. 
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Section  2 

CLUTTER  PROBLEMS  IN  RADAR  AIR  DEFENSE  SYSTEMS 


Clutter  will  affect  the  division  of  radar  functions  in  an  air  defense 
complex  in  a  manner  dependent  on  the  characteristics  of  the  clutter  and  on  the 
particular  function  (detection  or  tracking).  One  of  the  major  air  defense 
systems  in  use  today  is  the  surface-to-air  missile  system  (SAM)  and  the  role 
that  clutter  plays  in  affecting  air  defense  systems  effectiveness  is  illustrated 
in  the  following  discussion  of  the  problem  SAM  systems  encounter  in  defending 
against  low-altitude,  high-speed  aircraft  attacks.  The  program  plan  given  in 
Section  10  covers  many  facets  of  the  clutter  problem.  Because  the  problem  is 
such  a  large  one,  the  establishment  of  realistic  bounds  is  required  to  focus 
data  acquisition,  at  least  initially,  in  the  critical  areas  of  air  defense. 

This  section  also  provides  a  discussion  of  thesi  critical  areas.  Available 
data  strongly  indicate  that  Line-Of-Sight  (LOS)  coverage  limitations  (caused 
by  terrain  and/or  vegetative  masking)  are  primary  factors  in  determining  the 
effectiveness  of  a  ground  based  (or  truck  mounted)  SAM  system  against  low 

it 

altitude  high  speed  aircraft  targets  in  most  environments.  Appropriate 
guidance  technology  and  hardware  are  well  advanced  and  potential  threats  to 
future  SAM  systems  can  be  expected  to  have  a  very  good  capability  to  exploit 
the  characteristics  of  the  environment  by  proper  planning  and  execution  of 
low  altitude  high  speed  missions.  Future  tactical  aircraft  threats  should  be 
anticipated  with  capabilities  for  flying  in  close  proximity  to  the  ground 
(minimum  clearance  of  200  ft  or  so)  at  speeds  appropriate  for  the  environment. 
For  any  particular  site,  the  interactions  between  anticipated  offensive  and 
defensive  weapons  system  capabilities,  the  tactics  employed,  and  their  effects 
on  target  availability  and  clutter  existence  should  be  evaluated.  Such  a 
study  should  include  defense  site  quality  as  a  function  of  sensor  elevation 
above  ground  level  deployment  density,  and  their  effect  on  expected  overall 


For  this  discussion,  the  optical  LOS  can  be  assumed.  The  effect  of  ridge 
diffraction  and  multipath  to  provide  an  increase  in  the  effective  LOS  is 
probably  small  although  this  premise  is  subject  to  experimental  verification. 
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defense  system  effectiveness.  Knowledge  of  these  factors  would  greatly  assist 
in  quantifying  the  problems  to  be  solved  by  low  altitude  air  defense  systems. 

On  a  qualitative  basis,  consideration  of  the  interactions  discussed 
above  indicate  that  for  SAM  systems  deployed  in  relatively  flat  to  rolling 
terrain,  both  target  availability  and  clutter  existence  (as  a  function  of 
range  from  the  site)  will  be  greatly  different  depending  on  whether  or  not 
the  SAM  sensor  is  grourd  based  (vehicle  mounted)  or  elevated  above  objects 
contributing  to  the  local  masking  (such  as  small  terrain  variations,  trees, 
buildings  etc.). 

The  effects  of  local  masking  may  cause  significantly  decreased  LOS 
ranges  to  both  target  and  clutter  sources.  It  is  expected  that  non-elevated 
sensors  deployed  in  temperate  zones  would  often  be  subject  to  initial  LOS  to 
target  at  ranges  of  1-3  n.mi.  As  potential  clutter  sources  are  masked  beyond 
the  range  of  the  object  causing  the  local  mask,  the  target  should  be  available 
for  detection  in  a  clutter-free  state  in  many  situations.  Because  of  the  short 
LOS  ranges  involved,  such  SAM  systems  would  be  required  to  have  very  short 
reaction  times  (approximately  6-8  seconds)  to  be  effective. 

SAM  systems  employing  a  sensor  elevated  above  the  local  mask,  on  the 
other  hand,  are  presented  with  significantly  different  situations  when  deployed 
in  the  same  environment.  Target  LOS  ranges  on  the  order  of  10-15  miles  would 
be  common  and  clutter  might  exist  over  a  significant  portion  of  the  target 
path.  While  the  target  might  be  available  for  detection  in  a  clutter  free 
situation,  both  clutter  and  multipath  effects  could  be  expected  to  be  important 
for  a  large  portion  of  the  target  track.  From  the  SAM  sensor  design  standpoint, 
the  clutter  problem  is  more  difficult  for  the  elevated  sensor  because  of  the 
longer  ranges  required.  The  availability  of  the  target  at  long  ranges  does, 
however,  greatly  increase  the  allowable  SAM  reaction  time. 

Experience  has  shown  that  accurate  determination  of  LOS  to  target 

for  ground  based  defenses  requires  optical  measurements  (or  their  equivalent 

(2) 

from  stereo  pair  photos)  of  the  objects  causing  local  masking;  It  is  expected 
that  target  LOS  determined  from  topographical  map  data  would  be  adequate  for 
cases  where  the  sensor  was  elevated  well  above  the  local  mask. 
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Accurate  determination  of  LOS  to  clutter  sources  is  predictable  only 
when  significant  mask  angles  are  involved  (valleys,  etc.)*  It  is  recommended 
that  field  survey's  using  a  helicopter  borne  clutter  mapping  radar  (PPI)  be 
used  in  evaluating  the  LOS  to  clutter  sources  for  elevated  sensors. 

It  is  assumed  that  maps  of  radar  LOS  as  a  function  of  target  altitude 
above  the  local  terrain  can  be  obtained.  These  maps  will  govern  the  maximum 
allowable  interval  between  successive  searches  in  various  directions  for 
defense  against  radially  entering  targets.  That  is,  regions  where  small  LOS 
ranges  to  target  exist  will  be  searched  frequently  in  order  to  obtain  detection 
within  the  SAM  reaction  time. 

Provision  of  adequate  interval  of  LOS  exposure  by  reducing  vegetation 
masking  effects  does  not,  alone,  necessarily  allow  the  defensive  system  to  be 
effective  because  its  radar  may  not  be  able  to  acquire  and  track  the  low- 

altitude  target  which  is  in  close  proximity  to  the  terrain.  The  terrain 

produces  clutter  which  typically  affects  the  radar  system  detection  capability. 
In  general,  the  radar  cross  section  of  an  attacking  tactical  aircraft  may  be 
on  the  order  of  10  to  20  dBsm  with  potential  reduced  radar  cross  sections  of 

future  aircraft  being  lower  than  0  dBsm.  On  the  other  hand,  the  radar  cross 

section  of  terrain  and  cultural  features  per  square  meter  of  surface  area  may 
vary  from  -30  dB  to  perhaps  +30  dB,  with  -15  dB  being  a  typically  encountered 
average  value.  For  practical  SAM  radar  systems,  the  area  of  the  pulse  packet 
on  the  terrain  results  in  a  backscattered  radar  signal  from  the  terrain  which 
is  typically  very  large  compared  to  the  return  from  the  target  (assuming  the 
terrain  in  the  vicinity  of  the  target  is  not  masked),  especially  for  long- 
range  detection  of  low-altitude  tactical  aircraft  attacking  radially  over 
"choppy"  terrain.  Using  a  clutter  LOS  map,  those  areas  where  terrain  is  masked 
can  be  determined  and  detection  of  visible  targets  m  these  regions  will  be 
relatively  easy  due  to  the  reduced  amplitude  of  the  clutter.  A  clutter  signal 
will  be  present  even  if  the  terrain  in  the  same  range  gate  as  the  target  is 
masked,  because  of  unmasked  terrain  in  the  antenna  sidelobes  or  in  a  part  of 
the  main  lobe.  The  clutter  amplitude  will,  however,  be  significantly  reduced 
if  clutter  in  the  main  lobe  is  masked.  The  target  visibility  will  be  governed 
by  the  target  altitude  as  described  by  the  radar  LOS  maps  previously  discussed. 
Thus,  search  of  an  area  where  clutter  is  masked  will  provide  high  detectability 
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of  targets  having  radar  LOS.  The  availability  of  LOS  to  target  must  be 
evaluated,  however,  in  terms  of  an  assessment  of  minimum  altitudes  which  can 
be  achieved  by  potential  threats. 

To  provide  the  required  target  detection  (and  tracking)  capability 
in  regions  where  the  terrain  is  not  masked,  it  is  necessary  to  make  use  of 
other  characteristics  in  addition  to  the  amplitude  of  the  radar  return  signal . 
As  the  target  typically  has  a  relatively  high  radial  velocity  (if  it  is  a 
threat  to  the  SAM  system) ,  the  Doppler  shift  of  the  return  from  the  target  is 
typically  used  to  provide  an  additional  "dimension"  which  aids  in  the 
separation  of  the  target  from  clutter.  Experience  has  shown,  however,  that  in 
many  environments,  the  Doppler  separation,  azimuth  beamwidth  and  range 
resolution  capability  are  not  sufficient  to  permit  effective  detection  of 
low-altitude  targets  at  long  ranges.  This  situation  requires  (1)  that  the 
reaction  times  of  the  overall  defense  system  be  reduced  so  that  shorter 
ranges  of  detection  are  permitted  without  impairing  SAM  system  effectiveness, 
(2)  that  the  size  of  the  radar  resolution  cell  be  reduced,  or  (3)  that 
additional  characteristics  of  the  radar  return  be  employed  which  improve 
separation  of  the  target  return  from  the  clutter  signals.  Any  characteristic 
of  the  target  return  which  is  significantly  different  from  the  clutter  return 
alone,  whether  on  an  instantaneous  or  a  statistical  basis,  whether  natural 
or  induced,  potentially  offers  a  means  for  improving  the  capability  of  future 
radar  systems  to  detect  targets  in  clutter.  These  additional  characteristics 
include  such  things  as  frequency  sensitivity,  phase  coherence,  fluctuation  and 
scintillation  rates,  target  size,  polarization,  angle-of-arrival  variations 
and  knowledge  of  the  masking  and  clutter  characteristics  at  a  particular  site. 

A  much  improved  knowledge  of  the  characteristics  of  both  clutter  and  target - 
clutter  interaction  effects  is  needed  to  permit  efficient  design  and  use  of 
multidimensional  radar  sensing  for  the  development  of  more  effective  SAM 
systems  and  to  assess  their  performance  capabilities.  These  data  are  needed 
for  all  environments  in  which  the  SAM  system  will  be  required  to  operate.  A 
program  plan  for  experimentally  obtaining  the  required  clutter  data  and  methods 
for  data  reduction/correlation  and  use  in  simulation  is  given  in  Section  10. 
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For  the  short  and  long  range  conditions  mentioned  above,  the  clutter 

signal  amplitude  will  probably  cover  a  60  dB  dynamic  range.  During  a  clutter 

measurements  program,  as  described  in  Section  10,  it  would  be  desirable  to  use 

a  sensor  radar  providing  clutter  cross  sections  which  approach  typical  target 

2 

radar  cross  sections  (1  m  ).  A  synthetic  aperture  system  with  an  effective 
resolution  of  1  milliradian  and  a  pulse  length  of  20  nanoseconds  will,  for 

2 

much  of  the  terrain,  provide  effective  clutter  cross  section  approaching  1  m 
(assuming  a  nominal  20  dB  terrain  reflectivity) .  Variation  in  the  transmitted 
pulse  length  by  a  factor  of  100-500  (through  the  use  of  pulse  compression 
techniques)  should  provide  a  good  representative  distribution  of  effective 
clutter  cross  sections  and  associated  target -to-clutter  ratios. 

The  potential  existence  of  multiple  transmission  paths  to  a  low- 

altitude  target  (or  any  target  at  low  elevation  angles)  can  impair  the  target 

detection  range  by  causing  fading  of  the  target  return  and  widening  of  the 

apparent  Doppler  spectrum  of  the  target.  In  addition,  multipath  can  also 

seriously  degrade  the  accuracy  of  target  tracking  radars  whenever  the  target 

is  separated  from  the  clutter  by  less  than  an  antenna  beamwidth  (which  is  the 

typical  case  for  low-altitude  targets).  In  particular,  the  multipath  signals, 

when  added  vectorial ly  to  the  direct  radar  returns,  cause  variations  in  the 

apparent  angle  of  arrival,  and,  depending  on  the  particular  situation,  the 

apparent  angle  of  arrival  may  appear  to  come  from  angular  regions  outside  the 

angular  extent  of  the  target  and  the  clutter  for  a  significant  percentage  of 
(3) 

the  time.  These  errors,  together  with  signal  fluctuations  and  fading,  greatly 
complicate  the  tracking  of  low-altitude  targets.  It  should  be  noted  that  the  • 
multipath  returns  from  low-angle  targets  typically  exist  near  the  Doppler-shifted 
frequency  of  the  target  and,  therefore,  the  multipath  effects  cannot  be  rejected 
on  the  basis  of  Doppler  separation  alone.  It  is  expected  that  the  use  of 
additional  characteristics  such  as  previously  described  will  be  required  to 
allow  the  desired  improvements  in  tracking  of  targets  in  clutter.  Study  of 
multipath  effects  is  included  in  the  program  plan  of  Section  10. 

A  similar  situation  prevails  for  low-level  missile  systems.  Here 
an  illuminator  is  used  to  reflect  energy  from  the  incoming  target.  The 
missile  "sees"  this  energy,  and  tracking  circuits  within  the  missile  force 
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it  into  collison  course  with  the  target.  However,  for  low  flying  vehicles, 
the  illuminator  also  illuminates  ground  clutter  so  that  the  resultant 
illumination  pattern  is  modulated  as  a  function  of  target  position.  Simul¬ 
taneously,  the  signal  reflected  from  the  target  also  illuminates  ground  clutter, 
forming  a  diffuse  "image."  The  combination  of  both  multipath  perturbations 
place  a  lower  bound  on  system  effectiveness,  preventing  system  deployment 
against  targets  penetrating  below  a  given  altitude.  However,  lack  of 
detailed  information  on  clutter  multipath  prevents  an  accurate  assessment  of 
the  lower  bound,  short  of  actual  firings. 

The  foregoing  discussion  has  briefly  indicated  the  complexity  of  the 
SAM  system  design  problem  and  the  importance  of  target  and  clutter  characteristics 
knowledge  in  arriving  at  a  suitable  SAM  design.  Although  the  masking  problem 
could  often  be  solved  by  elevating  the  site  above  the  local  mask  (especially 
in  nearly  flat  terrains),  clutter  effects  would  still  restrict  the  effectiveness 
of  the  SAM  systems  in  many  environments. 
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Section  3 
CLUTTER  MODELS 


A  clutter  model,  for  the  purposes  of  the  discussion  below,  is  a 
methodology  to  determine  the  clutter  and  clutter/target  interactions  in  a 
radar  system.  The  form  of  a  clutter  model  can,  in  part,  suggest  methods  and 
techniques  for  clutter  simulation  and  in  some  cases  a  model  and  a  simulation 
technique  have  become  synonymous.  The  representation  of  clutter  return  as  a 
narrow-band  Gaussian  process  and  use  of  an  additive,  filtered  Johnson  noise 
source  for  clutter  simulation  purposes  are  usually  referred  to  interchangeably. 
This  representation  can  lead  to  erroneous  results  because  the  clutter  may  be 
only  partially  simulated.  For  example,  the  spectrum  of  a  noise  source  can  be 
adjusted  so  that  the  simulated  clutter  has  the  appropriate  correlation  interval 
as  a  function  of  range  (time  delay)  on  a  single  radar  scan  but  then  the  cor¬ 
relation  between  clutter  at  the  same  range  but  on  successive  radar  scans  may 
be  incorrect  (see  Section  7.2  for  a  further  discussion).  The  model  itself 
may  be  incorrect;  i.e.,  the  signal  from  a  Gaussian  process  is  unbounded, 
whereas  the  clutter  signal  amplitude  is  necessarily  bounded  and  thus  not 
Gaussian.  In  the  following  section  the  major  applications  of  clutter  models 
and  two  suggested  generic  forms  of  clutter  models  are  described. 

3.1  CLUTTER  MODEL  APPLICATIONS 


It  is  useful  to  consider  potential  applications  of  a  clutter  model 
in  order  to  evaluate  the  degree  of  complexity  required  in  the  model  and  in  poten¬ 
tial  simulation  techniques. 


In  the  past,  the  most  frequent  use  of  clutter  models  has  been  to 
determine  average  radar  performance  criteria,  such  as  the  probabilities  of 
false  alarm  and  detection  and  RMS  track  angle  error,  for  various  types  of 


terrain  and  radar  system  parameters.  Most  of  the  past  clutter  measurements 
were  directed  to  this  application  and  data  were  obtained  on  average  clutter 
backscatter  cross  section,  probability  distributions,  and  clutter  spectra  for 
/arious  types  of  terrain  and  radar  system  configurations. 

3.1.2  Engagement  System 

In  an  engagement  system  application  of  a  clutter  model,  the  time 
history  of  system  performance  must  be  evaluated  in  a  dynamic  situation  in  which 
the  target/clutter  radar  information  is  a  system  input  on  which  future 
action  is  to  be  based.  An  example  of  this  type  of  application  is  a  missile-borne 
radar  seeker  for  which  it  may  be  desirable  to  determine  system  performance  for 
general  classes  of  terrain  or  in  some  cases  for  a  very  specific  site. 

The  requirements  on  the  clutter  model  imposed  by  an  engagement  system 
application  are  generally  more  severe  than  those  imposed  by  an  average 
performance  evaluation.  For  example,  in  an  engagement  system  one  may  be  inter¬ 
ested  in  the  peak  clutter  amplitude  in  studies  of  loss  of  angle  or  range  track. 
The  peak  clutter  amplitude,  achieved  during  only  a  small  portion  of  the  total 
period  of  concern,  can  vary  considerably  without  significantly  affecting  the 
average  clutter  amplitude  or  average  system  performance. 

3.1.3  Comparative  Evaluation  of  Radar  Systems 

It  may  be  necessary  to  critically  compare  radar  performance  for 
various  systems  having  different  design  configurations  (frequency  band,  pulse 
shape,  antenna  patterns  and  polarization,  data  processing)  in  a  prescribed 
clutter  and  target  environment.  In  this  case  it  is  necessary  for  the  clutter 
model  to  have  sufficient  flexibility  that  physical  differences  in  illuminator 
and/or  observer  radar  parameters  can  be  represented  to  good  accuracy  by 
corresponding  changes  in  the  characteristics  of  the  clutter  and  target  signals. 
Included  in  this  application  is  the  use  of  the  model  for  radar  system  design 
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wherein  the  performance  capabilities  of  specific  radar  techniques,  such  as 
electronic  beam  scanning,  MTI,  Doppler  filtering,  etc.,  could  be  optimized  on 
paper  or  in  laboratory  simulations  without  requiring  recourse  to  field 
experiments. 

3.2  GENERIC  CLUTTER  MODELS 

In  this  section  are  considered  the  basic  forms  a  radar-clutter  model 
might  take.  Classes  of  models  are  presented  rather  than  final,  detailed  models 
factors  that  must  be  considered  in  a  sophisticated  clutter  model  are  described 
here  and  discussed  in  greater  detail  in  Section  7.  It  is  assumed  that  only  a 
relatively  refined  clutter  representation  is  of  interest;  very  little  attention 
is  therefore  paid  to  the  use  of  additive  band-limited  noise  to  represent  a 
complete  clutter  signal,  even  though  such  simulation  is  adequate  for  some 
applications. 

It  seems  reasonable  to  consider  two  general  classes  of  models;  any 
model  must  lie  in  one  of  these  classes  or  else  must  contain  a  combination  of 
elements  belonging  to  both  classes.  These  classes  are  henceforth  called 
"open-loop"  models  and  "closed-loop"  models. 

3.2.1  Open -Loop  Model 

In  Figure  3-1  is  shown  the  basic  form  of  an  open-loop  clutter  model. 
Here  a  Clutter-Signal  Generator  produces  a  signal  with  characteristics  similar 
to  the  one  received  from  the  clutter  for  the  specific  radar  being  evaluated. 
Similarly,  the  Target-Signal  Generator  produces  a  signal  with  characteristics 
similar  to  that  received  from  the  target  (note  that,  although  only  clutter 


OUTPUT 


Figure  3-1  OPEN-LOOP  CLUTTER  MODEL 
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models  are  considered,  the  technique  for  incorporating  the  target  signal 
characteristics  cannot  be  ignored).  These  two  signals  are  then  added  and  the 
sum  fed  to  the  radar  processor  and  the  desired  measures  of  performance  are 
computed. 


3.2.2  Closed-Loop  Model 

In  Figure  3-2  is  shown  the  basic  form  of  a  closed-loop  clutter  model. 
Here  a  signal  corresponding  to  that  transmitted  by  the  radar  is  processed  to 
obtain  signals  having  characteristics  of  the  clutter  and  target  signals. 
Provision  is  made  for  the  mutual  interaction  of  target  and  clutter  to  modify 
both  signals,  as  shown  by  the  coupling  between  clutter  and  target-signal  form¬ 
ation  processes.  The  resulting  signals  are  added  and  fed  to  the  radar  pro¬ 
cessor.  The  indirect-generation  method  is  therefore  more  nearly  like  the 
actual  radar  situation  (in  which  the  original  signal  is  modified  by  reflection 
from  the  target  and  clutter  scatterers)  than  is  the  direct-generation  method. 


Figure  3-2  CLOSED-LOOP  CLUTTER  MODEL 

With  the  open-loop  model  the  characteristics  of  the  clutter  signal, 
such  as  probability  density  functions,  spectrum,  polarization  characteristics 
and  correlation  functions  among  radar  observables,  must  be  known  as  a  function 
of  the  radar  parameters  and  terrain  type.  Most  of  the  previous  work  in  the 
clutter  field  has  resulted  in  use  of  some  form  of  "open-loop"  model  for  radar 
system  analysis.  For  example,  it  is  conventional  to  assume  a  Gaussian  probability 
distribution  for  the  clutter  with  average  power  given  by  the  integral  of  the 
terrain  reflectivity,  two-way  antenna  pattern  and  pulse-shape  product  over  the 
terrain  surface.  The  clutter  correlation  function  as  seen  by  the  radar  is 
determined  from  the  internal  "motion"  clutter  spectrum,  typically  assumed  to 
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have  a  Gaussian  shape,  and  from  the  effects  of  antenna  pattern  and  platform 
motion,  range  gate  movement,  and  various  system  instabilities. 

In  the  closed-loop  clutter  model,  the  scattering  processes  which 
result  in  the  clutter  return  are  emphasized.  Most  clutter  analyses  usually 
begin  with  some  form  of  closed-loop  model.  By  assuming  the  terrain  satisfies 
a  set  of  restrictions  regarding  spatial  stationarity  (in  a  statistical  sense), 
roughness  scale,  and  roughness  curvature,  one  obtains  a  tractable  analytical 
problem.  Terrain  of  concern  in  radar  systems,  however,  does  not  satisfy  the 
required  set  of  assumptions  at  (for  example,  roughness  must  usually  be  assumed 
either  very  large  or  very  small  with  respect  to  wavelength)  typical  radar 
frequencies.  Thus,  phenomena  such  as  depolarization,  local  masking,  and 
multiple  scattering  cannot,  at  the  present  time,  be  included  analytically. 

Hiis  does  not  mean  that  use  of  the  closed-loop  model  will  completely  represent 
the  physical  processes  mentioned  above,  but  rather  that  consideration  of  clutter 
as  a  scattering  process  in  a  closed-loop  model  can  serve  as  a  guideline  for  the 
planning  and  data  interpretation  of  future  clutter  measurements  programs. 


Section  4 

PROBABILITY  DENSITY  FUNCTIONS  FOR  A  SMALL  NUMBER  OF  SCATTERERS 


4.1  GENERAL 

If  ground  clutter  is  to  be  realistically  represented,  it  is  neces¬ 
sary  that  the  clutter  signal  have  statistical  properties  closely  related  to, 
if  not  identical  with,  those  of  the  signal  that  would  be  received  from  the 
radar  system.  Perfert  similarity  is  not  likely  to  be  necessary,  particularly 
in  view  of  the  wide  range  of  clutter  types  that  can  exist.  Two  of  the 
characteristics  that  must  be  represented  with  reasonable  accuracy  are  the  power 
level  and  amplitude  distribution.  The  clutter  correlation  function  also  must 
be  determined  in  some  cases,  especially  if  the  radar  to  be  evaluated  uses 
integration  or  clutter  cancellation  techniques  to  achieve  MTI .  Spectral 
characteristics  will  also  be  important  if  a  coherent  radar  is  to  be  considered. 

In  this  section,  probability  density  functions  of  the  in-phase,  and 
also  quadrature  phase,  components  of  the  clutter  signal  and  of  the  envelope 
of  the  clutter  signal  are  considered  with  emphasis  on  the  configuration  involving 
a  small  number  of  scatterers.  For  phase-coherent  detection,  one  is  interested 
in  the  probability  density  function  of  the  in-phase  and  quadrature  signals, 
whereas  for  envelope  detection,  the  envelope  probability  density  function  is 
required.  The  few-scatterer  configuration  is  significant  because  the  statistics 
of  the  clutter  signal  will  vary  widely  when  only  a  few  scatterers  are  dominant. 

It.  is  important  in  experimental  investigations,  such  as  those  described  in  the 
program  plan  of  Section  10,  to  determine  whether  this  effect  is  likely  to  be 
significant. 

The  most  common  clutter  representation  is  based  on  application  of 
the  Central  Limit  Theorem  to  show  that  the  in-phase  and  quadrature  components 
of  the  clutter  signal  are  approximately  independent  Gaussian  random  functions 
and  that  the  envelope  of  the  clutter  signal  will  therefore  be  approximately 
Rayleigh  distributed.  This  assumption  is  valid  for  many  cases  of  interest. 
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For  example,  scattering  from  most  meteorological  targets  (raindrops,  snowflakes) 
or  from  a  relatively  uniform  rough  surface  will  produce  a  signal  envelope  having 
a  Rayleigh  amplitude  distribution. 

However,  many  forms  of  terrain  and  man-made  clutter  will  not  satisfy 
the  requirements  for  application  of  the  Center  Limit  Theorem.  Whether  they  do 
depends  on  the  radar  resolution-cell  size.  As  the  radar  resolution  increases 
(size  of  resolution  cell  decreases)  it  is  expected  that  the  clutter  signal  will 
deviate  to  a  larger  degree  from  a  Rayleigh  amplitude  distribution  because  of  the 
reduced  number  of  effective  scatterers.  Two  of  the  requirements  that  must  be 
met,  for  application  of  the  Central  Limit  Theorem  are^  :  1)  there  must  be  many 
scatterers  contributing  to  the  clutter  signal,  and  2)  the  mean-square  Value  of 
the  field  from  any  one  of  these  scatterers  must  be  much  less  than  the  sum  of  the 
mean-square  values  of  the  fields  from  all  of  the  scatterers.  Although  these 
facts  are  well  known,  there  seems  to  be  little  information  available  relative 
to  how  many  scatterers  are  required  or  how  much  smaller  than  the  total  mean- 
square  value  any  individual  mean-square  value  must  be.  Asymptotic  expressions 
indicate  that  for  more  than  five  equal  amplitude  scatterers  the  Gaussian  and 
Rayleigh  statistics  are  reasonably  accurate.  However,  the  total  failure  of 
the  asymptotic  expression  when  there  are  two  equal  scatterers  raises  questions 
concerning  its  accuracy  for  five  to  ten  scatterers.  For  this  reason,  the 
probability  densities  for  in-puase  and  quadrature  components,  as  well  as  those 
for  the  envelope,  of  the  clutter  signal  received  from  a  two-scatterer  model 
and  from  a  four-scatterer  model  have  been  computed  for  several  specific  values 
of  relative  scatterer  amplitude.  These  probability  densities  are  compared 
with  those  computed  from  various  commonly  used  approximations  (depending  on 
the  situation)  to  determine  how  accurate  the  approximations  are  when  there 
are  a  small  number  of  scatterers. 

In  all  of  the  work  reported  here  it  is  assumed  that  the  relative 
phase  of  the  signal  received  from  each  scatterer  is  uniformly  distributed  over 
2 rr radians  and  that  the  relative  phases  of  the  scatterers  are  statistically 
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independent.  This  assumption  is  valid  if  the  placement  of  scatterers  within 
the  resolution  cell  of  the  radar  is  assumed  random  and  the  random  process 
giving  rise  to  the  scattered  signal  is  considered  as  the  ensemble  of  possible 
scatterer  placements.  For  a  given  set  of  scatterers  in  a  resolution  cell  the 
assumption  will  also  be  satisfied  if  the  relative  phase  among  scatterers  is 
independent  and  changes  (due  to  scatterer  movement,  internal  scatterer 
variations  or  changes  of  radar  frequency)  with  a  variance  much  larger  than 
radians . 


It  is  possible  that  this  assumption  will  not  be  satisfied  if  the 
clutter  at  a  given  radar  frequency  from  only  a  single  resolution  cell  is 
considered  and  if  the  relative  scatterer  movements  are  small  relative  to  a 
wavelength.  Probability  density  functions  can  be  derived  for  this  case  if 
there  are  many  scatterers  present.  This  configuration  can  be  included  in 
the  closed-loop  clutter  model  (Section  7)  and  its  presence  in  experimental 
data  may  be  detectable  from  the  clutter  correlation  functions  versus  time 
and  frequency  (Section  5) . 

4.2  PROBABILITY  DENSITY  FUNCTIONS 


The  clutter  signal  received  from  V  scatterers  can  be  written 
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where  the  and  ^  are  the  relative  amplitude  and  phase  of  the  signal  from 
the  ^  scatterer  and  /?  is  the  envelope.  The  sum  signal  S  can  now  be  exp- 
pressed  in  terms  of  an  in-phase,  *,  and  quadrature  component,^,  given  by 
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For  many  radar  and  clutter  situations,  ^  can  be  assumed  to  be  uniformly 
distributed  as  previously  discussed.  The  can  be  considered  to  be  fixed 
constants  (as  they  are  here)  or  to  be  random  variables  having  identical  distri¬ 
butions  (as  they  are  in  Reference  5).  In  Appendix  t  t'ir‘  probability  density 
functions  for  the  in-phase  and  quadrature  signals  {  f  ft)  and  /of  y)  )  and  for  the 
envelope,  are  determined  by  numerical  integration.  These  probability 

functions  are  exact  to  within  the  numerical  approximations  used  in  performing 
the  computations.  The  probability  density  />  ft)  appropriate  for  phase  coherent 
detection  analysis  is  compared  in  Appendix  A  to  two  frequently  used  approxi¬ 
mations  to  />(t)  for  the  cases  of  two  and  four  )  independent 

scatterers  and  for  several  specific  values  of  relative  scatterer  amplitude. 

The  approximate  functions  are  the  Gaussian  density  function  f^f*)  and  an 
asymptotic^  expansion  y>  {*)  which,  as M-eo,  converges  to  a  Gaussian  function. 


The  calculated  envelope  probability  density, f(R),  is  also  compared 
in  Appendix  A  to  the  Rayleigh  density,  fJR)  ,  to  an  asymptotic  expansion 
p  (R)  which  converges  to  a  Rayleigh  density^  1  for  largeV,  and  to  the 

'T't 


...  (6) 


Nakagamiv'''  m-distribution  function  (X)  .  Graphs  of  these  probability  functions 
and  the  analytical  expressions  used  can  be  found  in  Appendix  A. 


4.3  DISCUSSION  OF  RESULTS 


Computations  of  probability  densities  for  in-phase  (or  quadrature) 
components  and  for  amplitudes  of  two  or  four  randomly  phased  scatterers  have 
been  made  and  the  results  are  presented  in  Appendix  A.  For  the  in-phase 
(or  quadrature)  component  the  probability  density ,  f>(t) ,  is  non-Gaussian  for 
two  vectors  regardless  of  their  relative  amplitudes.  This  is  also  true  for 
four  scatterers  of  unequal  amplitudes.  For  four  scatterers  of  equal,  or  nearly 
equal,  amplitudes  the  probability  density  is  very  nearly  Gaussian.  The 
asymptotic  formulation  for  the  Central  Limit  Theorem,  as  given  in  Equation  A-10 
of  Appendix  A,  is  not  much  more  accurate  than  the  unmodified  Gaussian  probability 
density  function  for  four  scatterers.  This  fact  makes  questionable  the  generally 
used  assumption  that  the  Center  Limit  Theorem  is  very  nearly  valid  for 
five  or  more  random  variables. 
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Computations  of  the  probability  density  of  the  envelope,  fW), 
indicates  that  the  Rayleigh  approximation  is  fairly  good  for  four  equal  scat- 
terers  but  becomes  poor  if  the  scatterers  pre  unequal  in  amplitude.  Again  the 
asymptotic  formulation  does  not  seem  particularly  good.  The  Nakagami  m- 
distribution  is  a  very  good  one  when  there  is  one  strongly  dominant  scatterer 
with  three  much  smaller  ones,  and  it  appears  to  be  a  better  approximation  for 
all  of  the  four-scatterer  cases  considered  here.  For  two  scatterers,  neither 
the  Rayleigh  nor  the  Nakagami  m-distribution  is  very  good,  but  the  latter  is 
the  more  nearly  accurate,  especially  when  one  scatterer  is  much  larger  than 
the  other. 


f7) 

For  six  equal  scatterers,  results  by  Greenwood  and  Durand v  have 
been  used  to  indicate  that  all  of  the  approximations  are  reasonably  valid, 
(/})  being  more  accurate  than  (4)  ,  which  in  turn  is  more  accurate  than 


From  the  trends  found  for  the  cases  considered  in  this  investigation, 
some  conclusions  regarding  the  envelope  probability  density  function,  ftf), 
can  be  drawn: 

1.  For  two  scatterers,  as  would  be  expected,  f  (*i)  is  not  well 
represented  by  Rayleigh  or  Nakagami  m-distribution  functions; 
probably  the  latter  is  a  better  approximation. 

2.  For  four  scatterers,  />W)is  quite  well  represented  by  the 
Nakagami  m-distribution  for  a  wide  range  of  scatterer  amplitudes, 
whereas  the  Rayleigh  distribution  is  a  reasonable  approximation 
only  when  the  scatterers  are  equal,  or  very  nearly  equal,  in 
amplitude. 

The  configuration  of  a  relatively  few  significant  scatterers  in  the 
radar  resolution  cell  has  been  emphasized  because  this  case  seems  likely  to 
arise  for  some  types  of  terrain  with  the  higher  resolution  being  used  today  and 
because  this  configuration  and  its  effects  on  clutter  measurements  and  data 
interpretation  has  been  largely  ignored  in  previous  measurements.  The  many 
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scatterer  case  is  justified  for  certain  types  of  terrain  and  examples  of  both 

the  many- scatterer^  and  few-scatterer^  cases  can  be  found  in  clutter 

measurements  data.  High  resolution  C-band  measurements  (0.1  sec)  in  conjunction 

with  site  surveys  performed  in  England^^  indicated  major  scatterers  (cross 

section  *  100  ^  )  were  man-made  metal  structures  occurring  with  a  density  of 
_  2 

20/mi  .  The  maximum  cross  section  measured  with  lower  resolution  (4  /*,  sec 
pulse)  was  approximately  100  *r2with  a  50  percentile  of  10  •mt.  For  this  particular 
example  the  few  major  scatterer  configuration  would  probably  predominate  in  de¬ 
termining  the  clutter  characteristics  for  a  large  portion  of  the  measured  data. 

Because  of  the  large  effect  the  occurrence  of  the  few  major  scatterer 
case  has  on  clutter  probability  density  functions  and  correlation  functions 
(discussed  in  Section  5)  and  the  implications  with  respect  to  clutter  simulation 
it  is  important  to  be  able  to  isolate,  identify  and  determine  the  scattering 
characteristics  of  major  scatterers  during  the  performance  of  clutter  measure¬ 
ments.  The  use  of  a  high  resolution  coherent  radar  is  probably  necessary. 


Section  5 

CLUTTER  CORRELATION  FUNCTIONS 


5 . 1  GENERAL 

Correlation  functions  are  of  interest  in  the  analytical  and  experi¬ 
mental  investigation  of  clutter  and  its  effects  on  radar  systems  because  they 
are  a  means  of  describing  the  dependence  of  one  random  variable  on  another. 

In  the  design  of  radar  systems,  correlation  functions  frequently  are  introduced 
when  performance  criteria,  such  as  clutter  cancellation  in  an  MTI  system,  are 
evaluated  using  an  average  mean  square  analysis . ^ ‘ ^  In  the  analysis  of  clutter 
measurement  data,  dependencies  of  one  parameter  on  another,  such  as  between 
the  orthogonal  polarization  components  of  the  clutter  signal,  can  be  deter¬ 
mined  by  evaluating  their  correlation  functions.  In  clutter  simulation  it  is 
necessary  to  ensure  the  correct  form  of  dependence  (correlation)  among  the 
simulated  clutter  parameters  for  accurate  evaluation  of  radar  systems,  especially 
wnen  clutter  rejection  techniques  such  as  MTI,  Doppler,  or  polarization  ratio 
cancellation  are  employed.  Thus  correlation  functions  play  a  major  role  in  all 
aspects  of  the  clutter  problem. 

The  discussion  in  this  section  is  directed  to  the  correlation 
function  relating  signal  intensity  to  time  delay,  frequency  variation  and 
resolution  cell  movement.  Signal  intensity  is  only  one  of  the  clutter 
observables  (see  Section  8)  aid,  for  radar  system  evaluation,  one  is  interested 
in  a  much  broader  class  of  correlation  functions.  Parameters  such  as  the  cor¬ 
relation  between  polarization  components  of  the  clutter  signal,  between 
amplitude  and  p.iase  and  between  angle  of  arrival  at  different  spatial  locations 
are  important  characteristics  of  clutter.  Many  of  these  clutter  characteristics 
are  presently  not  subject  to  analytical  determination  and  a  clutter  measurements 
program  as  described  in  Section  10,  embodying  an  approach  as  outlined  in 
Section  8  is  required. 
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In  this  section  the  correlation  function  for  the  signal  intensity 
(squared  envelope)  of  the  clutter  signal  is  derived  for  an  arbitrary  number 
of  scatterers  in  the  resolution  cell  of  the  radar.  The  probability  density 
function  of  the  clutter-signal  amplitude  was  discussed  in  Section  4.  This 
information  is  basic  to  any  statistical  characterization  of  the  clutter  signal: 
any  simulator  to  be  used  for  radar  evaluation  should  be  capable  of  matching 
reasonably  well  the  amplitude  characteristics  of  the  clutter  signal.  Unspeci¬ 
fied  as  yet,  however,  and  not  considered  in  the  previous  section,  is  the  rate 
at  which  the  clutter  signal  changes  with  time  due  to  motion  of  the  scatterers, 
and  the  change  in  the  clutter  signal  if  either  the  frequency  of  the  radar  or  the 
observed  region  (set  of  scatterers)  is  changed.  If  a  radar  employs  some  form 
of  clutter  cancellation  technique,  such  as  MTI,  then  the  performance  of  the 
radar  system  will  vary  depending  on  the  amount  that  the  clutter  signals  are 
correlated  pulse-to-pulse,  scan-to-scan,  or  intrapulse.  Consequently,  if 
the  radar  is  to  be  evaluated  using  a  simulator,  the  simulator  must  produce  a 
clutter  signal  having  appropriate  time,  frequency,  and  spatial  correlation 
(although  in  a  specific  case  only  some  of  these  properties  might  be  important) . 

5.2  INTENSITY  CORRELATION  FUNCTION 

Analyses  of  correlation  properties  of  signals  arising  from  large 
numbers  of  scatterers  have  been  made  by  numerous  authors;  some  of  these 
analyses  will  be  referred  to  subsequently.  In  this  section  are  considered 
the  correlation  properties  of  clutter  signals  obtained  from  smaller  numbers 
of  scatterers,  specifically  for  cases  for  which  the  in-phase  and  quadrature 
components  of  the  received  signal  are  not  Gaussian  (i.e.,  when  the  conditions 
of  the  Central  Limit  Theorem  are  not  satisfied  by  the  random  process) .  The 
relationships  between  these  results  and  the  results  obtained  in  the  many- 
scatterer  cases  are  examined  as  well. 

It  is  shown  in  Appendix  B  that  the  correlation  function  for  signal 
intensity  (squared  envelope)  can,  in  general,  be  expressed  as  the  product  of 
three  correlation  functions,  one  for  frequency  shift,  one  for  time  (scatterer 
motion  effect),  and  one  for  displacement  of  the  radar  resolution  cell.  It 
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is  not  possible  at  the  present  time  to  verify  experimentally  the  complete 
factorization  of  the  clutter  correlation  function  as  described  above.  The 
results  of  this  section  are  therefore  intended  to  demonstrate  the  method  of 
factorization  of  the  correlation  function,  with  delineation  of  the  assumptions 
made  regarding  the  characteristics  of  the  scatterers,  for  the  case  of  an 
arbitrary  number  of  scatterers  in  the  radar  resolution  cell.  Discussion  of 
measurements  techniques  applicable  to  determine  the  correlation  function  are 
given  in  the  program  plan,  Section  10,  and  the  effects  of  clutter  correlation 
in  simulator  design  are  discussed  in  Section  7. 

The  correlation  function  of  signal  intensity  (envelope  squared)  is 
dependent  on  time  (because  of  scatterer  movements),  on  frequency  of  the  observing 
radar,  and  on  position  of  the  radar  resolution  cell  (which  governs  the  partic¬ 
ular  set  of  scatterers  which  are  observed).  Assume  that  two  clutter  measure¬ 
ments  are  performed  at  times  t  and  ^  at  frequencies  of  “>f  and  toz  respectively. 
For  the  measurement  performed  at  t,  a  set  of  scatterers  ^  ,  governed  by  the 
size  and  position  of  the  range  and  angle  resolution  cell  of  the  radar  at  , 
is  observed.  At  the  set  of  scatterers  observed  is  which  may  include 
all,  some  or  none  of  the  set  of  scatterers  depending  on  the  position  of 
the  radar  resolution  cell  at  (and,  possibly,  scatterer  movement  between  tf 
and  ) .  The  intensities  measured  at  times  z)  and  are  defined  as -T,  ,  and 
.  For  any  specific  case  of  scatterer  configuration  it  is  possible  to 
determine  -/)  and -/^  and  evaluate  the  variations  between  them  as  caused  by 
time,  frequency,  and  resolution  cell  changes.  This  result,  however,  is  so 
specific  that  it  would  be  of  relatively  limited  application  in  radar  system 

analysis  or  design.  Another  approach  is  to  determine  the  average  manner  in 

(4) 

which-/)  and  -/)  are  related.  One  of  the  averages  frequently  used  J  to  study 
the  dependence  among  variables  is  the  correlation  function,  &  ,  which  is  defined 
as : 


G(t,  ,  tz  ,  cj,  ,  ,  Nf  ,  /V2;=  <Tt  (5_1) 
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where <  > denotes  an  average  over  all  statistical  parameters  and  *1,t  is  the 
variance  defined  as 

*1,2  ~  /  *  3.  '  ^  h,  P  ^  >  (5-2) 

If  If  and  I2  are  uncorrelated,  5=0,  whereas  if  they  are  highly  correlated, 

Qzt.  Experimentally,  the  averages  are  approximated  by  forming  the  mean  of 
results  of  a  large  number  of  independent  measurements^.  It  will  be  assumed 
that  the  clutter  signal  is  stationary,  which  of  course  must  be  verified  for 
any  set  of  clutter  data.  Variations  of  clutter  signal  statistical  charac¬ 

teristics  with  time  can  be  caused  by  such  factors  as  rainfall  or  snowfall, 
atmospheric  effects  (ducting)  and  season.  These  effects  must  be  identified 
in  the  evaluation  of  clutter  data  by  meteorological  observations  (see 

(12) 

Section  8)  and  by  estimates  of  stationarity  from  the  clutter  data  itself.  J 
Assuming  stationarity,  the  correlation  function  can  be  shown  to  be  dependent 
only  on  the  time  delay,  T  =  £  -  tf  ,  between  the  measurements  rather  than  on 
the  actual  times,  tf  and  ^  .  The  effects  of  changes  in  the  environment  (hourly  - 
daily  -  seasonal)  on  radar  performance  are  small  over  the  radar  signal  process¬ 
ing  times  (of  the  order  of  seconds).  If  clutter  characteristics  change  with 
time  due  to  the  effects  mentioned  above,  the  efficacy  of  various  clutter 
rejection  techniques  will  also  vary  and  this  is  most  readily  evaluated  by 
treating  clutter  as  quasi -station ary;  that  is,  evaluate  system  performance 
assuming  stationarity  for  the  variety  of  clutter  characteristics  which  can 
occur.  It  is  extremely  important  in  a  clutter  measurements  program,  as  discussed 
in  Section  10,  to  correlate  meteorological  conditions/data  with  electromagnetic 
observable  data  to  determine  the  range  of  clutter  conditions,  their  frequency 
of  occurrence  and  causal  effects. 

In  Appendix  B,  the  correlation  function  of  Equation  5-1  is  derived 
and  shown  to  be  factorable  into  the  form: 

G(7,  U,y  eJtt  N  t  G(N)  G(r)  (5-3) 
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where  6(h)  is  the  correlation  function  which  expresses  the  effect  of  radar 
scanning  (if  and  are  identical,  G(n)  -  1  ).  Further  discussion  of  5(n) 
is  given  later  in  this  section.  6(7 )  is  the  correlation  function  which  gives 
the  effect  of  receiving  pulses  a  time  T  apart  and  6(&  oj)  describes  the  effect 
of  changing  frequency  by  4  eJ  ,  where  4  .  Note  that  stationarity  has 

been  assumed  in  the  derivation  of  Equation  5-3.  Further  assumptions  made  in 
the  derivation  are: 

1.  effects  of  multiple  scattering  and  local  masking  are  small 
(only  single  bounce  scattering  is  included) 

2.  amplitude  and  phase  of  the  individual  scatterers  is  independent 
of  frequency 

3.  relative  scatterer  locations  are  statistically  independent 
and  the  scatterers  are  randomly  distributed  over  a  region 
large  with  respect  to  wavelength 

4.  net  scatter  movement  over  the  period  will  not  remove  the 
scatterer  from  the  radar  resolution  cell 

5.  scatterer  velocities  are  statistically  independent. 

Specific  expressions  for  6  (7)  and  &(Ato)  are  given  for  uniformly  distributed 
scatterer  locations  in  Appendix  B,  Equation  B-21  or  B-22  for  uniform  or  Gaussian 
velocity  distributions,  respectively.  Multiple  scattering  and  local  masking 
makes  the  amplitude  of  the  scatterers  dependent  on  their  relative  locations. 

This  can  occur,  for  example,  if  one  scatterer  masks  the  radar  illumination  of 
another  or  if  the  scattered  signal  (bistatic)  from  one  scatterer  to  a  second 
has  approximately  the  same  amplitude  as  the  direct  radar  illumination  of  the 
second  scatterer.  If  the  scatterers  are  relatively  far  apart  with  respect 
to  their  size,  these  effects  may  be  small.  For  cases  such  as  trees  in  a 
dense  forest,  local  masking  and  multiple  scattering  are  significant  and  only 
the  tree  tops  may  contribute  as  significant  scatterers.  Multipath  scattering, 
for  example,  due  to  the  terrain  between  the  radar  and  resolution  cell,  was 
also  not  included  in  the  analysis  (see  Section  7.6). 
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Variations  of  the  amplitude  and  phase  of  the  scattering  obstacles 
in  a  radar  resolution  cell  with  changes  in  frequency  are  caused  by  mutual  phase 
cancellation  effects  among  the  scattering  centers  (assuming  a  composite 
scatterer) ,  changes  of  material  characteristics  (6/T)  with  frequency  and 
size  dependent  effects  (Rayleigh/optical  scattering  regions  and 
resonance).  Certainly  all  of  the  above  factors  are  important  in  evaluating 
relative  radar  performance  at  two  different  frequency  bands.  For  frequency 
variations  of,  say  10  percent,  only  phase  cancellation  effects  within  a  com¬ 
posite  scatterer  and  resonance  effects  would  be  of  importance  in  causing 
frequency  dependence  in  the  individual  scatterers. 

The  assumption  that  the  scatterers  are  distributed  over  a  region 
large  with  respect  to  wavelength  is  consistent  with  the  resolution  cell  size 
of  most  radars.  The  correlation  function  of  Equation  5-3  is  the  result  obtained 
when  scatterer  positions  can  be  considered  independent  for  successive  measure¬ 
ments  of  -T  and-Z^  .  If  only  a  single  set  of  scatterers  is  considered  and  the 
relative  scatterer  locations  are  constrained  to  small  movements  about  some 
mean  location,  the  complete  factorization  as  shown  in  Equation  5-3  will  not 
apply.  The  correlation  function  for  this  case  can  be  derived  by  extension  of 
the  analysis  given  in  Appendix  B. 

The  fourth  assumption,  that  scatterers  will  not  move  outside  of  the 
resolution  cell  within  the  time  7,  is  usually  satisfied  because  major  terrain 
scatterers  are  constrained  to  relatively  small  displacements  about  some  mean 
position. 


TTie  assumption  of  statistically  independent  scatterer  velocities  is 
reasonable  considering  the  variations  in  wind  conditions,  which  are  the  major 
source  of  scatterer  movement,  over  a  typical  resolution  cell. 

It  should  be  noted  that  the  closed-loop  clutter  simulator  discussed 
in  Section  7  is  not  constrained  by  the  assumptions  described  above  and  also 
made  in  the  derivation  of  Equation  5-3.  Thus,  the  sensitivity  of  the  clutter 
characteristics  and  correlation  functions  could  be  evaluated.  Because  open- 
loop  simulation,  discussed  in  Appendix  E,  is  not  as  closely  related  to  basic 
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scattering  characteristics  as  closed-loop  simulation,  it  is  difficult  to  as¬ 
certain  whether  open-loop  simulation  is  constrained  by  the  assumptions  1  through 
5  without  further  study  of  their  effects  on  clutter  signal  characteristics. 

5.3  DISCUSSION  OF  RESULTS 

Numerous  references  are  available  in  which  correlation  functions  have 
been  derived  for  the  frequency-shift,  moving-scatterer,  and  moving-gate 
situations.  In  most,  if  not  all,  of  these  references,  the  many-scatterer 
assumption  is  made.  In  this  section,  some  of  these  earlier  results  are 
considered  and  compared  with  those  obtained  here. 

First,  consider  the  effects  of  a  frequency  shift;  assume  that  the 

observed  region  remains  the  same  and  that  scatterer  velocity  is  zero. 
fl31 

Goldstein v  '  derives  an  expression  for  the  mean-square  difference  in  signal 
power  when  there  is  a  change  in  frequency  .  His  result  contains  a  misprint: 
in  his  equation  134,  the  exponential  2  should  be  inside,  not  outside,  the 
square  bracket.  Although  his  (pz'pt)  is  not  the  same  as  the  correlation 
function  we  found  in  the  preceding  section,  the  similarity  in  significance  and 
functional  form  can  be  seen.  McGinn  and  Pike  (Reference  12,  p.  60)  claim  that 
Goldstein's  result  is  seriously  in  error.  This  contention  does  not  appear  to 
be  justified.  'Hie  term  arising  from  equality  of  the  four  subscripts  has  been 
left  out,  as  they  say,  but  this  term  is  zero  anyway.  Failure  to  satisfy  McGinn 
and  Pike's  electrodynamic  requirement  (that  there  be  no  net  dc  level  in  the 
pulse)  should  be  relatively  unimportant  for  any  pulse  more  than  a  few  wave¬ 
lengths  long  or  for  any  instantaneous  signal  bandwidth  less  than  about  20 
percent.  Contributions  arising  from  the  net  dc  level  are  analogous  to  the 
terms  dropped  from  Equation  B-18  and  used  in  obtaining  Equation  B-19 
(Appendix  B) .  McGinn  and  Pike  say  that  some  Monte  Carlo  checks  indicate  an 
error  "as  great  as  10^"  in  Goldstein's  result;  whether  they  were  using  the 
corrected  form  of  his  equation  134  is  not  stated,  nor  is  the  quantity  equal  to 

10°.  On  the  basis  of  work  presented  here,  the  objection  raised  by  McGinn  and 

(14) 

Pike  thus  appears  invalid.  Wallace  also  derives  a  correlation  function 
for  frequency  shifting  assuming  insignificant  target  motion  during  the  times 
involved.  His  result,  obtained  for  the  case  of  many  scatterers,  is  the  same 
as  that  presented  here. 
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Next,  consider  the  moving-scatterer  case.  Fleisher^^ 
many  moving  scatterers  and  derives  a  correlation  function  similar  to  that 
derived  here.  The  methods  of  derivation  differ  significantly  but  his  equation  19 
reduces  to  our  form.  Fleisher  also  shows  that  the  correlation  function  for 
signal  intensity  (which  we  found  in  the  preceding  sec  ion)  is  the  square  of 
the  correlation  function  for  envelope  voltage  level  (?)  if,  and  only  if,  -V  is 
large  and  no  single  scatterer  dominates.  These  requirements  are  simply  those 
for  the  Central  Limit  Theorem  to  apply.  Therefore,  the  correlation  function 
of  the  envelope  (?)  is  the  square  root  of  the  correlation  of  the  intensity  (O 
if  the  envelope  is  Rayleigh,  but  as  *vt  nothing  can  be  said  of  the  correlation 
function  of  H  for  other  distributions. 


Finally,  consider  the  effect  of  a  moving  range  bin;  the  corresponding 
correlation  function  S(y)  war  found  in  Equation  B-14  (Appendix  B) .  There  is  no 
requirement  that  there  be  many  scatterers  present.  Consider,  as  a  special 
case,  a  situation  in  which  there  are  many,  homogeneously  distributed,  equal - 
amplitude  scatterers  in  the  region  observed  by  the  moving  range  bin.  Equation 
B-14,  for  this  case,  reduces  to: 

t 

6  (AS)  *  — ^  (5-4) 

where  M  is  the  total  number  of  scatterers  and  4^  is  the  number  of  scatterers 
common  to  both  measurements.  This  result  agrees  with  Rogers when  the 
appropriate  change  of  variables  is  made:  he  is  concerned  with  envelope  rather 
than  intensity,  so  his  correlation  function  is  the  square  root  of  that  derived 
here  (recall  Fleisher's  result  stated  above),  and  he  deals  in  relative  shift 
of  range-gate  po  ition  rather  than  in  numbers  of  scatterers.  His  result 
assumes  range  bins  that  are  not  extremely  short,  as  does  ours. 


It  is  also  of  interest  to  consider  how  fast  6CMJ  approaches  the 
value  it  will  have  for  many  scatterers  as increases.  Suppose  that  all 
scatters  have  unit  amplitude  and  that  the  scatterers  are  homogenously  dis¬ 
tributed;  then,  we  have:  . //y  i  \ 

__  (Tjit-Jf) 


G(n) 


N(N-l) 


1  ~  N 


(5-5) 
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If  we  assume  the  scatterers  to  be  homogeneously  distributed,  represents  the 

fraction  of  the  range  bin  that  is  common  to  the  first  and  second  measurements; 

it  is,  therefore,  reasonable  to  fix  ,  which  is  a  geometrical  parameter,  and 

to  let  M  (and,  thus,  V,  )  increase.  In  Figure  5-1  are  shown  some  examples  from 

b  f  Ml  \  * 

which  it  can  be  seen  that  the  limiting-case  value  of  \-jf)  is  approached 
relatively  slowly  when  ^  is  small.  It  is  interesting  to  note  that  when 

/if 

M,  >  4  or  5  the  value  of  G(M)  is  nearly  equal  to  the  -  °°  value,  at  least 

Ml 

for  the  range  of  values  of  -j*  considered  here.  This  result  is  plausible  from 
inspection  of  Equation  5-5. 

For  situations  involving  unequal  scatterers,  especially  when  one,  or 
a  few,  of  the  scatterers  dominate,  G(n )  does  not  reduce  to  an  easily  interpreted 
form.  The  situation  becomes  especially  hard  to  interpret  when  there  are 
relatively  few  scatterers.  On  the  basis  of  the  analysis  which  led  to  it,  it 
seems  reasonable  to  expect  Equation  B-14  (Appendix  B)  to  be  accurate  for  any 
set  of  scatterers;  the  assumptions  that  were  made  relative  to  negligible  terms 
in  our  derivation  of  in  no  way  affect  G(n).  The  only  effect  neglected 

in  our  derivation  arises  from  scatterer  motion  across  boundaries,  which  should 
be  small  for  terrain  clutter. 

In  the  preceding  section  and  the  discussion  above,  it  has  been  shown 
that  the  correlation  functions  for  target  motion  (time  correlation) ,  frequency 
shift,  and  range-bin  shift,  under  the  set  of  assumptions  previously  described, 
are  independent  quantities  that  are  simply  multiplied  together  to  obtain  the 
overall  correlation  coefficient.  It  must  be  remembered,  however,  that  when 
only  a  few  scatterers  are  present,  measurements  are  difficult  to  interpret, 
because  measured  quantities  will  in  such  cases  tend  to  depart  further  from 
their  mean  values  than  when  there  are  many  scatterers  (for  an  example  of  such 
difficulties  arising  in  the  field  of  radar  meteorology,  see  References  17  and 
18). 
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Figure  5-1  CORRELATION  FUNCTION  G(N)  FOR  RADAR  RESOLUTION  CELL 
MOVEMENT  VERSUS  NUMBER  OF  SCATTERERS  (N) 
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Section  6 
DOPPLER  EFFECTS 


In  the  two  previous  sections  (4  and  5)  the  probability  density 
functions  and  the  correlation  functions  were  discussed  for  clutter  signals 
arising  from  a  (not  necessarily  large)  number  of  scatterers.  The  probability 
density  function  for  the  amplitude  of  the  signal  (envelope)  was  shown  in 
Section  4  to  deviate  strongly  from  the  Rayleigh  density  function  when  there 
are  few  scatterers,  particularly  when  one  scatterer  dominates  the  rest.  In 
Section  5,  correlation  functions  relating  squared  signal  amplitudes  (i.e.,  the 
outputs  of  square- law  detectors)  were  found  to  describe  effects  of  frequency 
shift,  scatterer  motion,  and  range-gate  motion.  As  pointed  out  in  Section  5, 
the  probability  densities  for  linearly  detected  signal  amplitudes  can  be 
converted  into  probability  densities  for  square-law-detected  signal  amplitudes. 
The  inverse  operation  for  correlation  functions  (i.e.,  finding  the  correlation 
function  for  a  linearly  detected  signal  from  that  for  a  square-law-detected 
signal)  is  not  practical,  in  general,  except  in  the  many-scatterer  case  for 
which  the  probability  density  function  of  the  linearly  detected  signal  is 
Rayleigh. 


Results  in  Section  4  were  obtained  for  both  coherent  and  noncoherent 
radars,  whereas  Section  5  was  concerned  primarily  with  noncoherent  radars, 
although  many  of  the  results  of  Section  5  apply  to  a  coherent  radar  as  well. 

If  coherent  (Doppler)  radars  are  of  concern,  then  the  simulated  clutter  signal 
must  have  basically  the  same  amplitude  and  phase  spectral  characteristic  as 
the  actual  clutter  signal  would  have.  Effective  simulation  for  the  Doppler- 
radar  case  will  be  valid  for  a  noncoherent  radar,  but  the  converse  is  not 
necessarily  true. 

th 

Consider  a  set  of  A/  scatterers.  The  k  scatterer  is  at  range  at 
time  t^O  and  has  a  velocity  component  along  the  radar  line  of  sight  (a^ 
positive  for  a  velocity  that  increases  ) .  The  radar  operates  at  a  (radian) 
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frequency  cJe  (frequency  So  -  j^.)  with  free-space  wavelength  xo=j~,Cm  velocity 
of  propagation.  The  X  **  scatterer  has  a  scattering  amplitude  that  is  pro¬ 
portional  to  the  square  root  of  the  radar  cross  section  of  that  scatterer.  It 
is  assumed  that  electromagnetic  interaction  between  scatterers  in  the  radar 
resolution  cell  is  negligible;  for  most  cases  of  interest,  this  assumption  is 
reasonable.  Multipath  effects,  not  considered  here,  are  however  generally 
important  and  are  discussed  in  Section  7.6.  It  is  also  assumed  that  remains 
constant  during  a  measurement.  The  received  signal  is 

n  [  r  f  u  t  1 

E(t)  =  Z  \  co s  2  ± - !Z-  cJ0-cJ0t  (6_j) 

*./  L  c  J 


The  following  treatment  is  for  integration  over  all  time  and  embodies  the 
assumption  that  the  velocities^  and  amplitudes  Ak  are  constant  for  all  time. 
Practical  restrictions  will  impose  restrictions  on  observation  time  and  involve 
velocities  varying  in  some  degree.  These  restrictions  cause  individual  spec¬ 
tral  lines  to  broaden.  The  Fourier  transform  of  E  (t) is 


"  n  f  a(*>r*  ■3-u»rk 

v.  1 _ 4  '  Af\  ' 


/f*/  (  *• 

This  power  spectrum,  defined  from  o*oto  v*oo  ,  is  then 


(6-2) 


k--t  ° 


(6-3) 

from  which,  converting  to  conventional  frequency,  we  obtain 

° 

(6-4) 


In  deriving  Equation  6-4  it  is  assumed  that  all  u ^  are  distinct;  if  several 
scatterers  have  identical  u ^  values  (i.e.,  =  u*  =  ...)  then  must  be 
interpreted  as  the  phasor  sum  of  the  individual  scattering  amplitudes 
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Note  that  nowhere  in  the  derivation  of  Equation  6-4  was  N  assumed  to 
be  large;  it  is  only  necessary  that  N  remain  fixed  during  the  time  of  the 
measurement.  It  is  seen  from  Equation  6-4  that  ^Yfjhas  values  at  discrete 

t  U.  ft 

frequencies  that  depend  upon  scatterer  velocity  u ^  .  The  quantity  is 
simply  the  Doppler  frequency  shift  associated  with  the  k  scatterer.  If  all 
scatterers  had  the  same  amplitude  =  A  for  all/c)  then  P(f)  would  have  the 
same  shape  as  the  velocity  spectrum.  If  the^  are  unequal ,  p(f)  is  a  weighted 
velocity  spectrum. 

Similar  results  have  been  obtained  by  Rogers  (Reference  18,  p.  12) 
for  the  large- N  case.  His  result  is  a  continuum- spectrum  form  analogous  to 
our  Equation  6-4. 

The  basic  concept  on  which  a  doppler  radar  separates  the  returns  from 
scatterers  having  different  velocities  along  the  line  of  sight  is  presented  in 
the  above  discussion.  Two  limitations  must  be  considered.  First,  note  that  in 
practice  the  power  spectrum  P(-P)  is  not  centered  on  f0  but  on  some  much  lower 
frequency  obtained  from  one  or  more  heterodyning  operations .  A  spectrum  folding 
limitation  frequently  arises  due  to  the  multiple-line  structure  in  the  transmitted 
waveform  —  the  typical  PRF  harmonics  for  instance.  Also,  if  the  lower  frequency 
is  less  than  max  \—rK  folding  may  occur  and  negative  velocities  will  be  confused 
with  positive  velocities.  In  the  limit,  HP)  might  be  centered  on  0  frequency  (a 
situation  that  would  arise  if  the  signal  were  beat  against  a  local  oscillator 
operating  at  )  and  then,  for  a  system  employing  a  Doppler  frequency  measure¬ 
ment,  all  negative  velocities  would  be  reinterpreted  as  positive  velocities. 

In  such  systems  the  positive  or  the  negative  part  of  the  velocity  spectrum  can 
be  found  only  if  the  velocity  spectrum  is  symmetrical  (Reference  18,  p.  16). 

A  second  limitation  on  measurement  of  scatterer  velocities  is  the  resolution 
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limit  imposed  by  a  finite  pulse  length;  uncertainty  in  velocity  measurements 

C 

is  of  the  order  of  (zfcr  )  *  where  Te  is  the  effective  length  of  the  pulse 
(TgSsT  ,  the  pulse  duration,  if  the  pulse  is  approximately  rectangular  and 
multiple  pulse  integration  is  not  employed)  (Reference  19,  p.  203).  However, 
if  a  long  train  of  coherent  pulses  are  used,  one  can  obtain  an  arbitrarily 
high  degree  of  velocity  resolution  against  "ideal"  scatterers  (for  high  ~  , 
velocity  resolution  would  be  independent  of  pulse  length  but  dependent  on 
total  observation  time) . 

To  what  extent,  then,  must  Doppler  characteristics  be  simulated? 
Certainly,  if  a  coherent  (pulse-Doppler)  radar  is  to  be  evaluated  using  a 
simulated  clutter  signal,  it  is  necessary  that  the  spectral  characteristics  of 
the  clutter  signal  match  reasonably  well  those  of  the  clutter  being  simulated. 

If  the  clutter  signal  were  simulated  by  a  wide-band  noise  source,  the  spectrum 
of  the  clutter  signal  would  be  determined  by  the  bandpass  characteristics  of 
the  receiver.  The  apparent  velocity  spectrum  of  the  clutter  would  then  be  as 
established  by  the  receiver  characteristics  independently  of  the  actual  clutter 
characteristics.  This  inadequacy  is  typical  of  a  simulator  that  attempts  to 
use  a  noise  source  to  represent  a  clutter  signal;  further  discussion  of  this 
situation  is  given  in  Section  7.  Accurate  simulation  of  the  clutter-signal 
spectrum  is  of  greatest  importance  if  the  radar  to  be  evaluated  uses  signal  tem¬ 
poral  characteristics  to  discriminate  between  target  and  clutter.  For  example, 
suppose  a  Doppler  radar  employed  characteristics  of  the  phase  or  amplitude  of 
the  cross  section  versus  time  of  aircraft  or  missiles  to  enhance  detection  in 
the  presence  of  clutter.  It  would  obviously  be  necessary  in  this  case  to  use 
much  more  representative  clutter  signals  for  simulation  purposes  than  could 
be  generated  from  filtered  white  noise. 

In  the  foregoing  discussion  it  has  been  (inherently)  assumed  that  the 
radar  is  stationary  and  that  each  scatterer  has  a  (momentarily)  constant  velo¬ 
city  component  along  the  radar  line  of  sight.  If  the  radar  is  in  motion, 
additional  Doppler  shifts  result.  Although  these  effects  can  be  described  by 
Equation  6-4  by  including  radar  motion  in  the  definition  of  u ^  (which  thus 
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becomes  the  relative  velocity  between  the  radar  and  the  & scatterer) ,  it 
may  at  times  be  convenient  to  consider  the  Doppler  shift  as  arising  from  two 
motions:  first,  the  gross  motion  of  the  radar  toward  a  fixed  reference  near 
the  clutter-generating  region,  and,  second,  the  motions  of  individual  scatter- 
ers  relative  to  the  first  reference.  For  ground  clutter,  the  radar-motion 
effect  dominates  in  the  case  of  fixed  scatterers  (rocks,  ground  contour, 
buildings,  tree  trunks,  etc.);  the  second  effect  arises  when  there  is  motion 
of  scatterers  relative  to  the  ground  (moving  vehicles,  wind-moved  foliage,  etc.). 
The  radar-motion  component  is  geometrical  and  can  be  computed,  approximately, 
using  simple  trigonometry.  Because  (typically)  ground  clutter  is  observed 
over  a  range  of  elevation  and  azimuth  angles  the  clutter  signal  is  spread  over 
a  spectral  width  dependent  on  the  angular  width  and  orientation  (with  respect 
to  the  velocity  vector)  of  the  radar  beam.  Scatterer  motion  produces  additional 
Doppler  shifts  that  may  be  important,  depending  on  the  system.  Additional 
spectral  components  can  arise  from  antenna  side lobes  that  illuminate  terrain 
at  angles  (and,  consequently,  at  relative  velocities)  much  different  (usually) 
from  those  prevailing  along  the  main  lobe. 

If  a  clutter  simulation  is  to  be  used  to  evaluate  a  coherent  (e.g., 
pulse-Doppler)  radar,  it  is  necessary  that  the  spectral  characteristics  of  the 
simulated  clutter  signal  match  those  of  the  actual  clutter  signal.  It  is 
shown  here  (Equation  6-4)  that  the  power  spectrum  of  the  clutter  signal  is  a 
weighted  form  of  the  velocity  spectrum  regardless  of  the  number  of  scatterers 
contributing  to  the  signal  provided  multiple  scattering  effects  are  small. 
Separation  of  velocity  into  two  parts,  one  for  radar  motion  relative  to  a 
fixed  point  on  the  ground  and  one  for  motion  of  scatterers  relative  to  the 
ground,  may  prove  useful  for  some  forms  of  simulation,  since  radar  motion  is 
thus  specified  independently  of  scatterer  motion. 
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Section  7 

CLUTTER  SIMULATION 


7.1  GENERAL 

In  the  following  section  the  simulation  of  clutter  and  clutter/target 
interactions  for  evaluation  of  radar  systems  are  discussed.  Clutter  simulation 
techniques  are  based  on  some  form  of  clutter  model  -  that  is,  on  a  method  of 
representing  the  characteristics  of  the  clutter  signal  and  the  changes  in  the 
clutter  due  to  changes  in  illuminator  and/or  observer  radar  parameters.  It  is 
assumed  that  only  a  relatively  refined  clutter  simulation  is  of  interest,  and, 
thus,  little  attention  is  given  to  the  use  of  additive  filtered  noise  to 
represent  a  complete  clutter  signal,  although  some  of  the  limitations  of  this 
approach  are  discussed. 

The  generic  forms  of  open-loop  and  closed-loop  clutter  models  as 
described  in  Section  3.2,  also  see  Figures  3-1  and  3-2,  are  used  to  categorize 
simulation  techniques  for  clutter  signals.  The  use  of  additive  noise  to 
simulate  clutter  is  an  open-loop  simulation;  that  is,  the  properties  of  the 
clutter  signal,  probability  density  function,  average  power  and  spectrum  in 
this  case,  must  be  known  (or  assumed  known)  for  the  particular  radar  of  interest. 
The  clutter  signal  having  the  desired  characteristics  is  formed  and  added  to 
the  target  signal  (see  Figure  3-1)  for  evaluation  of  the  radar. 

7.2  BASIC  CONSIDERATIONS 

For  either  open-loop  or  closed-loop  clutter  models,  two  implementations 
as  a  simulator  are  possible:  either  analog  or  digital  computation  might  be 
used  (or,  possibly,  hybrid  computation,  although  this  technique  will  not  be 
discussed  explicitly  here).  There  are,  therefore,  four  possible  simulation 
techniques  discussed  in  this  section:  open-loop  simulation,  analog  and  digital, 
and  closed-loop  simulation,  analog  and  digital.  First,  however,  some  basic 
properties  of  these  techniques  are  considered. 
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7.2.1  Open  Loop  versus  Closed  Loop 


Closed- loop  simulation  has  significant  advantages  over  open- loop 
simulation  and  its  use  for  all  but  the  simplest  of  simulators  seems  assured. 
The  principal  advantage  is  that  a  well-designed,  closed-loop  simulator  would 
function  for  a  wide  variety  of  radar  signals  without  internal  adjustment, 
whereas  in  an  open-loop  simulator  it  would  be  necessary  to  adjust  both  the 
clutter-signal  and  the  target  signal  generators  for  any  change  in  radar  param¬ 
eters.  Both  techniques  require  about  the  same  amount  of  computation,  in  spite 
of  the  apparent  relative  simplicity  of  the  open-loop  simulator,  because  of  the 
need  for  additional  computation  for  proper  adjustment  of  the  signal  generators 
This  situation  becomes  even  more  serious  when  it  is  necessary  to  include  such 
effects  as  masking  and  multipath. 

7.2.2  Analog  versus  Digital 

The  use  of  analog-computation  simulation  is  quickly  found  to  have 
significant  disadvantages  when  it  is  compared  with  digital-computation  simul¬ 
ation.  The  main  disadvantage  is  a  requirement  for  a  large  number  of  signal 
generators  with  appropriate  controls  for  amplitudes  and  phases  (and,  possibly, 
frequency  as  well) .  Frequency  scaling  to  conventional  analog-computer 
frequencies  is  not  practical,  except  in  special  circumstances  (such  as  investi 
gation  of  pulse-shape  distortion  for  a  very  limited  number  of  pulses),  because 
of  the  attendant  increase  in  simulation  time  (thus  if  frequency  is  scaled  by 
10  ,  computation  time  is  increased  by  10°) .  It  therefore  may  be  necessary  to 

operate  at  least  part  of  the  simulator  at  high  frequencies;  although  feasible, 
one  of  the  conveniences  of  normal  analog  computation  is  lost. 

Because  the  digital-computation,  indirect-generation  simulator 
appears  to  be  the  most  flexible  type,  this  form  is  discussed  in  the  next 
section  in  some  detail.  Considerations  involved  in  inclusion  of  polarization 
effects,  masking  effects,  multipath  effects,  and  bistatic-system  effects  into 
the  closed-loop  simulator  are  included.  A  discussion  of  closed-loop  analog 
simulation  and  open-loop  simulation  (analog  and  digital)  is  included  in 
Appendix  E. 
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7.3 


CLOSED- LOOP  SIMULATION  -  DIGITAL  COMPUTATION 


As  noted  in  Section  7.2,  the  combination  of  digital  computation  with 
closed- loop  simulation  of  clutter  (and  target)  signals  offers  a  maximum  of 
flexibility.  Exact  computation  of  waveform  of  radar  signals  scattered  from 
clutter  is  obviously  impossible.  For  radars  and  pulses  of  commonly  used  band- 
widths,  useful  results  can  be  obtained  using  straightforward  approximations, 
provided  the  scatterers  producing  the  clutter  signal  satisfy  certain  restrictions 

1.  Scattering  is  independent  (interaction  among  scatterers 
can  be  ignored);  this  approximation  is  reasonable,  but  in 
some  cases  it  may  not  be  satisfied.  Inclusion  of  multiply 
scattered  signals  is  possible  in  the  simulation  and  is 
discussed  in  Section  7.3.4. 

2.  Scattering  from  a  particular  scatterer  is  not  a  rapidly 
varying  function  with  frequency.  Thus,  scattering  from 
highly  resonant  (within  radar  band)  objects  (e.g.,  dipoles) 
may  need  modified  treatment,  although  again  there  is  no 
difficulty  in  including  such  effects  in  the  simulation  if 
the  appropriate  parameters  are  known  (see  Section  7.3.5). 

3.  Acceleration  of  scatterers  is  small  enough  so  that,  during 
the  time  involved  in  a  range  scan,  the  velocity  of  a 
scatterer  can  be  assumed  constant. 

Removal  of  these  restrictions  is  discussed  after  the  more  simple  case,  in  which 
they  are  satisfied,  has  been  treated. 

7.3.1  Computation  of  Clutter  Signal  Envelope 

First,  assume  that  restrictions  1  through  3  are  satisfied  and  that 
the  clutter  signal  arises  from  discrete  scatterers  in  the  antenna  beam.  For 
the  moment,  assume  that  the  antenna  does  not  scan  and  that  there  are  N  signifi¬ 
cant  scatterers  along  the  ground  in  the  range  of  interest.  If  the  transmitted 
radar  signal  is: 
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*  c,  (t)  cos  ( cj0  i  ) 


(7-1) 


where  ef(i)is  an  envelope  function  that  varies  slowly  relative  to  the  (radian) 
carrier  frequency  oJ0  ,  then  it  can  be  shown  (Appendix  C)  that  the  signal 
received  from  the  N  scatterers,  if  they  are  stationary,  is: 

e,f*-rk)cas\<Jc(t-Tk)]  (7-2) 

where  is  the  amplitude  of  the  return  from  the  /c scatterer  and  7 ^  is  the 
(2-way)  time  delay  corresponding  to  propagation  from  the  radar  to  the 
scatterer  and  back.  Therefore,  the  appropriately  scaled  and  time-shifted 
envelopes  can  be  added  together,  provided  the  relative  RF  phases  at  the  carrier 
frequency  are  taken  into  account.  Equation  7-2  holds  regardless  of  the  band¬ 
width  of  the  signal  (i.e.,  it  is  not  actually  necessary  for  ef(t)  to  vary 
slowly  relative  to  u)Q  for  Equation  7-2  to  hold). 

7.3.2  Receiver  Filtering  Effects 

If  the  receiver  has  a  bandpass  characteristic  that  distorts  the 
received  waveform,  as  is  commonly  the  case,  then  a  modification  to  the  above 
result  (and  that  of  Appendix  C)  is  required.  Provided  ef{t)  varies  slowly 
enough  relative  to  oJ0  so  that  spectral  folding  is  not  significant,  and  provided 
the  receiver  bandpass  is  symmetrical  about  u)0 ,  the  received  waveform  can  be 
found  from  Equation  7-2  with  being  replaced  by  -  ?^),  where  ejt) 

is  the  waveform  obtained  when  (f),  the  envelope  function,  is  passed  through 
the  low-pass  equivalent  of  the  bandpass  filter  in  the  receiver  (i.e.,  the 
filter  characteristic  obtained  by  shifting  the  bandpass  from  its  center 
frequency  in  the  receiver  to  0).  If  the  bandpass  filter  (receiver)  is  not 
symmetrical  about  o)0  ,  a  more  complicated  analysis  becomes  necessary.  The 
nonsymmetrical  part  of  the  filter  response  now  produces  a  quadrature  phase 
signal  component.  Two  equivalent  low-pass  filters  are  also  now  required,  one 


Spectral  folding  is  unimportant  for  most  radars;  only  in  an  extremely  short- 
pulse,  high-resolution  radar  will  it  be  significant. 
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to  find  amplitude  of  the  in-phase  component,  and  the  other,  the  amplitude  of 
the  quadrature  component.  Although  some  additional  computation  is 

required,  there  is  no  fundamental  difficulty  introduced  by  an  offset  bandpass 
filter,  and  the  basic  method  discussed  can  still  be  applied. 


7.3.3 


>ler  Effects 


Another  modification  to  Equation  7-2  is  required  when  the  scatterers 
are  moving  and  thereby  producing  different  Doppler  shifts.  If  the 
scatterer  has  velocity  component  uk  along  the  line  of  sight  (positive  in  the 
direction  away  from  the  radar) ,  and  if  Tk  now  represents  the  delay  associated 
with  the  k  scatterer  at  t  =  0,  Equation  7-2  becomes: 

tT(t)  =Eak  et(i'Tk  “  T*  *)**’  \£Jo(f~  t'etJork\  (7-3) 

k*f 

No  longer  is  the  straightforward  addition  of  properly  phased  contributions 
possible,  because  the  phase  shift  associated  with  the  carrier  now  is  itself 
time  varying  (at  the  Doppler- frequency  rate).  Each  of  the  pulses  added  together 
is  thus  at  a  different  frequency  (assuming  distinct  u ^  values).  The  envelope 
is  therefore  time  varying  in  a  much  more  complicated  way  than  when  all  u^  =  0. 

We  can  write  Equation  7-3  as: 


a-k\e'&-rk-irit\t2T.  ’-%*■*)• 
k=f  [  J  A*2 


•COS  2JT^  (Uk~um  )+CJo  fa  ~  7-m)  •  C0S 

J  J  L  -J 


(7-4) 


where  <f)  (t)  is  a  phase  modulation,  the  envelope  being  the  quantity  in  curly 
brackets  (the  exponent  1/2  is  dropped  if  a  square-law  detector  is  used) .  For 
scatterers  more  than  a  pulse  width  apart,  the  interaction  term  drops  out.  If 
there  is  a  large  constant  velocity  superimposed  on  the  ^values  (as  when 
ground- clutter  measurements  are  made  from  an  airborne  radar) ,  Equation  7-4 
should  be  modified  by  letting  u>0  be  the  carrier  frequency  Doppler  shifted  by 
the  constant  velocity;  u ^  then  are  variations  in  the  velocity  relative  to  the 
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reference  value  used  to  establish  u)e  .  Note  that  even  in  this  case  computation 
need  not  be  at  the  RF  rate,  since  the  envelope  of  the  output,  which  is  being 
computed,  varies  at  a  much  lower  rate,  in  general.  Although  this  method  is 
not  exact,  because  of  the  ambiguity  introduced  by  one's  choice  of  carrier 
frequency  in  the  derivation,  Equation  7-4  should  give  useful  results  if  some 
care  is  exercised.  A  more  exact  method  would  be  to  compute  e.r(t)  from 
Equation  7-3  and  use  some  form  of  (numerical)  envelope  detector.  Possibilities 
include:  1)  the  finding  of  peaks  of  the  waveform  and  assigning  these  values  to 
the  envelope  and  2)  the  procedure  of  squaring  e-r(t)oT  finding  its  absolute 
value,  if  a  linear  detector  is  desired,  and  passing  the  result  through  a 
(numerical)  low-pass  filter.  Such  methods  are  more  realistic  but  require  far 
more  computer  time. 

The  discussion  in  the  preceding  paragraph  applies  to  an  incoherent 
radar  that  envelope  detects  the  received  signal.  A  Doppler  radar  uses  filters 
to  separate  components  of  different  frequencies,  and,  hence,  the  portions  of 
the  signal  received  from  scatterers  having  different  velocities  (Section  6). 

varies  much  more  slowly  than  (J0,  one  can,  to  a  first  order,  assign 
those  scatterers  for  which  J  lies  in  a  particular  Doppler-filter  band¬ 

width  to  that  channel  of  the  receiver.  Each  channel  will  then  have  a  signal 
formed  by  a  subset  of  the  N  scatterers.  If  the  signals  in  the  Doppler  channels 
are  then  detected  incoherently,  the  above  methods  can  be  applied  separately  in 
each  channel.  If  e,(t)  is  short  (contains  few  cycles  of  &)0)  there  is  a 
broadening  of  the  spectrum  associated  with  each  u ^  value  that  must  be  taken 
into  account;  this  broadening  is  a  cause  of  the  Doppler  ambiguity,  an  uncer- 
tainty  in  velocity  measurement  of  the  order  of  .  ,  where  f  *  —  and  Te  , 

the  effective  pulse  length,  is  of  the  order  of  the  pulse  duration  if  the  pulse 
is  approximately  rectangular  (Reference  19,  p.  203).  Inclusion  of  these  effects 
in  the  clutter  simulation  could  take  at  least  two,  and  probably  more,  forms. 
First,  one  could  compute  the  RF  waveform  and  use  numerical  filtering  to  sort 
out  which  signals  appear  in  which  channel.  This  method  is  accurate  for  any 
spectral  characteristics  but  involves  a  large  number  of  computations. 
Alternatively,  analytical  expressions  for  the  amount  of  spectral  spreading 
could  be  developed  and  used  to  modify  the  first-order  approximation  in  which 
no  spreading  was  assumed.  Of  greatest  significance  here  is  the  fact  that  the 
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form  of  simulation  may  have  to  be  changed,  depending  on  the  form  of  processing 
involved;  the  alternative  would  be  to  produce  the  complete  waveform,  which 
would  be  usable  in  any  processor,  but  >*hich  would  be  expensive  in  computer 
time.  It  is  clear  that  further  work,  as  described  in  the  program  plan,  is 
required  before  an  optimum  technique  can  be  recommended  for  simulation  of 
clutter  for  a  coherent  (pulse-Doppler)  radar. 

7.3.4  Multiple  Scattering 

The  discussion  so  far  has  been  based  on  the  assumption  that  restric¬ 
tions  1  through  3  at  the  beginning  of  this  section  are  satisfied.  Let  us 
consider  the  effect  of  relaxing  these  restrictions.  Restriction  1  can  be 
violated  in  at  least  two  ways.  If  interaction  of  scatterers  is  a  simple, 
multiple-bounce  effect  whereby  the  wave  striking  one  scatterer  is  diffracted 
to  another  and  thence  back  to  the  radar,  a  straightforward  extension  of  the 
methods  discussed  above  is  possible.  It  is  only  necessary  to  introduce  a 
fictitious  scatterer  in  addition  to  each  member  of  an  interaction  pair;  this 
fictitious  scatterer  has  a  delay  corresponding  to  the  total  propagation  path 
and  an  amplitude  that  depends  on  the  appropriate  bistatic  cross  sections  of 
the  intersecting  scatterers  and  the  range  between  them.  The  difficulty  in 
implementing  this  computation  lies  with  the  assignment  of  parameters,  not  in 
the  computation  itself.  If,  on  the  other  hand,  the  scatterers  are  close  enough 
that  they  are  electromagnetically  coupled,  no  simple  model  will  suffice.  It 
is  probable  that,  in  this  case,  the  interacting  pair  could  best  be  treated  as 
a  single  scatterer  with  (possibly)  a  rapidly  varying  frequency  response.  This 
situation  then  would  correspond  to  a  violation  of  restriction  2. 

7.3.5  Scatterer  Frequency  Dependence 

If  a  scatterer  has  a  scattering  amplitude  that  varies  rapidly  with 

frequency,  in  violation  of  restriction  2,  the  system  response  can  be  found  by 

assuming  the  frequency- response  characteristic  of  the  scatterer  to  be  combined 

* 

with  the  filtering  characteristic  of  the  receiver.  The  resulting  (generally) 

* 

It  is  necessary,  of  course,  that  a  suitable  scattering  description  of  the 

scatterer  be  available. 
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nonsymmetrical  filter  can  then  be  treated  using  in-phase  and  quadrature 
equivalent  low-pass  filters  as  discussed  above.  For  each  such  scatterer,  then, 
two  distorted  pulse  envelopes  are  fou.  1,  one  being  an  envelope  for  an  in-phase 
carrier,  and  the  other,  an  envelope  for  a  quadrature  carrier.  The  resultant 
envelope  of  the  received  signal  can  be  computed  as  before. 

7.3.6  Scatterer  Acceleration 


Failure  to  satisfy  restriction  3  in  a  noncoherent  radar  causes  little 
difficulty  if  the  RF  waveform  is  to  be  computed  and  detected.  If  Equation  7-4 
is  to  be  used,  variable  u can  be  included  straightforwardly  with  (probably) 
little,  if  any,  error.  Restriction  3  is  important  if  a  coherent  (pulse-Doppler) 
radar  is  involved,  because  the  change  in  velocities  will  cause  smearing  of  the 
Doppler  spectrum.  The  feasibility  of  including  variable  uk  in  the  computation 
depends  upon  the  type  of  computation  being  made,  which  in  turn  depends,  as 
noted  above,  upon  the  use  to  be  made  of  the  simulation.  If  the  time-varying 
RF  waveform  is  computed  from  Equation  7-3,  ^  can  be  made  a  variable  with 
little  additional  effort.  If  analytical  expressions  for  spectral  spreading 
are  used  to  modify  the  first-order  approximation,  additional  modification  to 
compensate  for  spectral  spreading  resulting  from  nonconstant  velocities  must 
be  made.  Thus,  to  some  extent  at  least,  the  restrictions  can  be  relaxed  at 
the  expense  of  added  computations. 

7.3.7  Range  Gate 

We  have  thus  seen  that,  if  some  reasonable  restrictions  are  met, 
relatively  straightforward  computation  will  yield  a  simulated  clutter  signal, 
er(t),  as  a  function  of  time,  provided  suitable  values  can  be  found  for  the  a ^ 
and  parameters.  If  the  receiver  is  range  gated,  er(t)  can  be  simply 
multiplied  by  a  time  function  corresponding  to  the  range  gate;  in  its  simplest 
form,  the  range  gate  accepts  values  of  c.r(t) for  t  within  specific  limits  and 
makes  &r(t)  zero  for  times  not  within  these  limits  (rectangular  range  gate). 

If  there  are  several  7^  values  in  each  range-resolution  cell,  the  (simulated) 
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received  clutter  signal  will  have  reasonable  tine,  frequency,  and  spatial  (in 
the  range  direction)  correlation  functions  (Section  5),  so  long  as  appropriate 
parameters  are  assigned  the  variables.  If  some  scatterers  are  far  enough  to 
the  side  of  the  radar  line  of  sight,  they  may  be  less  strongly  illuminated 
because  of  the  antenna  pattern.  This  effect  can  easily  be  accounted  for  by  a 
modification  of  a,*  for  each  scatterer  to  compensate  for  the  antenna  gain  in  its 

_  U  .  f 

direction.  Note  that  a must  also  contain  an  r ^  (or,  analogously,  T ^  ) 

factor  to  correct  for  the  range  to  the  scatterer;  this  correction  is  unimportant 
if  a  relatively  narrow  range  interval  far  from  the  radar  is  cf  interest  but 
becomes  very  important  if  many  ranges,  both  near  and  far  away,  are  included 
in  the  investigation. 

7.3.8  Antenna  Scanning 

Next,  suppose  the  radar  antenna  scans  in  azimuth.  As  it  scans,  new 
scatterers  enter  the  beam  and  old  ones  leave  it  (at  a  fixed  range) .  One 
straightforward  simulation  procedure  is  to  assign  angular  locations  to  scatterers 
over  the  range  to  be  covered  by  the  (azimuth)  scanning  antenna  and  to  use 
factors  related  to  the  antenna  gain  to  modify  the  scattering  coefficient  of  each 
scatterer,  depending  on  its  angular  location.  To  provide  a  reasonable 
approximation  to  the  actual  spatial  correlation  properties  of  the  signal  during 
azimuth  scan,  several  scatterers  must  occur  per  beamwidth  in  each  range- 
resolution  length.  If  the  antenna  is  to  scan  over  wide  angles,  a  much  larger 
number  of  scatterers  must  therefore  be  assigned  than  was  the  case  for  a  non¬ 
scanning  antenna.  Computer-storage  problems  may  arise  because  of  the  large 
number  of  scatterers  whose  amplitudes  and  phases  must  be  stored. 

Note  that  regardless  of  the  manner  of  representing  the  transmitted 
and  received  radar  signals  (i.e.,  whether  just  envelopes  or  the  entire  RF  wave¬ 
form  are  computed),  the  same  representation  of  the  scatterers  suffices:  a 
reflection  coefficient  and  location  roust  be  stored  for  each  scatterer.  The 
only  change  in  scatterer  information  required  by  changes  in  radar  type  is  in 
the  frequency  dependence;  if  the  simulation  is  to  be  used  for  radars  operating 
at  different  frequencies,  the  stored  reflection  coefficients  must  then  be 


functions  of  frequency.  Another  modification  might  be  desirable  in  the 
computation:  if  the  radar  had  a  large  resolution  cell,  one  might  use  fewer 
scattering  elements  (with  correspondingly  increased  amplitudes)  than  he  would 
for  a  radar  having  a  small  resolution  cell.  Again,  it  is  clear  that  this  con¬ 
sideration  affects  the  computational  procedure  used  but  not  the  clutter 
information  needed,  since  enough  scatterer  parameters  must  be  included  to 
produce  realistic  correlation  behavior  with  the  highest  resolution  radar  to  be 
investigated. 

7.4  SIMULATION  OF  POLARIZATION  CHARACTERISTICS  OF  RADAR  CLUTTER 

Thus  far  the  polarization  properties  of  the  clutter  signal  have  been 
ignored.  If  a  radar  using  ono  particular  polarization  were  always  to  be 
simulated,  the  foregoing  would  be  an  adequate  clutter  model.  However,  many 
current  and  proposed  radars  use  multiple  or  variable  polarizations,  and  for 
such  radars  the  simple  clutter  model,  wherein  a  single  reflection  coefficient 
(possibly  a  function  of  frequency)  permits  computation  of  that  portion  of 
the  received  signal  arising  at  one  scatterer,  is  not  sufficient. 

If  polarization  properties  of  radar  scatterers  (whether  target  or 
clutter  sources)  are  to  be  included  in  a  scattering  model,  the  reflection 
coefficient  for  each  scatterer  must  be  expressed  as  a  matrix  rather  than  as 
a  single  number  (Reference  19,  pp.  560-566;  Reference  21,  Appendix  E) .  Three 
amplitudes  and  two  relative  phases  (in  addition  to  another  phase  associated 
with  Xk)  are  required  to  describe  the  backscattering  from  a  target  or  clutter 
scatterer  (the  eight  elements  of  the  matrix,  four  amplitudes  and  four  phases, 
are  reduced  because  of  reciprocity  considerations  and  the  removal  of  one  phase 
angle  that  can  be  associated  with  an  effective  range  to  the  scatterer) . 
Computer  storage  requirements  are  therefore  greatly  increased  over  those  for 
simulation  when  clutter  polarization  is  neglected. 

Computation  of  waveform  envelopes,  or  of  complete  RF  waveforms,  is 
only  slightly  complicated  (although  computation  time  may  be  significantly 
increased)  by  the  requirement  that  polarization  information  be  included.  TVo 
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transmitted  signals,  one  corresponding  to  the  vertically  polarized  portion  of 

the  transmitted  wave,  and  the  other,  to  the  horizontally  polarized  portion, 

must  now  be  used  (although  in  the  special  cases  of  vertically  or  horizontally 
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polarized  transmission  one  of  them  may  be  zero) .  These  signals  will  be  in 
phase  if  linear  polarization  is  used,  but  they  will  be  out  of  phase  if  ellipti¬ 
cal  or  circular  polarization  is  used.  It  is  therefore  necessary  to  use  in-phase 
and  quadrature  components  with  a  separate  envelope  function  for  each.  Each 
component  of  the  signal  is  then  operated  on  by  the  scattering  matrix;  the 
amplitude  and  phase  of  each  component  are  thus  modified.  The  signal  at  the 
receiving  antenna  is  the  sum  of  the  signals  from  all  of  the  scatterers;  the 
summation  procedure  and  the  computation  of  envelopes,  etc.,  are  the  same  as 
were  described  previously,  except  for  the  need  to  include  in-phase  and 
quadrature  components  (which  were,  in  fact,  required  in  some  cases  even  when 
polarization  effects  were  neglected)  and  the  need  for  separate  horizontal-  and 
vertical-polarization  components. 

The  principal  difficulty  associated  with  the  polarization-inclusive 
clutter  model  is  the  assignment  of  appropriate  elements  to  the  scattering 
matrices.  Obviously  these  quantities  cannot  be  computed  from  basic  scattering 
theory  for  any  but  the  simplest  of  scatterers.  Measurements  of  scattering 
matrices  of  all  of  the  individual  scatterers  of  the  myriad  varieties  that  can 
produce  ground  clutter  are  also  not  feasible.  A  more  practical  approach  to 
the  problem  is  to  measure  scattering  matrices  of  some  typical  clutter-producing 
objects,  to  use  the  parameters  so  obtained  as  a  guide  in  assigning  coefficients 
to  the  scatterers  as  a  group  in  the  simulation,  and  then  to  compare  the  results 
obtained  from  the  simulation  with  actual  clutter  measurements.  In  this  way  it 
might  be  possible  to  obtain  realistic  clutter  simulation,  including  the  polari¬ 
zation  properties  of  the  clutter,  without  an  excessive  number  of  measurements 
being  required. 

it 

In  this  discussion,  vertical  and  horizontal  polarization  have  been  adopted  as 
reference  polarizations.  Any  orthogonal  pair  of  polarizations  might  be  chosen; 
right-  and  left-circular  polarizations  are  sometimes  a  more  convenient  choice. 
The  choice  of  a  particular  pair  to  simplify  the  discussion  is  not  intended  to 
preclude  other  choices. 
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7.5 


MASKING  EFFECTS 


In  the  discussion  so  far,  it  has  been  tacitly  assumed  that  the  radar 
directly  observes  the  clutter  elements  (and  the  target,  if  any)  at  all  points 
along  the  scattering  surface.  Although  this  condition  is  frequently  fulfilled 
for  an  airborne  radar  (though  by  no  means  always,  especially  in  mountainous 
terrain)  it  is  seldom  fulfilled  for  a  ground-based  radar.  Masking  of  clutter 
and  targets  by  clutter  can  be  produced  by  a  wide  variety  of  objects.  In  rough 
terrain,  valleys  may  be  masked  by  hills,  as  illustrated  in  Figure  7-1.  Ihe 
shaded  regions  are  not  observed  by  the  radar.  Clutter  sources  on  the  ground 
in  these  regions  do  not  affect  the  radar  because  they  are  not  observed  by  it; 
(it  is  generally  assumed  that  significant  diffraction  into  such  regions  is 
unlikely).  A  radar  target  is  equally  unobservable  in  these  regions  except  for 
multipath  and  diffraction  effects.  Masking  can  also  be  produced  by  vegetation 
or  buildings  between  the  radar  and  more  remote  clutter  sources. 


Figure  7-1  TERRAIN  MASKING 

Thus  it  is  quite  possible  that  the  radar  will  receive  clutter  returns 
from  some  range/azimuth  cells  and  not  from  others;  in  Figure  7-1  no  clutter 
returns  would  be  received  from  those  ranges  corresponding  to  the  shaded  regions. 
The  clutter  simulation  techniques  discussed  above  would  produce  clutter  returns 
at  all  ranges.  To  permit  inclusion  of  terrain  masking  effects  (not  vegetation 
or  man-made  masking),  it  is  necessary  only  to  compute,  from  topographic  infor¬ 
mation,  which  scatterers  are  masked  and  to  assign  zero  scattering  amplitude  to 
them.  Therefore,  no  change  in  the  clutter  model  or  calculations  employed  in 
the  simulation  is  required.  Computation  of  shadowed  regions  from  contour 
information  (from  a  topographic  map)  is  feasible,  and  has  been  included  in 
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digital-computer  simulations  at  CAL^22^  to  evaluate  gross  masking  effects 
(deep  shadows).  Note  that  such  computations  are  based  on  geometry  of  the 
region  involved  and  for  a  ground-based  radar  this  set  of  geometrical  param¬ 
eters  could  be  computed  once  at  the  start  of  the  computations  and  the  result  , 
stored  for  use  in  a  number  of  situations  involving  the  same  basic  terrain.  The 
discussion  in  Appendix  D  illustrates  a  method  for  evaluating  the  gross  effects 
of  masking  and  clutter. 

When  masking  is  produced  locally  by  trees  or  manmade  objects  that 
are  not  included  in  the  topographic  map,  inclusion  of  these  masking  effects 
is  more  difficult,  especially  if  masking  is  by  a  thin  stand  of  trees  that  only 
partially  masks  the  region  beyond.  In  the  latter  case,  reasonable  approximations 
based  on  selected  measurements  should  permit  realistic  simulation  of  the  masked 
clutter  although  accurate  data  would  probably  require  optical  measurements 
performed  in  the  field.  If  target -scattered  signals  are  to  be  included  in  the 
simulation,  the  same  masking  information  found  for  clutter  can  be  used  to 
indicate  when  the  target  is  being  masked. 

7.6  MULTIPATH-PROPAGATION  EFFECTS 

In  the  previous  discussion  it  has  been  assumed  that  target  and  clutter 
sources  are  independent.  This  assumption  implies  that  signals  received  from 
clutter  sources  are  unchanged  by  the  presence  of  a  target,  and  that  the  signal 
from  the  target  is  not  affected  by  the  rough  terrain  between  radar  and  target. 
These  independently  scattered  and  received  signals  are  added  at  the  radar 
receiver. 


Situations  in  which  there  is  interaction  between  terrain  and  target 
can  also  arise,  as  shown,  for  example,  in  Figure  7-2.  Note  that  the  target  in 
Figure  7-2  can  also  represent  a  clutter  scatterer  for  which  multipath  effects 
on  the  clutter  return  would  be  evaluated.  Several  possible  propagation  paths 
are  shown  here:  1)  along  the  paths  RT-TR,  2)  over  the  paths  RG-GT-TR  wherein 
the  radar  wave  is  reflected  from  the  ground  to  the  target  and  back  to  the 
radar,  3)  the  reciprocal  path  RT-TG-GR,  4)-  the  doubly  reflected  path  RG-GT-TG-GR, 
and  5)  over  the  path  RT-TS-ST-TR  where  the  bistatically  scattered  energy  from 
the  aircraft  (path  TS)  is  reflected  at  near  normal  incidence  from  the  ground 


Figure  7-2  MULTIPATH  EFFECTS 

back  to  the  aircraft,  again  suffers  bistatic  scattering  and  then  proceeds  to 
the  receiver.  For  a  low-flying  aircraft  (T),  the  differential  time  delays 
will  frequently  not  be  sufficient  to  permit  separation  in  range  of  the  echoes 
over  these  different  paths.  The  changing  interference  resulting  from  the 
changes  in  path  length  cause  the  target  signal  to  fluctuate  rapidly  even  if 
the  target  would  produce  a  steady  signal  in  the  absence  of  ground  reflections. 
Multipath  interference  also  perturbs  angle  tracking:  since  part  of  the  signal 
now  appears  to  arise  from  point  G,  the  apparent  angle  of  arrival  can  fluctuate 
over  a  much  wider  range  than  it  would  if  the  target  were  in  free  space.  Note 
that  multipath  via  the  route  RT-TS-ST-TR  arrives  at  the  receiver  from  the 
target,  so  that  this  path  by  itself  does  not  produce  large  angle-of-arrival 
errors.  It  does,  however,  cause  large  amplitude  fluctuations  through  inter¬ 
ference  with  the  direct  path,  RT-TR,  which  may  then  allow  the  TG-GR  path  to 
cause  an  increased  angle  of  arrival  error.  If  the  simulation  is  to  be  capable 
of  including  angle  of  arrival  effects,  it  will  be  necessary  to  include  angle 
information  with  the  signal  components  at  the  receiving  antenna  and  to  gather 
this  data  when  performing  experimental  programs.  For  example,  it  might  be 
necessary,  in  simulating  a  monopulse  system,  to  retain  sum  and  difference 
channels  in  the  receiver.  Angle-of-arrival  effects  would  also  be  important 
for  CON  scan,  TWS,  and  other  angle-of-arrival  measurement  techniques  --  the 
simulation  of  angle-of-arrival  (i.e.,  amplitude  and  phase  fronts)  should 
properly  be  part  of  the  "clutter"  simulation. 

For  specular  scattering  from  relatively  smooth  terrain,  there  is  only 
one  point  G,  easily  found  from  geometrical  considerations,  that  is  especially 
important.  It  is  the  point  at  which  the  ray  to  the  (inverted)  image  r'  from 
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the  radar  intersects  the  ground.  For  the  path  RT-TS-ST-TR,  it  is  the  point  S 
directly  under  the  target  that  is  important. 

For  angles  near  grazing  over  smooth  terrain,  the  reflection  coefficient 
is  very  nearly  -1;  at  larger  angles,  analytical  expressions  depending  upon  wave¬ 
length,  polarization,  angle  of  incidence,  and  physical  properties  of  the  material 
in  region  G  are  available.  It  is  thus  possible  to  include  the  effect  of  ground 
reflection  by  including  an  extra  simulated  clutter  or  target -returned  signal  in 
the  direction  RG  and  with  amplitude  and  phase  determined  by  the  reflection 
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coefficient  of  the  terrain  near  G  and  the  radar  cross  section  of  the  target. 

If  the  ground  is  gently  rolling,  two  or  more  specular-reflection 
points  analogous  to  G  may  occur.  In  this  case  the  reflection  coefficient  of 
the  terrain  may  have  to  be  modified  to  take  into  account  the  spreading  of  the 
beam  produced  by  curvature  of  the  reflecting  surface  near  the  scattering 
region.  The  simulated  target  signal  now  must  include  additional  components, 
such  as  are  described  above,  for  each  of  the  specularly  reflecting  points. 
Inclusion  of  such  effects  will  increase  computation  time  and  complexity, 
especially  if  several  targets,  rather  than  one,  are  involved,  but  no  serious 
difficulties  should  be  introduced. 

So  far  in  this  section  the  discussion  has  been  limited  to  relatively 
smooth  terrain  that  produces  specular  reflections  in  the  direction  of  the 
target.  It  is  very  possible  that  the  terrain  may  be  rough  enough  to  produce 
significant  amounts  of  diffuse  scattering  and  only  negligible  specular 
scattering;  there  can  also  be  an  intermediate  range  of  roughnesses  for  which 


For  propagation  path  RG-GT-TG-GR,  the  target  radar  cross  section  is  needed  at 
a  slightly  different  angle  from  that  used  for  the  RT-TR  propagation  path. 

For  propagation  paths  RT-TG-GR  and  RG-GT-TR,  the  bistatic  target  cross  section 
is  required.  For  small  bistatic  angles,  which  will  usually  be  required  for 
low-flying  aircraft  (for  which  clutter  problems  are  most  serious),  the 
bistatic  cross  section  can  be  related  to  the  monostatic  cross  section 
(Reference  23).  The  special  case  RT-TS-ST-TR  entails  cross  section  for  a 
90-degree  or  larger  bistatic  angle,  which  is  frequently  much  greater  than 
for  the  monostatic  case. 


significant  amounts  of  both  specular  and  diffuse  scattering  occur.  A  general 
discussion  of  specular  and  diffuse  scattering  is  available  in  Reference  24; 
for  a  discussion  of  closely  related  effects,  see  also  Reference  23.  Increasing 
the  roughness  in  the  vicinity  of  a  specular-reflection  point  reduces  the 
reflection  coefficient  in  the  specular  direction.  As  a  rough  estimate, 
specular  reflection  can  be  considered  important  when  the  effective  roughness 
parameter  5  n^-  is  less  than  unity,  where  Ah  is  the  standard  deviation  of 

the  distribution  of  heights,  y  is  the  grazing  angle,  and  *  is  the  radar  wave¬ 
length  (Reference  24,  p.  246).  In  a  simulation,  it  is  probably  reasonable  to 
include  a  factor  related  to  the  roughness  parameter  to  reduce  the  reflection 
coefficient  as  roughness  increases;  such  factors  depend  upon  the  statistical 
distribution  of  heights,  but  are  of  the  same  order  of  magnitude  for  several, 
quite  different  distributions  (Reference  24,  Figure  12.1).  Slight  roughness 
can  therefore  be  treated  in  much  the  same  manner  as  was  smooth  terrain. 

As  terrain  roughness  increases,  the  specularly  reflected  wave  becomes 
less  and  less  important  in  its  effect  on  the  radar  return  from  the  target  (or 
targets),  whereas  diffuse  scattering  becomes  more  important.  Inclusion  of 
diffuse  scattering  is  much  less  straightforward  than  was  specular  scattering, 
because  propagation  here  occurs  along  many  different  paths  with  random  phasings 
because  of  the  many  scatterers  that  now  contribute  to  the  signal  at  the  target, 
and,  on  the  return  path,  to  the  signal  at  the  receiver.  Diffuse  scattering 
normally  arises  in  a  region  (the  "glistening  surface")  that  is  centered  on  the 
specular- reflection  point  and  increases  in  size  as  the  surface  roughness 
increases.  For  small  surface  roughness,  most  of  the  diffuse  scattering  arises 
near  the  center  of  this  region,  whereas  for  large  surface  roughness  most  of  the 
diffuse  scattering  arises  near  the  ends  of  this  region  (Reference  24,  pp.  249- 
277). 


The  basic  clutter  model  described  in  earlier  sections  of  this  report 
can  be  used  to  include  multipath-propagation  effects.  The  ground-clutter 
elements  that  contribute  to  diffusely  scattered  target  illumination  can  be 
represented  by  bistatic  scattering  coefficients  analogous  to  the  monostatic 
scattering  coefficient  a ^  used  to  describe  backscattering  from  clutter  elements. 
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The  scatterer  location  relative  to  the  radar  can  still  be  given  by  7,  provided 
care  is  used  to  reinterpret  this  2-way  time  delay  when  one-way  situations  are 
being  considered.  Initially,  the  phase  associated  with  the  scatterer  can  be 
assigned  by  a  random  process.  Once  the  location  and  bistatic  scattering 
properties  of  the  scattering  element  have  been  chosen,  the  calculation  of  the 
received  signal  along  each  of  the  paths  associated  with  that  scattering 
element  is  possible,  provided  the  bistatic  cross  section  of  the  target  is  known 
or  can  be  approximated.  Since  the  bistatic  angle  at  the  scattering  element 
changes  as  the  target  location  changes,  the  scatterer  properties  must  be  known 
(or  approximated)  over  the  range  of  bistatic  angles  that  can  occur.  The 
summation  of  received  signals  along  all  of  the  paths  associated  with  many 
(simulated)  scatterers  can  now  be  used  to  represent  the  effect  of  target -clutter 
interaction  (multipath  effect) .  Direct  target  return  and  clutter  backscatter 
are  simulated  as  before  and  those  signals  added  to  the  signal  arising  from 
interaction.  Also  to  be  added  are  the  specularly  reflected  interaction  signals 
discussed  at  the  beginning  of  this  section. 

By  using  a  large  number  (of  the  order  of  the  number  used  for 
simulation  of  clutter  backscatter)  of  simulated  clutter-interaction-producing 
scatterers,  one  can  probably  achieve  quite  realistic  signal  characteristics 
even  though  highly  accurate  bistatic  information  on  clutter  scatterers  is  not 
available.  Inclusion  of  masking  effects  follows  from  the  technique  discussed 
in  Section  7.4.  Comparison  of  interaction-induced  effects  on  the  received 
signal  as  measured  in  the  field  (see  Section  8  for  a  discussion  of  measurement 
techniques)  and  as  simulated  using  the  procedure  described  here  would  be 
required  before  good  estimates  of  required  accuracies  (and  number  of  scatterers 
required  in  the  simulation)  could  be  made. 

If  polarization  properties  of  the  scattering  process  are  to  be 
included  in  the  simulation,  scattering  matrices  must  replace  scattering 
coefficients,  as  indicated  in  Section  7.3.  There  is  some  doubt  as  to  whether 
the  bistatic  scattering  matrix  is  symmetric:  it  is  symmetric  for  perfectly 
conducting  targets,  but  may  not  be  for  targets  of  other  materials  (References 
20  and  22) .  For  the  present  application  the  only  effect  would  be  to  raise 
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the  number  of  parameters  in  the  scattering  matrix  from  five  to  seven.  In  any 
event,  it  is  obvious  that  a  great  deal  of  information  must  be  stored  and  that 
forward  scattering  measurements  (see  Section  7)  are  required  for  generating 
realistic  values  for  the  scatterer  parameters. 

7.7  BISTATIC  SYSTEM  CONSIDERATIONS 


It  is  also  possible  that  a  simulator  would  be  given  the  requirement 
to  include  a  bistatic  system  simulation,  that  is,  to  be  able  to  produce  a 
realistic  clutter  (and  target)  signal  at  a  receiving  antenna  when  the  trans¬ 
mitting  antenna  is  elsewhere.  The  principal  modification  required  in  the 
clutter  simulation  is  the  need  for  bistatic  scattering  coefficients  (or 
matrices)  in  place  of  the  monostatic  scattering  coefficients  (or  matrices) 
that  were  used  above.  Similarly,  of  course,  bistatic  target  information  would 
have  to  be  used. 

Inclusion  of  masking  effects  is  slightly  more  complicated  in  this 
situation,  because  masking  can  occur  on  the  transmitter-to-scatterer  path  or 
on  the  scatterer-to-receiver  path.  A  program  capable  of  including  masking  for 
the  bistatic  radar  is  suitable  for  use  with  monostatic  radars,  of  course, 
although  the  converse  is  not  necessarily  true.  As  the  bistatic-radar  program 
must  be  more  complicated,  however,  and  will  require  more  computer  time, 
separate  monostatic  and  bistatic  programs  might  well  be  desirable. 

Extension  of  the  target-clutter-interaction  signal  simulation  dis¬ 
cussed  in  Section  7.5  to  include  bistatic  radar  is  straightforward,  since 
bistatic  scattering  coefficients  are  nei  led  to  compute  the  interaction  effects 
anyway.  Programming  of  interaction  effects  should  therefore  probably  be  done 
directly  for  the  bistatic  situation,  if  bistatic  systems  are  ever  to  be 
investigated  using  the  simulation. 
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7.8 


CONCLUSIONS  -  CLUTTER  SIMULATION  TECHNIQUES  AND  PROBLEMS 


Several  basic  forms  of  clutter-signal  simulators  have  been  discussed. 
The  most  flexible  general  type  appears  to  be  the  closed- loop  form  using  digital 
computation.  This  type  of  simulation  can  include  polarization  effects,  masking, 
clutter-target  interactions,  and  bistatic-radar  effects.  Inclusion  of  all  of 
these  would  however  require  a  complicated  computer  program.  Such  a  program 
can  probably  best  be  developed  in  stages  rather  than  as  a  single  effort.  The 
necessary  input  data  for  such  a  computer  program  must  be  obtained  from  clutter 
measurements  as  described  in  Section  10. 

Other  simulations,  including  direct-generation  digital  and  analog- 
computer  simulations,  have  also  been  considered  (Appendix  E).  These  forms 
have  various  advantages  and  disadvantages  but  in  general  seem  considerably 
less  satisfactory  than  indirect-generation  digital  simulation.  An  exception 
occurs  if  a  very  simple  simulation,  embodying  few  refinements,  is  to  be  used 
for  evaluation  of  a  basic  radar:  here  either  analog  or  digital  computation  of 
a  direct-generation  simulator  may  permit  adequate  operation  at  lower  cost. 


Section  8 

RADAR  CLUTTER  MEASUREMENTS  TECHNIQUES  AND  PROBLEMS 


8.1  GENERAL 

The  goal  of  a  clutter  measurements  program  should  be  to  provide  the 
systems  designer  with  inputs  which  will  guide  radar  design  and  allow  meaning¬ 
ful  and  quantitative  simulation  of  the  effects  of  clutter  on  radar  systems. 

The  following  sections  will  discuss  some  of  the  important  elements  related 
to  the  clutter  measurements  program. 

Observables 

The  primary  electromagnetic  clutter  observables  are  as  follows: 

1.  amplitude 

2 .  phase 

a.  spatial 

b.  temporal 

3.  polarization 

These  observables,  as  measured  by  a  radar  sensor  system  are  depen¬ 
dent  on  the  parameters  of  the  radar  (resolution  cell  size,  frequency  and 
polarization),  on  the  terrain  characteristics,  and  on  the  position  of  the  radar 
relative  to  the  terrain  (aspect  angle) .  Most  measurement  programs  in  the  past 
have  collected  data  on  1,  2b  and  3  separately,  or  1  and  2b  in  combination. 
Rarely  have  all  three  been  collected  simultaneously  to  allow  correlations  to 
be  made  among  the  variables.  Additionally,  the  instrumentation  systems  used 
did  not  provide  a  significant  parametric  variation  of  pulse  packet  size  to 
determine  its  influence  on  the  measurement.  As  a  result,  the  clutter  measure¬ 
ments  are  tied  to  a  particular  set  of  sensor  parameters  (i.e.,  pulsewidth, 
beamwidth,  and  so  on).  A  table  of  sources  of  clutter  measurement  data 
including  type  of  data  is  included  in  Appendix  F,  along  with  a  reference  list. 
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Many  of  the  past  limitations  have,  of  course,  been  brought  on  by  the 
limitations  of  the  data  acquisition  system.  On  the  other  hand,  at  this  point 
it  should  be  realized  that  a  comprehensive  measurements  program  will  require 
a  complex  data  collection  and  data  analysis  and  must  be  accompanied  by  a  cor¬ 
respondingly  complex  environmental  (physical)  sensing  and  describing  system. 

The  latter  system  which  is  meant  to  relate  the  measurements  to  the  environment 
will  require  sensors  other  than  the  measuring  radar  to  obtain  the  necessary 
data.  Although  it  may  be  unnecessary  to  employ  this  comprehensive  system  in 
all  future  measurements,  it  is  of  pressing  importance  to  gather  some  pilot 
quantities  of  data  in  this  comprehensive  and  well-correlated  manner  in  order 
to  guide  subsequent  measurements  efforts.  To  attempt  to  skimp  on  measurement 
parameters  at  this  stage  would  be  false  economy  indeed. 

Figure  8-1  is  a  very  simplified  flow  diagram  of  a  possible  overall 
approach  to  the  problem.  The  various  elements  include  the  data  sensing  and 
collection  system,  the  environmental  describing  functions,  the  various  pro¬ 
cessing  systems  which  reduce  the  data  to  levels  suitable  for  the  given  task 
and  the  correlation  system  which  relates  the  essential  data  elements  (amplitude, 
phase  and  polarization)  to  the  functions  which  describe  the  environment. 

Finally,  the  correlation  system  produces  a  series  of  graphs  which  as  a  function 
of  the  parameters  indicated  (in  the  box  labeled  "parameter  variation")  produce 
a  probability  of  obtaining  a  given  value  of  reflectivity,  spectral  bandwidth, 
clutter  angle-of-arrival  variation  (related  to  spatial  phase)  and  correlation 
functions  among  the  variables  as  a  function  of  the  terrain  area  in  question. 

The  following  paragraphs  will  discuss  the  various  elements  indicated 
in  Figure  8-1. 

8.2  ENVIRONMENT  AND  ITS  DESCRIPTION 

Typically,  one  of  the  least  defined  factors  in  past  clutter  measure¬ 
ments  programs  has  been  an  accurate  quantitative  description  of  the  environment 
being  measured.  Since  both  natural  and  cultural  targets  may  be  included,  a 
good  knowledge  of  the  distribution  of  various  types  of  scatterers  should  be 
obtained. 
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Thus,  accurate  maps  and  photographic  (stereo)  data  are  needed  to 
locate  accurately  potential  major  scatterers  or  low  return  areas  within  the 
area  to  be  measured.  If  a  high  resolution  radar  is  used  in  conjunction  with 
a  map  of  the  area  (including  terrain  slope),  then  a  correlation  between  the 
cartographic  map  and  the  radar  map  can  be  made  in  terms  of  the  measured  radar 
amplitude  and  phase  data  and  the  terrain  slope,  vegetation  and  distribution  of 
major  scatterers  or  non-scattering  regions.  The  allowed  coarseness  of  the  maps 
depends  on  whether  the  targets  are  being  viewed  for  high  altitude  and  long 
range  (typically  low  grazing  angle  and  negative  mask  angle)  or  from  ground-based 
positions  (positive  mask  angle  and  short  ranges) .  In  the  former,  local  masking 
is  typically  not  a  problem  and  thus  the  accuracy  of  information  on  vegetation 
and  other  potential  obstructions  is  not  as  critical  as  for  the  ground-based 
system.  On  the  other  hand,  for  the  ground-based  systems,  local  masking  which 
obscures  both  the  target  and  potential  clutter  sources  must  be  considered  and 
included. 

At  present  the  most  accurate  way  to  obtain  local  masking  data  is 
optically.  This  procedure  is  very  tedious  and  it  would  be  better  if  a  series 
of  PPI  photographs  could  be  generated  as  a  function  of  height  above  the  site 
(e.g.,  by  use  of  helicopter  mounted  radar)  to  get  a  measure  of  the  true  radar 
mask.  Conceptually,  the  analysis  of  such  "PPI  maps"  could  be  automated  to 
obtain  the  depth  of  mask  as  a  function  of  range  and  azimuth  angle.  This 
technique  could  speed  up  the  preparation  of  near  range  portions  of  radar  coverage 
diagrams  for  ground  based  sites,  especially  when  data  with  respect  to  effects 
of  radar  elevation  is  required.  Since  we  are  more  concerned  here  with  "on-off" 
type  data  (i.e.,  is  or  is  not  the  terrain  visible?),  the  radar  characteristics, 
except  where  diffraction  effects  occur,  are  not  extremely  important. 

8.3  SENSORS 

Because  of  the  desire  to  obtain  the  clutter  observables  and  their 
interrelationships,  any  sensor  planned  for  a  future  clutter  measurements 
program  must  be  multifunctional,  i.e.,  it  must  be  able  to  collect  amplitude, 
phase  and  polarization  data  simultaneously.  Furthermore,  radar  systems  of 
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the  future  will  work  in  environments  where  the  radar  pulse  packet  size  may 
vary  from  a  size  equivalent  to  the  target  dimensions  to  many  times  greater  than 
the  target  size.  Thus  the  measurements  radar  sensor  must  have  inherent  resolution 
to  cover  a  wide  range  of  pulse  packet  variations.  For  airborne  backscatter 
measurements,  a  synthetic  aperture  radar  may  be  a  good  candidate  because  of  its 
variable-resolution  potential.  Elevation  and  azimuth  monopulse  or  an  array  of 
receiving  antennas  will  allow  collection  of  spatial  information  (amplitude, 
phase  and  angle-of-arrival) .  If  a  multimode  antenna  feed  is  used,  polarization 
diversity  information  can  also  be  collected. 

Finally,  since  the  clutter  return  varies  with  transmitter  wavelength, 
different  carrier  frequencies  must  be  provided.  Theoretically,  the  focussed 
synthetic  aperture  radar  resolution  is  independent  of  range  and  frequency.  Thus, 
such  a  device  seems  ideal  to  determine  the  effects  of  pulse  packet  size  on  the 
clutter  characteristics  described  above. 

For  ground-based  measurements,  it  may  be  possible  to  utilize  a 
helicopter-  or  truck-mounted  synthetic  aperture  which  would  allow  very  low 
altitude  (from  essentially  ground  level  on  up)  data  collection. 

For  the  forward-scatter  measurements,  an  arrangement  of  two  aircraft 
(e.g.,  helicopters)  could  be  provided  with  one  of  the  aircraft  as  the  trans¬ 
mitter  illuminator  and  the  other  as  a  receiver.  Since  coherence  would  be 
desired  in  the  receiver  (relative  to  the  transmitter),  an  accurate  frequency 
standard  could  be  used  in  the  aircraft  and  a  locking  pulse  might  be  transmitted 
via  a  data  link  to  assure  phase  lock  between  the  transmitter  and  receiver. 
Alternatively,  stable  local  oscillators  in  conjunction  with  stable  platforms 
could  be  used  at  both  link  terminals  to  provide  relative  phase  coherency  over 
the  measurement  period  (see  Appendix  G  for  a  discussion  of  this  approach). 

With  this  arrangement,  amplitude,  phase,  angle-of-arrival  and  polarization  data 
including  forward  scatter  effects  could  be  received  by  the  moving  aircraft  of 
the  pair  and  processed  coherently  with  only  one  synthetic  aperture  receiver 
used  in  a  one-way  ranging  configuration. 
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Because  of  the  large  amount  of  data  to  be  collected  (all  ranges  and 
azimuths  for  a  given  height  of  radar) ,  the  data  must  be  collected  in  a  form 
readily  amenable  to  digital  computer  processing.  Here  again,  the  normal  use 
of  optical  and  photographic  techniques  by  synthetic  aperture  radars  allow  the 
storage  of  large  quantities  of  data  which  can  be  subsequently  read  by  flying 
spot  scanners  and,  potentially,  pattern  recognition  devices.  Furthermore, 
optical  analogues  of  various  correlation  schemes  are  normally  used  to  provide 
a  transform  from  the  spatial  to  the  frequency  domain  and  vice  versa.  With 
these  techniques,  amplitude,  phase,  polarization  and  various  correlation 
functions  could  be  ’’mapped"  for  a  given  area. 

8.4  DATA  REDUCTION 

The  first  level  of  processing  indicated  in  Figure  6  is  designed  to 
put  the  data  in  a  form  suitable  for  computer  processing  and  may  involve  an 
analog-to-digital  conversion  or  analog-to-analog  (such  as  tape  outputs  to 
optical  outputs  or  vice-versa).  At  this  stage  of  processing,  all  the  necessary 
identifying  data  should  have  been  incorporated  to  speed  second  level  correla¬ 
tion  processing. 

If  the  data  is  collected  on  film  which  might  be  the  case  if  synthetic 
aperture  radars  are  used  for  data  collection,  a  flying  spot  scanner  would  be 
used  to  recover  the  data  and  processing  time  can  be  traded  for  bandwidth. 

Because  optical  techniques  can  be  used  to  perform  spatial  to  frequency  transfor¬ 
mations,  consideration  should  be  given  to  using  these  techniques  to  save  digital 
computer  processing  times  since  joint  probability  or  correlation  function 
generation  takes  large  amounts  of  computer  time. 

A  major  problem  area  in  performing  the  second-level  correlation  pro¬ 
cessing  is  that  of  defining  a  parameter  or  group  of  parameters  which  will 
categorize  the  terrain  and  allow  the  generation  of  suitable  correlation  functions 
for  certain  terrain-describing  parameters.  Accomplishing  the  latter  would 
allow  one  to  predict  clutter  characteristics. 
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As  indicated  in  Figure  8-1,  terrain  masking  or  mask  angle  distributions, 
relief,  mean  slope,  etc.  are  all  potential  terrain  parameters  which  might 
be  used.  Probably  two  of  the  most  critical  parameters  will  be  mask  angle 
and  mean  slope.  The  former  is  important  because  it  indicates  the  lack  of 
target,  clutter,  or  both  and  is  one  of  the  primary  factors  in  the  difference 
between  an  average  clutter  cross  section  and  the  instantaneous  value  obtained 
in  practice.  Because  of  the  variability  of  terrain,  the  problem  of  obtaining 
a  single  or  a  few  generalized  descriptive  parameters  for  any  terrain  is 
exceedingly  difficult  and  is  a  definite  item  of  study  in  future  clutter  research 
and  measurements  programs. 


Section  9 

CLUTTER  PROBLEMS  -  INNOVATIONS  IN  RADAR  TECHNOLOGY 


9.1  ELECTRONIC  STEERING 

Electronically  steered  phased  arrays  have  the  capability  of  rapidly 
controlling  the  antenna  aperture  illumination  function.  This  provides  the 
ability  to  perform  rapid  scanning  among  beam  positions  through  electronic 
control  of  the  (usually  linear)  phase  taper  across  the  aperture.  In  this  manner, 
the  phased  array  can  be  programmed  to  rapidly  perform  search  and  track  proce¬ 
dures  over  the  extent  of  the  array  angular  coverage.  The  format  of  the  search 
and  track  procedures  can  also  be  varied  depending  on  the  particular  detection 
requirements.  The  primary  advantage  of  ground-based  electronically  steered 
phased  arrays  over  mechanically  steered  systems  is,  however,  the  capability 
to  perform  multiple  engagements  against  threats  from  different  directions . 
Applicable  detection/tracking  techniques  for  any  single  engagement  are  similar 
for  both  a  phased  array  and  mechanically  steered  antenna  and  thus  clutter  prob¬ 
lems,  in  this  regard,  would  be  similar  for  both  systems. 

One  form  of  electronically  steered  phased  array  uses  digital 
array  element  phase  control  which  provides  a  number  of  possible  phase  settings 
between  +180°  phase  shift.  The  instantaneous  array  bandwidth  is  governed  by 
the  aperture  size  for  a  maximum  allowable  phase  error  (which  in  turn  governs 
the  array  pattern  sidelobe  level)  across  the  array  aperture.  For  example, 
allowing  a  maximum  phase  error  of  0  results  in  an  instantaneous  fractional 
bandwidth  given  by 

Ai_  =  *  (9-D 

■f  Z  TT  q 

where  —  is  the  aperture  length  normalized  with  respect  to  wavelength  and  6  is 
the  angular  direction  of  the  beam  from  the  array  normal.  Thus,  for  an  aperture 
dimension  of  25  A  (  2.5°  beamwidth)  the  array  bandwidth  would  be  approximately 
0.7%  assuming  an  allowable  phase  error  of0=^ radians  and  a  scan  angle  of 
©=  45°.  While  this  bandwidth  is  sufficient  for  most  processing  techniques 
presently  used  for  clutter  rejection,  consideration  would  have  to  be  given  to 
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this  factor  for  techniques  such  as  frequency  hopping  which  can  be  used  to 
reduce  the  clutter  correlation  and  so  improve  target  detection. 

For  AMTI  systems  the  array  design  with  respect  to  antenna  sidelobe 
control  is  particularly  important  because  of  increases  in  the  clutter  spectrum 
bandwidth  due  to  platform  motion  and  resultant  decrease  in  clutter  suppression. 
The  requirement  of  low  antenna  sidelobe  levels  places  very  severe  requirements 
on  array  element  phase  and  amplitude  accuracy  and  this  factor,  when  compounded 
with  the  operating  environment,  can  lead  to  serious  design  problems. 

It  is  possible  to  consider  techniques  wherein  the  array  amplitude 
and  phase  aperture  distribution  are  controlled  in  order  to  maximize  the  target 
(signal)  to  clutter  ratio.  These  techniques  require  knowledge  of  the  angular 
distribution  of  the  clutter  characteristics  to  perform  the  optimization.  It 
may  be  possible  to  employ  such  optimizations  for  fixed  sites  where  detailed 
clutter  "maps"  can  be  generated  although  the  potential  advantages  must  be 
evaluated  with  respect  to  the  significant  increases  in  beam  control  and  compu¬ 
ter  complexity.  Clutter  measurements,  such  as  described  in  Section  8,  would 
be  required  to  determine  the  clutter  environment  for  a  given  location.  From 
these  measurements  a  clutter  model  for  use  in  indirect  clutter  simulation,  as 
described  in  Section  7,  could  be  obtained  and  used  to  evaluate  the  performance 
of  such  antenna  optimizations. 


9.2  SIGNAL  PROCESSING/DESIGN  AND  COMPUTER  TECHNOLOGY 


The  design  of  radar  signal  waveforms  and  data  processing  circuits 
for  clutter  reduction  purposes  requires  knowledge  of  the  signal  characteristics 
of  both  clutter  and  targets.  In  the  absence  of  detailed  knowledge  clutter 
has  conventionally  been  modeled  as  an  additive  narrow  band  Gaussian  process 
with  a  Gaussian  spectrum.  Effects  of  platform  motion,  antenna  pattern  and 
range  gate  movement  and  system  instabilities  on  the  clutter  spectrum  can  be 
included  by  well  known  relations.  In  some  cases  clutter  rejection  tech¬ 
niques  have  been  evaluated  using  a  white  noise  representation  of  the  clutter 
signal . 
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The  assumptions  described  above  concerning  the  nature  of  the  clutter 
signal  are  usually  made  on  the  basis  of  mathematical  tractability  and  are 
reasonable  for  initial  technique  investigations.  They  may,  however,  yield  quite 
incorrect  results  even  for  present-day  radars.  In  studies  of  future  clutter 
rejection  schemes,  where  increased  rejection  is  desired,  the  approach  outlined 
above  will  probably  not  be  sufficient  to  permit  a  realistic  evaluation  of 
increases  in  clutter  rejection  versus  increased  system  cost  and  complexity. 

The  increased  use  of  integrated  circuits  and  computer  technology  can 
be  foreseen  in  future  radar  systems,  especially  with  regard  to  the  use  of 
electronically  steered  phased  arrays  and  multidimensional  (amplitude,  frequency, 
polarization)  radar  detection  and  tracking.  The  detection  and  tracking  func¬ 
tions,  threat  assessment  and  allocation  of  priorities  will  probably  be 
automated  to  a  high  degree  in  these  systems  and  the  detection  procedure  may  be 
controlled  by  a  clutter/masking  map  for  the  terrain  in  the  vicinity  of  the  radar. 
It  is  evident  that  the  present  state-of -knowledge  concerning  clutter  is  in¬ 
sufficient  for  the  design  of  such  systems.  The  basic  problem  facing  the 
designer  is  that  clutter  cannot  be  represented  to  good  accuracy  in  simple 
functional  form  and  that,  even  for  a  fixed  radar  site,  the  clutter  character¬ 
istics  vary  spatially  and  with  time.  Thus,  to  realistically  evaluate  clutter 
rejection  techniques,  before  field  tests,  some  form  of  simulation  of  clutter 
is  required.  Present  simulators  using  filtered  noise  are,  however,  little 
better  (if  any)  in  evaluation  of  clutter  rejection  techniques  than  the  analytical 
analog.  In  some  cases  recordings  of  clutter  data  may  be  available  for  simulation 
purposes  from  a  radar  sensor  similar  to  that  being  considered;  generally  such 
data  are  available  only  after  system  field  tests  have  been  performed.  The 
only  solution  to  this  aspect  of  the  problem  is  to  perform  a  well-designed 
clutter  measurements  program  (see  Section  8)  and  use  the  resultant  data  in 
a  realistic  form  of  simulation  (see  Section  7) . 
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Section  10 

PROGRAM  PLAN  FOR  FUTURE  CLUTTER  RESEARCH 


The  need  for  further  research  to  determine  the  characteristics  of 
clutter  (principally  ground  clutter)  is  unquestioned  because  the  parameters 
which  enter  into  the  physical  process  of  clutter  signal  generation  are  not 
well  understood.  Thus,  there  are  many  undefined  relationships  which  exist 
among  the  physical  and  electromagnetic  characteristics  of  terrain.  In  general, 
there  is  a  need  for  quantifying  relationships  among  terrain  parameters  (specif¬ 
ically  terrain  cover  and  slope)  and  the  explicit  radar  variables -frequency, 
amplitude,  phase,  polarization  and  pulse  packet  size.  Furthermore,  there  is 
a  need  to  develop  techniques  for  using  clutter  data  in  the  simulation  of 
clutter  signals  and  clutter-target  interactions.  Thus,  the  objective  of  the 
clutter  research  program  plan  described  in  this  section  is  twofold:  a)  to 
provide  a  better  understanding  of  the  underlying  characteristics  of  clutter 
with  specific  application  to  the  air  defense  systems  problems,  and  b)  to 
provide  the  basic  clutter  data  and  techniques  for  employing  the  data  in  the 
design  of  simulators  and  realistic  evaluation  of  radar  systems. 

It  should  be  noted  that  the  development  of  new  clutter  suppression 
techniques  is  dependent  largely  on  obtaining  a  better  understanding  of  the 
nature  of  clutter,  that  is,  the  scattering  sources  which  combine  to  produce 
clutter  and  the  characteristics  of  their  spatial  and  temporal  variations. 

There  are  four  major  areas  of  development  which  appear  to  require 
significant  efforts.  A  general  listing  follows  and  a  discussion  of  possible 
approaches  follows  the  general  listing. 

1.  Data  Acquisition 

Develop  ,  with  associated  design  configuration  study, 
a  system(s)  for  the  field  collection  and  recording 
of  instantaneous  clutter  data  (both  backscatter  and 
forward  scatter)  from  various  terrains  in  a  manner  to 
allow  its  re-creation  for  data  reduction  and  analysis. 
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2.  Environmental  Description 

Develop  techniques  for  relating  clutter  measurements 
data  to  geomorphological  data  (i.e.,  terrain  type, 
relief,  slope,  etc.)  to  provide  clutter-signal  predic¬ 
tions  based  on  visual  inspection  of  the  terrain  for  a 
wide  variety  of  clutter  (cultural  and  noncultural)  to 
be  used  in  the  determination  of  defense  systems 
requirements  and  performance  evaluations. 

3.  Data  Reduction 

Develop  relationships  defining  the  probability  functions 
and  spectra  of  clutter  signals  and  the  correlation 
functions  of  the  clutter  signal  observables  as  functions 
of  electromagnetic  (frequency,  polarization,  phase,  pulse 
packet  size)  and  physical  (terrain  geometry,  number  of 
scatterers,  spatial  distribution,  and  so  on)  character¬ 
istics  .  Perform  data  interpretation  studies  to  associate 
clutter  measurements  with  scattering  processes. 

4.  Simulation 

Development  of  techniques  for  realistically  simulating 
clutter  in  the  design  and  evaluation  of  defense  systems. 

10.1  DATA  ACQUISITION 

In  future  clutter  measurements  program,  clutter  data  should  be 
collected  to  permit  the  correlation  among  the  various  radar  dimensions  of 
primary  interest,  i.e.,  frequency,  polarization,  amplitude,  phase  (spatial 
and  temporal) .  In  addition,  implicit  or  derived  dimensions  such  as  angle-of- 
arrival  (which  is  pertinent  to  the  target  location  and  tracking  problem) 
should  also  be  determined  and  correlated  with  the  above  dimensions. 
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Because  we  are  dealing  with  a  multivariable  problem  and  furthermore, 
because  advanced  radar  techniques  will  certainly  involve  multivariable  radar 
signal  processing,  information  on  the  interactions  among  the  variables  is 
necessary  to  define  the  behavior  of  clutter  and  clutter  plus  target.  Sufficient 
variation  in  the  radar  parameters  must  be  provided  to  determine  its  influence 
on  the  interactions  among  variables. 

One  of  the  most  complete  clutter  measurements  facilities  in  present 

use  is  the  four -frequency  (P,  L,  C  and  X)  airborne  radar  with  dual  linear 

(24) 

polarization  capabilities  available  at  NRL.  The  coherent  system  output 

could  theoretically  be  used  with  a  suitable  processor,  optical  or  otherwise, 
to  obtain  high  resolution  synthetic  aperture  data  provided  suitable  frequency 
stability  and  motion  compensation  capability  were  available.  As  it  is,  the 
pulse  packet,  which  can  be  varied  by  4:1  (0.2SyCi  sec  to  1.0,*  sec),  provides  a 
significant  variable  resolution  capability.  Clutter  conditions  should,  however, 
be  investigated  where  the  available  system  azimuth  resolution  (effective) 
approaches  the  target  dimensions  at  ranges  up  to  20  n.mi.  (see  Section  2). 

This  could  be  accomplished  by  using  the  synthetic  aperture  approach  mentioned 
above  or  by  performing  the  clutter  measurements  at  ranges  short  enough  to 
provide  the  desired  resolution.  Although  the  NRL  four- frequency  radar  is  the 
closest  existing  approximation  to  what  might  be  considered  an  ideal  system  for 
airborne  clutter  measurements,  the  system  is  limited  in  its  ability  to  collect 
data  which  include  simultaneous  measurements  of  the  variables  (amplitude, 
phase,  polarization,  and  angle-of-arrival  or  spatial  phase  and  amplitude)  as 
a  function  of  frequency,  polarization  and  pulse  packet  size  or  resolution. 

The  measurement  of  clutter  characteristics  at  or  near  the  ground 
surface  would  require  installing  the  radar  on  a  tower  (perhaps  mounted  on  a 
truck)  or  helicopter. 
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10.1.1 


Clutter  Measurements  Radar  Sensor 


Ideally,  a  clutter  measurements  radar  sensor  would  have  the 
following  characteristics: 

a.  Frequency  diversity,  probably  UHF  through  Ka  bands 
with  interpulse  and  intrapulse  diversity  (e.g.,  10% 
bandwidth)  to  measure  frequency  dependent  effects. 

b.  Polarization  diversity  (on  both  transmission  and 
reception) . 

c.  Spatial  amplitude  and  phase  measurements  over  an 
aperture . 

d.  Capability  for  both  coherent  and  non- coherent 
signal  processing. 

e.  Variable  resolution  or  pulse  packet  size. 

It  is  desirable  that  a  clutter  measurements  radar  have  a  resolution 
capability  approaching  the  target  dimensions  at  ranges  up  to  approximately 
20  n.mi.  This  would  require  a  pulse  width  of  20  nsec  and  a  beamwidth  on  the 
order  of  1  mr.  The  ability  to  decrease  the  resolution  (longer  pulse,  larger 
beamwidth)  by  a  factor  of,  say,  500  would  allow  collection  of  clutter  data 
over  a  wide  range  of  radar  parameters . 

The  primary  difference  between  the  clutter  measurements  radar 
described  here  and  those  previously  used  to  obtain  clutter  data  is  that 
several  clutter  observables  can  be  measured  and  recorded  simultaneously.  In 
this  mannei  the  correlations  among  the  observables  can  be  extracted  from  the 
data. 


Development  of  techniques,  such  as  described  in  Section  8,  are 
required  to  obtain  basic  data  on  the  characteristics  of  the  forward  scattered 
or  multipath  signals.  The  basic  problem  here  is  that  the  surface  area  on  the 
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ground,  effective  in  forward-scattering,  is  generally  extremely  large  and  the 
categorization  of  the  environment  becomes  difficult.  It  is  probable  that  the 
same  instrument  as  described  previously  in  this  section  for  backscatter 
measurements,  can  be  used  to  perform  the  forward  scattering  measurements.  A 
remote  illuminating  antenna  excited  with  a  stable  transmitter  would  provide  a 
system  suitable  to  perform  short  duration  coherent  measurements  (see  Section  8) . 

One  means  of  obtaining  the  high  resolution  described  above  while 
retaining  the  flexibility  of  variable  resolution  is  through  the  use  of  a 
coherent  synthetic  aperture  radar  employing  pulse  compression.  The  radar 
could  be  mounted  on  a  truck,  or  tower  for  measurements  near  the  ground,  or  in 
a  helicopter  or  aircraft  for  higher  altitude  clutter  measurements.  It  is  not 
believed  that  new  techniques  or  state-of-art  advancements  are  required  to 
develop  the  clutter  measurements  radar  with  the  suggested  capabilities. 

Before  it  can  be  recommended  however  that  a  radar  as  described  above 
be  developed  for  clutter  measurements  the  specific  studies,  discussed  in  the 
following  paragraphs,  would  be  required. 

10.1.2  Required  Studies 

A  study  should  be  performed  on  the  significant  differences  between 
clutter  data  obtained  from  a  variable  aperture  or  pulse  length  (real  resolution) 
radar  versus  that  obtained  from  a  synthetic  aperture  pulse  compression  radar. 
This  recommendation  is  made  because  although  the  systems  may  have  the  same 
effective  resolution  on  discrete  targets ,  the  effect  of  spatial  and  temporal 
clutter  signal  correlations  may  be  different,  especially  when  one  considers 
clutter  with  a  few  predominant  scatterers  among  many  smaller  ones  or  cases  with 
significant  multipath  present  in  the  path  between  the  radar  and  clutter. 

A  study  should  also  be  performed  on  antenna  techniques  to  determine 
the  spatial  amplitude  and  phase  characteristics  of  the  clutter  signal.  Infor¬ 
mation  derived  from  these  measurements  should  be  capable  of  determining  angle 
of  arrival  clutter  effects  on  monopulse,  phased  array,  con-scan,  beam-lobing  or 
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TWS  radars.  This  information  could  possibly  be  obtained  by  using  a  planar 
array  of  receiving  antennas,  two  crossed  linear  arrays  with  vertical  and 
horizontal  orientation  or  with  a  four-beam  monopulse  cluster.  Considerations 
such  as  effects  of  spatial  correlation  and  receiver  noise,  measurement  accu¬ 
racy  and  relative  cost/ complexity  should  be  studied  for  various  antenna 
configurations . 

10.1.3  Recommendations  for  Initial  Data  Acquisition  Radar 

The  radar  described  in  Section  10.1.1  would  require  extensive  devel¬ 
opment  and  consideration  should  be  given  to  a  radar  system  with  somewhat  more 
limited  capability  in  the  first  phase  of  the  measurements  program  to  allow 
ground  and  flight  tests  to  be  performed  at  an  earlier  data.  Hie  results  of 
the  studies  described  in  Section  10.1.2  are  required  before  detailed  specif¬ 
ications  of  the  initial  data  acquisition  radar  can  be  made.  The  characteristics 
described  below  are  representative  of  the  capabilities  which  would  probably  be 
required. 


The  clutter  measurements  (initial)  radar  would  be  suitable  for  truck 
or  tower  measurements  and  adaptable  to  airborne  operations  with  the  following 
basic  characteristics: 

a.  Frequency  coverage  L  band  through  Ku  band  preferably 
with  separate  RF  heads  through  a  common  IF  system 

b.  Polarization  diversity  -  either  linear  or  circular 
(but  not  simultaneously)  on  transmit  and  receive 

c.  Aperture  amplitude  and  phase  measurement  capability 

d.  Both  coherent  and  non-coherent  signal  processing 
(possibly  including  synthetic  aperture  data  collection) 

e.  Variable  pulse  length  (possibly  through  pulse 
compression)  with  a  20  ns  effective  pulse  length 
lower  limit 
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f.  Provisions  for  remote  phase  lock  or  stability  between 
main  radar  data  acquisition  device  and  remote  trans¬ 
mitter  for  forward  scatter  measurements  to  determine 
multipath  effects. 

10.2  ENVIRONMENTAL  DESCRIPTION 

An  accurate  quantitative  description  of  the  environment  in  which  the 
clutter  is  measured  will  be  required  to  allow  correlation  of  the  clutter  data 
with  the  environmental  factors.  Detailed  site  investigations  in  conjunction 
with  high  resolution  radar  maps  will  facilitate  the  isolating  of  major  scat- 
terers  and  provide  easier  interpretation  and  correlation  of  clutter  character¬ 
istics  in  terms  of  terrain  features . 

The  use  of  averages  or  maxiraa-minima  to  describe  primary  terrain 
parameters  such  as  reflectivity  and  spectral  bandwidth  will  typically  give 
either  an  overly  optimistic  or  pessimistic  estimate  of  system  performance 
depending  on  the  terrain,  system,  or  target  characteristics  involved.  For 
example,  if  the  spatial  breakup  of  clutter  caused  by  terrain  masking  -  line  of 
sight  -  interactions  is  not  taken  into  account,  assignment  of  an  average  value 
of  terrain  reflectivity  to  be  used  in  a  given  area  will  produce  too  high  a 
value  of  clutter  for  that  particular  region  and  could  predict  poor  performance, 
where  in  fact  the  target  was  in  the  clear.  Therefore,  it  is  necessary  to 
obtain  visibility  data  to  determine  spatial  clutter  breakup  and  LOS  availability 
for  use  in  clutter  data  reduction  and  interpretation.  In  the  absence  of  accu¬ 
rate  detailed  maps,  stereophotographs  should  be  taken  to  obtain  a  reference 
"image"  which  will  provide  an  accurate  picture  of  the  distribution  and  type  of 
scatterers  in  the  data  acquisition  radar  field  of  view.  In  addition,  because 
of  the  need  for  Doppler  spectra  of  clutter  and  its  dependence  on  wind  conditions, 
wind  velocity  measurement  should  be  included  in  the  environmental  sensing  system. 
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10.2.1  Recommendation  for  Environmental  Designation 

1.  For  various  proposed  air  defense  sites  around  the  country,  obtain 
visibility  (masking)  data  using  the  most  appropriate  method,  i.e.,  optical  or 
map  (200  meter  grid  spacing)  generated.  The  former  would  be  used  primarily  to 
obtain  local  masking  data  where  the  radar  cannot  be  readily  elevated.  The 
latter  approach  could  be  used  to  describe  relatively  flat  country  for  an 
antenna  elevated  above  the  local  mask. 

2.  Collect  aerial  stereophotographs  during  clutter  measurements  to 
allow  a  reference  "image"  to  be  generated  indicating  distribution  of  cultural 
and  non-cultural  targets. 

3.  For  airborne  clutter  measurements,  as  required  to  evaluate  missile 
seeker  systems,  collect  aircraft  altitude,  position,  velocity  and  heading 
information  to  obtain  information  on  the  relative  aircraft  to  ground  position. 

10.3  DATA  REDUCTION 

The  reduction  of  clutter  and  environmental  data  could  probably  be 
performed  using  a  high  speed  digital  computer  facility  although  optical  tech¬ 
niques,  such  as  described  in  Section  8,  should  also  be  studied  to  determine  if 
sufficient  flexibility  is  achievable  to  perform  the  many  types  of  data  proces¬ 
sing  which  would  be  desired. 

The  amplitude  and  phase  imaging  collected  by  the  data  acquisition 
system  could  be  read  out  with  a  flying  spot  scanner.  Pulse  height  (amplitude) 
and  spectrum  analysis  would  yield  range-azimuth  "maps"  (for  a  given  area). 

These  would  present  probability  distributions  of  clutter  amplitudes  and  power 
spectra  vs.  polarization,  carrier  frequency,  pulse  packet  size  and  elevation 
angle.  Hie  major  problem,  as  discussed  in  Section  8,  is  the  definition  of 
terrain  variables  for  the  correlation  of  clutter  data.  Simple  descriptions  of 
terrain  cover,  such  as  "hardwood  forest"  are  not  sufficient.  Further  delin¬ 
eations,  such  as  season  are  also  not  entirely  sufficient;  for  example,  a  recent 
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rainfall  can  significantly  affect  clutter  characteristics.  Certainly  factors 
such  as  the  ones  given  above,  and  in  addition,  terrain  masking,  relief,  mean 
slope,  and  the  density  and  characteristics  of  major  isolated  scatterers  will 
be  primary  variables.  Further,  studies  must  be  made  to  determine  terrain 
describing  functions  suitable  for  predicting  clutter  characteristics. 

It  is  recommended  that  data  interpretation  studies  be  performed  in 
conjunction  with  the  clutter  data  reduction  to  associate  the  observed  clutter 
signal  characteristics  with  scattering  processes.  These  analytical  studies  may 
include  laboratory  models  and  experiments  to  perform  investigations,  under 
controlled  conditions,  of  scattering  phenomenon  as  related  to  clutter  charac¬ 
teristics  . 

10.4  SIMULATION 

It  is  recommended  that  simulation  of  clutter  data  for  future  air 
defense  system  design  and  evaluation  be  along  the  guidelines  of  the  closed- 
loop  clutter  model  described  in  Section  7  .  Initial  investigations  of  the 
clutter  simulator  would  be  directed  to  comparison  of  simulated  and  actual 
clutter  characteristics  as  functions  of  pulse  packet  size,  and  terrain  type. 
The  required  clutter  parameters  for  implementation  of  the  model  should  be 
available  from  the  measurements  program. 

10.5  PROGRAM  PLAN  SUMMARY 

The  following  list  is  a  summary  of  the  suggested  program  plan  for 
future  clutter  research.  The  appropriate  section  numbers  are  indicated  for 
reference  purposes. 

1.  Study  effects  on  measured  clutter  data  as  obtained  by 
a  synthetic  aperture,  pulse  compression  radar  versus 
a  real  antenna  aperture  short  pulse  radar  (Section  10.1.2). 
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2.  Investigate  antenna  techniques  for  obtaining  spatial  amplitude 
and  phase  clutter  data  (Section  10.1.2). 

3.  Design  and  fabricate  an  initial  clutter  data  acquisition 
radar  using  the  results  of  items  1  and  2,  above 
(Section  10.1.3). 

4.  Examine  experimentally  techniques  for  determining  line-of- 
sight  availability  and  local  masking  effects  for  radar  system 
site  investigations  (Sections  2,8  and  10.2)  in  conjunction 
with  the  performance  of  clutter  measurements  using  the 
initial  data  acquisition  radar  (item  3,  above). 

5.  Study  means  of  readily  associating  terrain,  cultural  and 
man-made  objects  with  radar  clutter  data  to  aid  in  data 
interpretation  and  to  facilitate  clutter  prediction 
(Section  10.2). 

6.  Perform  an  investigation  of  environmental  descriptors 
and  associated  sensors  or  measurement  techniques  for 
various  types  of  terrain  to  provide  environmental  data 
during  the  clutter  measurements  program  (Section  10.2). 

7.  Investigate  data  storage  (magnetic  tape  and  optical)  and 
reduction  (digital  computer  and  optical)  techniques  for 
obtaining  various  probability  functions,  power  spectra 
and  (cross)  correlation  functions  of  the  clutter  and 
associated  environment  (Section  10.3). 

8.  Perform  study  and  initial  development  of  a  closed-loop 
digital  clutter  simulator  (Sections  7  and  10.4).  Compare, 
simulated  and  measured  (item  4,  above)  characteristics  of 
clutter  signals. 
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Appendix  A 

INVESTIGATION  OF  PROBABILITY  DENSITY  FUNCTIONS 


A.l  PROBABILITY  DENSITY  FUNCTIONS 

The  clutter  signal  received  from  V  scatterers  can  be  written: 

5  =  *e1'  =  ZAJ  e*J  (A-l) 

J»  t 

where,  Ay  and  are  the  relative  amplitude  and  phase  of  the  signal  from 
the  J*  scatterer  and  R  is  the  envelope.  The  sum  signal  S  can  be  expressed  in 
terms  of  an  in-phase,  x  ,  and  quadrature  component,  y  ,  given  by: 


A/ 


z  -  Re.  (S)  -  R  cos  9  *  £  Aj  cos  <f>j 

(A-2) 

y  *  Jm(3)  ■  RsmG  ■  ^  Aj  sm 

(A-3) 

The  (fy  values  are  assumed  to  be  uniformly  distributed  over  2  n  radians 
and  statistically  independent.  The  A-  values  are  considered  to  be  fixed  con¬ 
stants. 


Consider  first  the  quantities  x  and  y  and  their  probability  densities. 

*  m 

Since  the  characteristic  functions  of  x  and  y  are  identical,  we  find  that  v  : 

p(Z>'t/  [  %  (*i  "■>]  cos  nt  citr  (A-4) 

is  the  exact  expression  for  the  probability  density  function  corresponding  to 
X;  obtained  from  Equation  A-4  by  replacing  z  by  y  .  Integration  to 

* 

References  for  this  appendix  are  listed  on  the  last  page  of  Appendix  A. 
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obtain  p(x)  or  p(q)  is  not  easily  performed.  For  M*2l  and  A,  •  Az  the  result  can 
be  written  in  terms  of  a  Legendre  function  of  the  second  kind  of  order  -1/2, 
but  this  result  still  requires  computation  to  obtain  p(z)  and  in  addition  is 
only  one  very  special  case.  Consequently,  computer  programs  capable  of  com¬ 
puting  p(z)  as  a  function  of  z  for  fixed  parameters  A;  have  been  written  for 
/y*2and  ior  these  programs  are  discussed  later  in  this  Appendix  along 
with  results  of  computations  using  these  computer  programs. 


Equation  A-4  is  exact  for  any  set  of  constant  amplitudes  A ■  .  If 
N»t,  p(x)  must,  by  Central  Limit  Theorem,  become  Gaussian.  The  variance  for 
scatterers  is  AL  ;  hence  p(z)  becomes  the  Gaussian  distribution: 


Pe(x) 


(A— 5) 


in  the  limit  as  //-*■  00  . 

Suppose  At  -A  for  all  i  .  In  this  case  it  is  possible  to  write  an 
asymptotic  expression  for  p(z)  that  indicates  how  p(z)  approaches  pc(t)  for 
large  but  finite  values  of  /V  .  (Actually  the  A;  could  all  be  random  variables 
having  the  same  distribution;  for  the  present  investigation,  constant,  but  not 
necessarily  equal,  A !•  are  of  interest.)  Following  the  method  of  Reference  1, 
pp.  103-104,  we  obtain: 

(A-6) 

it  is  clear  that  for  //  large,  pgA(z)  is  a  good  approximation  for  values  of  x 
not  too  far  out  into  the  tails  of  the  distribution.  If  a  good  approximation 
is  needed  at  some  z»cr,  then  it  is  necessary  that  V  be  very  large.  Because 
of  the  very  rapid  rate  of  decrease  in  pe(z)  for  large  »  ,  errors  in  this  region 
of  the  distribution  are  normally  of  no  concern.  For  x  within  2  standard  devia¬ 
tions  of  zero,  the  correction  factor  in  Equation  A-6  is  small  even  for  values 
of  n  as  small  as  4  or  5.  There  is,  however,  reason  to  question  the  validity 
of  Equation  A-6  for  small  Af  ;  for  N=  2,  for  example,  the  correct  p(z)  and  pe/tCz) 
are  very  dissimilar  (see  Figure  A-l) . 
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Also  of  importance  is  the  probability  density  function,  p(t),  of  the 
amplitude,  R  ,  of  the  received  signal  defined  in  Equation  A-l.  The  probability 
density  function  is  stated  by  Nakagami  (Reference  2,  Equation  94)  to  be: 

p(R)  -  *  /  trJifv* )  IT Te(v'Ai)dv'  (A-7) 

Jo  *-/ 

This  result  is  given  without  proof  and  is  claimed  to  hold  without  restriction. 
Watson  (Reference  3)  proves  a  very  similar  result  in  a  discussion  of  a  two- 
dimensional  random  walk  (a  problem  completely  analogous  to  the  present  one) . 
Watson  proves  the  expression  for  the  probability  distribution  function,  and 
it  is  straightforward  to  show  that  the  derivative  of  his  result  is  exactly  the 
probability  density  function  given  in  Equation  A-7. 

For  the  special  case/V=  2,  A,  =42*4  ,  an  exact  expression  for 
can  be  written  in  a  simple  form: 

’  rr/lA*-**  ’  0  ~  *  <  2A  (A-8) 

-  O  ,  #<  O ,  R  ■>  ZA 

This  result  is  useful  in  checking  some  of  the  numerical  integrations  performed 
in  the  following  analysis. 


As  before,  the  probability  density  function  can  be  written  in  simpler 
form  if  Af»\.  In  this  case,  we  obtain  the  Rayleigh  density  function: 


p(R) 


R.  >  O 


(A— 9) 


-  o 


R  <  O 


where ; 


«  =  £  A- 


If  all  of  the  A  -  =A  ,  an  asymptotic  approximation  to  p(R)  can  be 
written;  we  have:  (Reference  1,  pp.  133-136) 
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*• m  ~  * K)  Vh  (£?-Hr*fV--] 


(A- 10) 


Beckmann^  notes  that  a  more  thorough  analysis  given  in  a  Russian  text  has  a 
factor  1/2  instead  of  the  factor  3/8  given  in  Equation  A-10  and  that  the  larger 
value  may  be  more  nearly  correct  because  of  the  inclusion  of  higher-order  terms 
in  the  Russian  work.  Derivation  of  the  next  higher  order  term  using  Beckmann's 
method  does  not  change  his  first-order  term,  however,  so  the  question  remains 
open.  It  should  also  be  noted  that  for  a  small  number  of  scatterers,  less 
than  about  fifteen,  increased  accuracy  cannot  be  obtained  through  use  of  the 
next  term  in  the  asymptotic  expansion.  As  is  shown  later  in  this  appendix,  a 
better  form  of  approximation  is  available  for  the  cases  of  interest. 

One  reason  for  questioning  the  use  of  the  asymptotic  formulations 
given  by  Equations  A-6  and  A-10  is  the  total  failure  of  pM(R)  when^  »  2: 
although  one  would  not  expect  good  agreement  in  this  case,  the  complete  differ¬ 
ence  in  the  resulting  functions  suggested  that  the  equation  might  not  be  valid 
except  for  quite  large  N  .  In  fact,  it  is  shown  later  that  these  equations 
are  relatively  accurate  for  4  (and  probably  for  N  -  3  as  well)  as  long  as 

A  .  ■  A. 

•» 


A  final  approximation  to  p(R)  has  been  evaluated  as  a  part  of  the 
work  reported  here.  This  approximation  is  the  Nakagami  m-distribution 
(Reference  3);  this  distribution  is  defined  by: 


M  ( R.m,  Cl)  - 


-3*z 

r(m)  e 


(A— 11) 


When/r?«  1  this  distribution  becomes  a  Rayleigh  distribution,  when  m  ■  1/2 
(the  smallest  allowable  value)  this  distribution  becomes  a  one-sided  Gaussian 
distribution.  As  m  increases,  the  distribution  becomes  more  and  more  sharply 
peaked.  Nakagami  states  that  for  practical  purposes  p(R)\ s  approximately  equal 
to  pM(A)  where: 
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Pa, S#}  *  M(X,o'n,Z.Jll) 


(A-12) 


where ; 


(A— 13) 


In  these  expressions  represents  the  mean-square  value  of  the  distribution 
of  the  At  (assumed  to  be  a  random  variable  having  a  distribution  given  by 
Equation  A-ll).  We  are  interested  in  constant  values  for  the  Ac  ;  hence 
and  il  •  =  A*  .  Thus  we  obtain : 

y  t  * 


o 


m 


(A— 14) 


and; 


(A— 15) 


For  4/*  2  (a  case  for  which  one  would  assume  that  the  Nakagami  m-distribution 
is  a  poor  approximation  to p(RY)  we  have: 


za; 


(A— 16) 


ZflL  =  At+At  (A-17) 

For  larger  values  of  /V  the  required  parameters  are  found  from  Equations  A-14 
and  A-15.  As  will  be  shown,  this  approximation  is  very  good  for  many  situations, 
even  when  N  is  as  small  as  4;  it  is  especially  useful  when  one  amplitude  is  much 
larger  than  the  others,  a  situation  in  which  the  Rayleigh  distribution  fails 
because  the  Central  Limit  Theorem  is  not  satisfied. 
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NUMERICAL  CALCULATION  OF  PROBABILITY  DENSITIES 


A. 2 


As  noted  in  Section  A.l,  the  exact  expressions  for  the  probability 
densities  (Equations  A-4  and  A-7)  cannot  in  general  be  written  in  closed  form. 
Examination  of  the  integrands  indicate':  relatively  slow  convergence;  further¬ 
more,  the  integrands  are  very  oscillatory,  so  that  small  increments  must  be 
used  in  a  straightforward  numerical  integration  procedure.  These  two  considera¬ 
tions  indicate  that  numerical  integration  alone  is  not  practical.  Instead,  it 
is  possible  to  make  use  of  asymptotic  expansions  of  some  of  the  Bessel  functions 
to  obtain  integrable  functions  in  the  large -argument  region  for  which  the 
approximation  is  good.  Numerical  integration  then  need  only  be  used  over  a 
limited  range  of  argument  values.  In  this  way  useful  accuracy  can  be  obtained 
with  a  minimum  of  computational  effort. 


mation: 


Inherent  in  all  of  the  CAL-deve loped  computer  programs  is  the  approxi¬ 


Ja(z) 


(A-18) 


For  x  greater  than  10  or  20  this  approximation  is  excellent,  as  can  be  seen 
from  Table  A-l.  Table  A_, 

BESSEL  FUNCTION  APPROXIMATION 


z 

J.  (*) 

10  2 

-0.249617 

-0.2S0016 

11*8 

+0.001967 

+0.00442 

19*6 

+0.1800407 

+0.180112 

21*2 

+0.002017 

+0.00212 

At  peaks  of  T0(z)  such  as  at r  =  10.2  and  x  =  19.6,  the  approximation  is  excellent. 
Agreement  is  much  poorer  near  zeros,  such  as  at  z  =  11.8  and  z  *  21.2,  but 
contributions  to  the  integral  are  less  important  in  these  regions. 

First,  consider  the  computation  of  p(x)  using  Equation  A-4  when/V=  2. 

The  integration  can  be  performed  in  two  parts,  viz: 
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p(x)  *  I,  *  It 


(A-19) 


where ; 


I,  = 


(vAt)  Jt(  trAt)  cos( vx)dtr 


(A-20) 


h  - 


J'a(v'Af)  T0(vAt)  cosi'ir#) eLv* 


(A-21) 


7f  was  found  using  numerical  integration;  the  program  uses  10-point  Gaussian 
quadratures  over  an  externally  specified  number  of  subintervals  of  irt .  The 
limit  Vi  must  be  chosen  such  that  vtA,  and  VfAt  are  both  greater  than  10  or  20 
so  that  the  asymptotic  approximations  to  the  Bessel  functions  can  be  used  to 
find  .  If  A,  or/4*is  small,  a  large  value  of  ir,  is  thus  required.  In  finding 
It  ,  it  is  then  necessary  to  break  t^up  into  enough  subintervals  so  that  the 
10-point  quadrature  is  sufficiently  accurate  within  each  subinterval.  This 
consideration  is  important  because  of  the  large  number  of  oscillations  the 
integrand  can  experience  before  v'has  reached  trf  .  This  effect  is  especially 
important  if  Af  and  4t  are  quite  different  in  value. 


Ig  can  be  approximated  closely  by: 

It  12a,  «/  A  \cos(v'(A,-Ai  ))  +Mn(vfAt+Ax))^  cos(vX)dv  (A-22) 


This  integration  was  performed  exactly  in  terms  of  sine  and  cosine  integrals. 
The  result  is: 

27rjb3r{,iWA'*A‘**). j)  -">]  - 

+  Ci  [Vr| A,-Az  +x«|]  +  Cl\vi  |/4,-A2-2j||j 


(A— 2  3) 


where ; 


SlU)  *  Cc(z)  ±  -f*- dt 
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/ 


Note  the  use  of  si(z)  =  Si(z)-rr/2  in  place  of  the  more  common  St(z)\  this  form 
is  used  because  it  simplifies  Equation  A-23  and  also  because  sifz)  is  obtained, 
along  with  Ci(z)  ,  from  a  standard  computer  subroutine. 

The  CAL  computer  program  (PR0B1)  computes  Tf  ,  (*-5),  and  yo^using 

numerical  integration  of  Equation  A-20  and  direct  numerical  evaluation  of 
Equation  A-23  along  with  Equation  A-19.  In  addition,  the  Gaussian  probability 
density  is  computed  from  Equation  A-5;  although  not  expected  to  be  accurate, 
the  Gaussian  density  function  is  useful  for  comparison  purposes. 


A  second  computer  program,  PR0B2,  was  written  to  permit  integration 
of  Equation  4  when  /V  =  4.  Again  the  range  of  integration  was  broken  up  into 
two  parts  analogous  to  Equations  A-20  and  A-21,  numerical  integration  was  used 
over  the  range  0  to  ,  and  an  exact  analytic  expression  for  the  integral  from 
v,  to  infinity  of  the  asymptotic  approximation  to  the  integrand  was  evaluated. 
Thus  the  integral  was  written  in  the  form  of  Equation  A-19  with: 


T,  -  jp  J  TofvAi  )Te(  vAt )  To  ( vAs)  To  ( v'Ao)  cosfirx  )oL * 
J  Ta  (vA, )  J,(vAt )  Jt  ( XrAj  )  J0  ( V'Ao  )  cos  ('V-x)d. 


(A -2  4) 
(A-25) 


Provided  vy  is  chosen  large  enough,  all  four  Bessel  functions  in  Equation  A-25 
can  be  replaced  by  their  asymptotic  values  (Equation  A-18)  to  form  a  new 
integral  This  integral  can  be  evaluated  analytically,  yielding: 


where ; 


^  =  4n 


,,L.  ,  [ ‘t3.( *  Bj  si(Bjv,-)\  +■ 

VmJ ITT.  ifr,  L  **-,  1  J 

-t  >J 


3,  -  A,-  Az  *  A  3  -  A++  * 
Bj  --  A,  -At  -  A  3  +*4+* 
3f  =  A,  +  AZ  -A3  -  A 4  +  Z 


(A-26) 
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B7  -- 

Af 

*  Az  +  A  3 

+  A4  +  *• 

ll 

A, 

+  A  z  +A3 

-  Af  t  X 

Bf/  - 

A, 

*■  A z  -Aj 

+  A4  +  * 

3,s  =■ 

A 

-  Ag  +  Aj 

+  A4  +* 

3,r  - 

-A, 

+  A  3  +  As 

+  As  4  x 

- 

4, 

-2x,  n 

--  f.  A.  7,  3,  //,  13,  15. 

This  computer  program  also  computes  the  Gaussian  probability  density 
that  should  approximate  p(x),  at  least  when  the  A;  are  nearly  equal,  and,  in 
addition,  the  asymptotically  corrected  Gaussian  probability  density  which  is 
expected  to  be  valid  only  when  A-  *A  for  all  four  i  . 

Computer  programs  PR0B1  and  PR0B2  thus  permit  computation  of  p(z)uhen 
there  are  2  or  4  scatterers.  Extension  to  6-  or  8-scatterer  cases  would  be 
straightforward  but  quite  complicated  because  of  the  very  large  number  of  terms 
that  would  appear  in  the  expression  for  1 note  the  increase  from  4  terms  in 
Equation  A-23  to  A-32  terms  in  Equation  A-26.  Computation  of  p(z)  when  N  is 
odd  is  not  practical  because  4+  then  is  irrational  and  cannot  be  integrated 
in  closed  form. 


Also  of  interest  is  the  probability  density  function  of  the  envelope 
(amplitude)  of  the  received  signal,  p(A) .  First,  consider  p(A)vhen  there  are 
2  scatterers  (/ V-  2  in  Equation  A-7) .  The  required  integral  is: 

p(A)-Af  JoC#v)  J0fA,v)  T0(Atir)  V'dv'  (A-27) 

Again,  asymptotic  expansions  of  Bessel  functions  for  v*  between  some  v*,  and 
infinity  will  be  used.  However,  since  p(R)  may  be  desired  for  quite  small  P, 
J0(Pv)is  retained  in  both  integrals;  asymptotic  expansions  are  used  only  for 
T0(At  tx)  and  To(Ai’f)  .  An  additional  change  over  the  method  of  integration  used 
above  was  required  here:  to  permit  analytic  integration  of  lz,  the  range  must 
be  from  0  to  infinity.  TTie  numerical  integration,  X, ,  from  0  to  2^  thus  must 
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be  modified  to  compensate  for  the  inaccuracy  of  the  asymptotic  approximations 
for  small  v' .  We  thus  have: 


r,  =  Rf'  { vT„  (v>A, )  T0(vA2  )  jsso ((At  +  A  £)v)  +■ 

+  cos  (\ At  -Az  I  v']  [  £  (Rv)dtr 


Zz  - 


rr-^AiAz 


(QS  +  QC) 


where; 


r«o 

QS  » J  Je  (Rr)  Sin  (( At  )v)aLv- 


O 

oo 

1 


R>  Af  *■  A  z 
R*  At  *  Az 

O  <  A  <  A  t  +■  Az 


r°° 

QC-J  J^(Rv)cos(\A/*-Az  \v)  dv' 


(A— 28) 


(A— 29) 


00  '  *  "  \A<~AH 

0  » 

Thus  Iz  would  give  p(R)\i  the  asymptotic  approximations  to  the  Bessel  functions 
were  exact.  The  numerical  integration  of  It  provides  a  correction  based  upon 
the  difference  between  the  actual  integrand  and  its  approximation  used  in  Jt  . 
The  derived  probability  density  is  then  given  by: 


p(R)  I,  +  Tz  (A- 30) 

As  in  the  programs  discussed  above,  it  is  necessary  to  do  the  numerical  inte¬ 
gration  of  Jf  in  a  sufficient  number  of  subintervals  to  provide  the  required 
accuracy . 


In  addition  to  the  computation  of  p(R)  discussed  above,  the  PR0B3 

* 

program  computes  the  Rayleigh  probability  density  (using  Equation  A-9  for  P*(R}). 


PR0B3  and  PR0B4  (to  be  discussed  next)  also  compute  p^  (A)  ,  but  in  their 
present  form  they  produce  a  correct  result  for  px/t(R)  only  if  all  AL  =  1. 
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When  n  -  4,  integration  of  Equation  A- 7  is  considerably  more  difficult 
because  of  the  appearance  of  a  pole  at  *  0  in  the  asymptotic  form  of  the 
integrand.  This  difficulty  is  circumvented  through  some  manipulations  of  vari¬ 
ous  terms  and  the  use  of  three,  rather  than  two,  subintegrals.  We  have: 


p(R)  =  I,  *  I*  +  I3 

I,  --  *£*%(**){  vT‘  W  Jo(A*  W  J*  (*s  *)  Jc  (A4  *•)  *• 

+  f  \^-(scn((p,  +AZ ) v-)+cos(lA,  -Afl*-))  ■  ( sm((A i+A4)v)t cos(\  Aj-A4j  W)]|  d. 


(A-31) 


V' 


(A-32) 


l2  , 


z  Jc  * 


j-  (1-cos(B,v'))+(1-cos(D2.V'))  +  sm(B37r)  i-Sin(B4  ir) 


(A-33) 


1,  -  kt[  <Lt 

J*r,  * 

where ;  _ l _ 

T  "  Trz^A~tAtA5A4 

B*  *  \A  t  Az  -  Aj  -  A+\ 

■  A f  -t-  A2  +  A3  +  A* 

B3  =  \A,  -At  |  *  A3  +  A+ 

3+  —  A3  +  A+  -  |  Af  -  At  | 

3s-  -  A ,  *■  Az  +  \a3  -  A4.  | 

si  ~  *-  Az  -  \a9-  A4\ 

By  *  |  A,  -At  |  +  |  A3  -  A4  | 

*8  \(\A,-Az\-\A3-At\) 


(A-34) 
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Integration  to  obtain  J,  again  was  performed  numerically,  using  10- 
point  Gaussian  quadratures  over  sufficient  subintervals  to  obtain  the  required 
accuracy.  Sir  integrand  values  are  not  computed  at  the  endpoints  of  an 
interval,  the  uAvAJJ.on  by  0  in  the  integrand  at  v-  0  causes  no  problem.  If  a 
numerical  integration  technique  requiring  computation  of  the  integrand  at 
*  -  0  were  used,  care  would  have  to  be  taken  to  include  the  appropriate  expres¬ 
sion  for  use  at  V*  0, 


Integration  to  obtain  I3  was  performed  in  one  or  two  parts,  depend¬ 
ing  upon  /}.  For  Rvy  less  than  10,  the  expression  was  integrated  numerically, 
using  10-point  Gaussian  quadrature,  from  Rv^to  10.  Integration  from  10  to 
infinity,  or  from  Rvy  to  infinity  when  XV, >10,  was  performed  using  an  asympto¬ 
tic  expression,  equation  11.1.29  on  page  482  of  Reference  4.  (For  all  of  the 
cases  so  far  computed,  13  has  been  negligible;  if  Is  were  to  be  computed  for 
very  small,  but  nonzero,  values  of  R  ,  significant  values  of  Z3  would  result.) 


Integration  to  obtain  Ix  requires  two  basic  integrals: 
J  —  Je(Rv)stn  (Bv)  cLi/'  1 


-Cos(3ir)^elv'3-  J 


JSfn(3)  , 

0  < 

R 

B 

B 

R 

sen  R  , 

1  < 

B 

W(i|T 

)  ,  0* 

0 

9 

1  < 

(A-35) 


(A-36) 


There  are  four  integrals  of  each  type  in  lz  .  Because  of  the  discontinuous 
behavior  of  these  integrals,  considerable  logic  was  required  in  the  computation. 
Computation  is  simplified  through  use  of  the  expression: 


cosh 


-/ 

x 


In 


(z  2zf) 


(A-37 
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The  Rayleigh  density  distribution,  pg(*),  from  Equation  A-9,  and  the 
asymptotically  corrected  Rayleigh  density  distribution,  P*a(*)  are  also  com¬ 
puted  by  the  PR0B4  program.  In  addition,  the  Nakagami  0-distribution  was  com¬ 
puted  using  Equations  A-ll  through  A-17.  Thus  it  is  possible  to  compare  the 
various  approximations  with  a  relatively  accurate  computation  of  p(P)  . 

A.  3  COMPARISON  OF  PROBABILITY  FUNCTIONS 


The  discussion  will  be  broken  into  two  parts  to  reduce  confusion. 
First,  the  behavior  of  p(x)  ,  the  probability  density  of  either  the  in-phase  or 
the  quadrature  components  of  the  received  signal,  is  discussed  for  both  2- 
scatterer  and  4-scatterer  cases.  Then  the  behavior  of  p(P)  ,  the  probability 
density  of  the  envelope  of  the  received  signal,  is  discussed  for  both  cases. 

Although  the  discussion  in  this  section  is  centered  on  radar 
scattering  from  randomly  located,  constant -magnitude  scatterers,  it  should  be 
remembered  that  the  results  are  much  more  general.  The  applicability  of  the 
Central  Limit  Theorem  when  there  are  only  2  or  4  random  variables  can  thus  be 
estimated  from  the  results  given  here.  Similarly,  any  problem  that  can  be 
made  analogous  to  a  two-dimensional  random  walk  can  be  examined  using  these 
results. 


First,  let  us  consider  p(x)  for  /V  »  2.  When  A ,  *  Az  «  1,  p(x)  is 
shown  in  Figure  A-l.  Note  the  very  rapid  drop  in  p(x)  for  x  near  2:  at  x  «  1.99 
p(X)  =  0.1601;  at*  =  2,  p(x)  «=  0.0808;  at  *  =  2.01,  p(x )*  0.0010.  Since  A*)*  0, 
2,  the  value  0.0010  found  at  *  =  2.01  can  be  considered  a  result  of  the 
numerical  approximations  used.  Also  shown  are  the  Gaussian  density  function, 
p6(t)  ,  and  the  asymptotically  modified  Gaussian  density  function,  p6A(x). 

Note  that  even  for  N  -  2  the  asymptotic  correction  is  quite  small;  it  is 
obviously  insufficient  to  give  a  reasonable  approximation  to  p(x) .  In  all  of 
the  figures  fo’r  p(x)  it  must  be  remembered  that  p(x),  pc(x)  ,  and  p60(x)axt  even 
functions  of  *  for  all  values  of//  and  Az  ;  only  the  positive- *  parts  of  the 
probability  densities  are  plotted. 
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QUADRATURE  PHASE  COMPONENT  PROBABILITY  DENSITY  FUNCTIONS  FOR 
TWO  EQUAL  AMPLITUDE  RANDOM  (UNIFORM)  PHASE  SCATTERERS 
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When  A ,  and  At  are  unequal,  the  pole  in  p  (x)  that  occurred  ati>  0 
for  A,  -  At  (see  Figure  A-l)  moves  to  X  m\A,-At  |  ;  because  of  the  symmetry  of 
p(x)  ,  a  pole  now  appears  at^«  -\Af-At\  as  well.  Several  p(K)  curves  have  been 
plotted  for  unequal  values  of  A,  and  At ;  in  each  case,  <r-  -j  (A  Z  +  A  Z )  has  been 

kept  equal  to  unity  to  permit  easier  comparison  of  the  results.  In  Figure  A-2a, 
Af  *  1.36,  At  =  0.40,  while  in  Figure  A-2b,  A,  =  1.41,  At  =  0.10.  The 
curves  shown  in  Figure  A-2  are  obviously  non-Gaussian  in  shape,  as  would  be 
expected.  Also  shown  in  Figure  A-2b  is  the  probability  density  for  A,  ■  -f2. 

At*  0.  This  limiting  case  occurs  when  the  second  scatterer  drops  to  zero 
amplitude,  leaving  only  one  randomly  phased  contributor.  When  Av*  0,  we  have: 


Pt(*) 


X*£  At 


(A -38) 


O 


A 


* 

f 


as  the  exact  probability  density  function.  Note  that  for  A,  =  •/ 2 ,  Az  •  0  the 
p,(x)  curve  is  nearly  the  same  as  the  p(x) curve  except  for  the  shift  in  pole 
location.  Thus,  for  values  of  At  between  0.1  and  0  it  should  be  easy  to 
estimate  quite  accurately  the  shape  of  p(x). 


We  have  seen  that  when  there  are  only  two,  randomly  phased  components 
pit)  bears  little  resemblance  to  a  Gaussian  probability  density  function.  Con¬ 
sider  next  the  situation  in  which  there  are  4  randomly  phased  contributors. 

In  Figure  A-3  are  shown  p(x),  pg(x)  ,  and  p6/t(x)  for*<>0  (as  noted  before,  the 
probability  densities  are  even  functions  of*).  It  can  be  seen  that  the 
probability  density  function  p(x)  is  now  very  nearly  Gaussian.  Also,  although 
the  asymptotically  corrected  Gaussian  distribution  is  evidently  a  closer 
approximation  to  p(x)  than  is  the  Gaussian  distribution,  not  that  pon(x)  does 
not  give  a  perfect  representation  of  p(z). 

Consider  some  cases  of  unequal  amplitude  scatterers.  When  A,=  Az  = 
0.8,  Aj  =  1,  and  A*  =  1.311S,  p(x)  is  still  nearly  Gaussian,  as  shown  in  Figure 
A-4.  If  two  scatterers  are  large  and  two  are  small,  the  p(x)  curve  departs 
more  radically  from  the  Gaussian  in  shape;  for  example,  in  Figure  A-5  are 
compared  p(t)  and  ps(z)  for  At  *  At  =  0.20,  Aj  =  A.  -  1.40.  When  one  scatterer 
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Figure  A-2  QUADRATURE  PHASE  PROBABILITY  DENSITY  FUNCTIONS  FOR  TWO  RANDOM 
(UNIFORM)  PHASE  SCATTERERS 
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strongly  dominates  oyer  the  other  three,  a  distinct  peaking  of  the  p(x)  curve 
occurs.  Figures  A-6  through  A-9  demonstrate  this  behavior.  In  Figure  A-8 
and  A-9  there  is  obviously  no  point  in  comparing  the  probability  density 
function  with  that  of  a  Gaussian  process.  Instead,  the  curves  for  probability 
density  if  A,  -  Az  =  Aj  =  0  and  A+  remains  at  its  original  value  are  given;  these 
curves  are  found  using  Equation  A-38  with  At  replacing  At  .  We  see  that,  as  for 
the  2-scatterer  case,  the  actual  probability  density  closely  approximates  that 
of  a  single,  randomly  phased  scatterer  when  the  other  scatterer(s)  becomes 
small  in  magnitude. 

In  Figures  A-l  through  A-3  and  A-4  through  A-9  it  has  been  shown 
that  p(x)  is  nearly  Gaussian  if  there  are  four  equal,  or  nearly  equal,  scatterers 
but  that  when  there  are  two  scatterers,  p(i)  is  non-Gaussian  for  all  sets  of 
scatterer  amplitudes  tried  so  far.  It  appears  probable  that  p(x)  will  never  be 
even  approximately  Gaussian  in  the  two-scatterer  case.  If  there  are  four 
scatterers  and  one  or  two  dominate  the  others,  a  non-Gaussian  probability 
density  results. 

Also  of  interest  in  clutter  analysis  is  the  envelope  probability 
density,  piP)  ,  that  corresponds  to  the  probability  densities  p(x)  discussed 
above.  The  situation  is  as  shown  in  Figure  A-10.  The  phasors  A, ,  Aj.,  Aj  ,  A * 
have  uniformly  distributed  phase  angles;  x  represents  the  horizontal  (or, 
equivalently,  the  vertical)  projection  of  the  sum  of  the  phasors,  while  R  is 
the  magnitude  of  the  phasor  sum.  R  thus  must  be  greater  than  (or  equal  to) 
zero,  where  as  x  can  be  positive  or  negative.  If  p(*)  and  p(y) are  Gaussian  and 
uncorrelated,  then  p(f?)  is  Rayleigh.  Thus  it  would  be  expected,  on  the  basis 
of  the  results  discussed  above,  that,  for  N  =  2,  p(rr)  will  not  be  Rayleigh, 
but,  for  N  =  4  and  scatterers  of  approximately  equal  magnitude,  p(R)  will  be 
approximately  Rayleigh.  This  expectation  is  correct,  as  will  be  shown  in  the 
following  discussion  and  figures. 
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Figure  A- 9  QUADRATURE  PHASE  PROBABILITY  DENSITY  FUNCTIONS  -  FOUR  SCATTERERS 
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Figure  A-10  ADDITION  OF  PHASORS 


First,  suppose  N  *  2  and  At  »  A2  =  1.  For  this  case.  Equation  A-8 
was  used  to  find  p(/t)  exactly.  This  equation  is  useful  for  checking  the  numer¬ 
ically  computed  values  of  p(f)  using  the  PR0B3  program  (see  Section  A. 2).  In 
Table  A-2  the  exact  plP)  is  compared  with  the  numerical  solution  for  several 
choices  of  r,  ( v,  is  the  upper  limit  in  the  numerical  integration  of  I, )  and 
m  {m  is  the  nuraber  of  sub  intervals  used  for  the  numerical  integration  I,  : 
increasing  m  ,  within  bounds,  increases  accuracy  but  also  increases  cost  of 
computations) . 

We  see  that  V)  =  10  leads  to  small  errors,  less  than  0.5  percent 
(entirely  negligible  as  far  as  a  plot  is  concerned),  and  that  »  20  leads  to 
very  small  errors;  m  *  0.25  v-  or  0.5  seems  to  be  satisfactory.  P  graph  of 
the  results  shown  in  the  first  two  columns  of  Table  A-2  is  given  in  cigure 
A- 11.  Also  shown,  for  comparison,  are  p/t(P)  ,  the  Rayleigh  density  distribu¬ 
tion,  and  pm (A),  the  asymptotically  corrected  Rayleigh  density.  The  actual 
p(p)  obviously  bears  no  resemblance  to  either  approximation,  as  was  expected 
from  the  non-Gaussian  character  of  p(x)  .  Note  that  in  this  and  the  remaining 
figures  p(P)~  0  for>f<0:  p(R) is  not  symmetric,  although  p(x)  is .  It  is  also 
obvious  from  Figure  A-9  that  p(A')-  0  for 

Next,  consider  the  case  of  A  At  •  From  Figure  A-10  we  see  that 
pip)-  0  for  A</Af-AtJ.  Suppose  A,  *  1.26  and  A2  s  0.63.  Then  p(R)  is  non¬ 
zero  only  for  0.63  <P<  1.89,  as  can  be  seen  in  Figure  A-12.  For/?<0.63, 

-3  -4 

the  PR0B3  computer  program  produces  outputs  of  the  order  of  10  or  10  ; 

these  values  are  a  result  of  the  numerical  approximations  used  and  indicate 


Figure  A-ll  ENVELOPE  PROBABILITY  DENSITY  FUNCTION  FOR  TWO  EQUAL 
AMPLITUDE  RANDOM  (UNIFORM)  PHASE  SCAT7ERERS 
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relatively  accurate  computations.  For  R>At+  A  t  the  program  always  returns 
zero.  The  probability  density  is  obviously  non-Rayleigh.  In  Figure  A-13  is 
shown  p(ft)  for  A,  =  1.36,  Aa  *  0.40;  here  P(r)  is  nonzero  only  for  0.9  <  ft  <  1.76. 
Hand  computations  of  the  appropriate  Nakagami  m-distribution  were  made  and  the 
resulting  p^ ( ft)  shown  on  Figure  A-13.  The  shape  of  fN  (ft)  distribution  is  wrong, 
but  it  at  least  peaks  in  the  region  of  p(ft)P  0  and  thus  indicates  where  finite 
p(it)  values  exist.  pt ^  (ft)  thus  could  be  used  as  a  first-order  approximation 
to  p(R)  if  the  shape  of  the  curve  were  not  critical.  In  Figurea  A-14  and  A-15 
are  shown  p(ft)  for  =  1.40,  A2  =  0.20,  and  for  Af  -  1.4.,  A2  =  0.10.  The 
trend  noted  before  continues,  because  the  range  of  possible  R  values  continues 
to  diminish. 

From  the  results  given  in  Figures  A-ll  through  A-15  we  see  that  for 
two  scatters  the  probability  density  of  is  non-Rayleigh  even  when  A  f  s  At  . 

From  the  one  case  for  which  the  m-distribution  approximation  ( (ft)  )  was 
computed  it  appears  that  p(ft)  is  better  (although  poorly)  approximated  by  the 
Nakagami  m-distribution  pN(ft)  than  by  the  Rayleigh  pR  (ft)  or  asymptatically 
corrected  Rayleigh  p  (ft)  distributions . 

Finally,  consider  p(ft)  when  H  -  4.  In  Figure  A-16  p  (r)  ,  f>n  (r)  , 
Pnt\  •  ant^  Ph  ^or  Pi  ~  ^2  ~  *3  ~  *+  *  1  are  compared.  The  asymptotically 
corrected  Rayleigh  probability  density  function ,  pRA(R)  ,  is  nearly  the  same  as 
the  Rayleigh  probability  density  function,  p  (r)  ;  both  functions  give  a 
reasonable  approximation  to  p(n)  .  The  Nakagami  m-distribution  curve, 
also  provides  a  reasonably  good  approximation  to  p  (r)  .  The  discontinuity  in 
shape  of  p(R) at  R  =  2  is  characteristic  of  p(R)  ,  as  will  be  seen  shortly.  Such 
discontinuities  in  shape  result  from  the  discontinuous  behavior  of  the  integrals 
given  in  equations  A-35  and  A-36;  consequently,  the  locations  of  discontinuities 
in  shape  can  be  predicted  from  knowledge  of  these  integrals  and  the  values  of 
Bj  ,  j  =  1,  ...,8  given  after  equation  A-34,  and  numerical  integration  can  be 
performed  at  more  closely  spaced  values  of  ft  in  the  region  of  a  discontinuity 
to  permit  the  shape  of  the  p  (ft)  function  to  be  accurately  found.  This 
technique  was  used  for  Figures  A-16  through  A-22. 

In  Figure  A-17  are  shown  p(ft) ,  p^(R),  and  p^(R)  when  Af  =  Az  *  0.8, 
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Figure  A-14  ENVELOPE  PROBABILITY  DENSITY  FUNCTION  FOR  TWO  SCATTERERS 
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A#  ■  1,  and  1.3115.  We  see  that  p(n) departs  more  from  the  Rayleigh 
approximation,  p (/\ ) ,  than  would  be  expected  on  the  basis  of  the  quite  close 
agreement  of  p(x)  and  the  Gaussian  probability  density  function  (Figure  A-10). 

On  the  basis  of  this  result  it  would  appear  that  }(*)  must  be  very  nearly 
Gaussian  if  p(/f)  is  to  be  well  approximated  by  a  Rayleigh  distribution.  For 
this  set  of  parameters,  with  all  scattering  amplitudes  of  the  same  order  of 
magnitude,  the  Nakagami  and  the  Rayleigh  distributions  are  different,  but 
either  one  is  a  reasonably  accurate  approximation  to  f>(R)  . 

In  Figure  A-18  is  given  p(R)  for  ■  Rt  ■  0.2,  4j  »  4m  ■  1.4;  pA)  for 
the  same  scattering  amplitudes  as  given  in  Figure  A-4.  A  tendency  towards 
the  two-spiked  p(4)  found  earlier  for  /V  •  2  can  be  noted  here.  For  this  situa¬ 
tion  neither  prf/O  nor  pv(4) is  a  very  good  approximation  to  p(X). 

Finally,  consider  the  situation  in  which  one  scatterer  strongly  dom¬ 
inates  three  equal  scatterers.  In  Figure  A- 19  it  can  be  seen  that  for  this 
situation  p#(R)  is  a  very  good  approximation  to  p(R)  ,  but  p/tU)  is  a  relatively 
poor  approximation.  The  Nakagami  m-distribution  iPn(X) )  much  more  accurately 
indicates  the  finite  width  of  p(fi) .  In  Figure  A-20  are  shown  p(X)  and  pM(X)  for 
a  more  strongly  dominant  scatterer.  The  Rayleigh  density  is  a  very  poor  approxi¬ 
mation  to  p(4)  for  the  cases  shown  in  Figures  A-20  through  A-22  (the  computed 
Rayleigh  approximation  would  be  the  same  in  each  of  these  figures  as  shown 
in  Figure  A-19).  The  Nakagami  m-distribution  becomes  better  as  the  strong 
scatterer  becomes  more  dominant. 

The  curves  of  Figure  A-22  shows  typical  behavior  of  the  Nakagami 
m-distribution  for  a  case  of  a  single  very  strong  scatterer  in  which p(p) is 
nonzero  over  only  a  small  part  of  the  range  of  ^  (and  for  which  the  Rayleigh 
distribution  is  an  extremely  poor  approximation) .  The  Nakagami  m-distribution 
follows  the  true  distribution  closely  (except  for  a  small  discrepancy  near  the 
central  peak)  and  drops  rapidly  outside  the  range  of  nonzero  p( ’/).  In  this 
case,  p(R)  is  zero  for  R>  2.1481  and  for  /«1. 8481.  For 4"-  2.2  and^»  1.5, 
values  closely  neighboring  these  limits,  p#(A)  has  already  dropped  to  0.0335 
and  0.9  x  10"15,  respectively.  The  drop  above  the  upper,  and  below  the  lower, 
limit  is  in  each  case  monotonic;  by  the  time  R  has  reached  2.5  in  the  upper 
regime,  p#(R)  has  dropped  to  0.528  x  10“50. 


118 


A. 4 


PROBABILITY  DENSITY  FUNCTIONS  FOR  V  -  4 


In  the  foregoing  results  were  presented  of  computations  for  two  and 
for  four  scatterers  of  like  or  different  scattering  amplitude.  As  noted  in 
this  Appendix,  extension  of  the  calcualtions  to  more  than  four  scatterers  would 
be  highly  complicated,  although  straightforward.  The  probability  distribution 
function  has  already  been  tabulated  for  6  ;V<  24  unit-amplitude  scatterers  by 
Freenwood  and  Durand  (Reference  5).  They  used  numerical  quadrature  on  the 
integrands  which,  f or  fi/t.  6  and  A ■  1,  converges  fast  enough  to  obviate  the 
necessity  for  analytic  integration  of  remainders. 

Their  results  for  the  probability  distribution  P(A)  can  be  translated 
into  our  probability  density  form  p(R)  by  using  finite  difference  differenti¬ 
ation.  In  Figure  23  is  shown  plR)  for  .V  ■  6;  p(R)  has  been  taken  from  Table  I 

of  Reference  5,  using  p(R  +0.2?)*.  - - to  approximate  the  derivative. 

Also  shown  in  Figure  A-23  are  pA(R)  ,  p/tt(A)*  and  p#(R).  It  is  seen  that  Pg_  (R) 
is  a  very  good  approximation  to  p(R )  in  this  case.  P/fA(*)  is  an  even  closer 
approximation  at  most  values  of  R;  the  curve  for  pM  (A)  cannot  be  shown,  but 
circles  indicate  some  of  the  points  that  lie  on  it.  The  Nakagami  m-distribution 
approximation,  pM(R)  ,  is  a  poorer  approximation  but  is  still  quite  accurate. 

One  important  factor  should  be  noted:  p(R)  in  Figure  A-23  lacks  any 
sharp  breaks  such  as  that  at  R  ■  2  in  Figure  A- 16.  Greenwood  and  Durand  indicate 
that  Rayleigh  showed  that  the  derivatives  of  the  distribution  function  P(R)  of 
order  j  (  2)  for /v' even  and  //- 3 )  for  a' odd  are  discontinuous  at  the 

points  N- 2k  *  k  *  0,1,2,3,111  Thus  for  /V  *  4,  the  density  function 
p(R)  -  P(R)  has  discontinuities  at  R  ■  2  and  at  R  ■  4.  Returning  to  Figure 
A-16,  it  is  seen  that  the  plot  of  p(R) is  continuous  but  has  a  discontinuity  in 
slope;  actually,  there  should  be  discontinuities  at  ^  *  2,4  but  these  do  not 
appear  because  of  slight  inaccuracies  resulting  from  the  numerical  integration. 
Using  Rayleigh's  result,  we  see  that  for  P  ■  6  or  7,  p(R)  is  continuous  but 
has  discontinuities  in  slope.  None  appear  in  Figure  A-23,  but  this  lack  may 
simply  be  a  result  of  the  lack  of  sufficient  tabulated  values  of  P(R)  to  permit 
a  sufficiently  fine-increment  plot  of  p(R) .  For  N >  7,  p(R)  and  its  slope 
are  continuous. 
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Appendix  B 

DERIVATION  OF  CORRELATION  FUNCTION 


Man/  situations  involving  stationary  or  moving  radar,  stationary  or 
moving  scatterers,  stationary  or  moving  (scanning)  range  bin,  etc.,  can  occur 
in  radar  system  analysis.  Rather  than  analyzing  several  special  cases,  a 
fairly  general  case  was  studied  and  the  special  cases  arising  from  it  are 
discussed  later. 

At  time  6  a  a  the  radar  signal  arises  from  /V,  scatterers  having 
scattering  amplitudes  c,  ,  ct  ,  ...,  c#f .  The  received  signal  is  thus  a  voltage: 

*,e'0'  -  £  Cjei¥j  (B-l) 


where  \ tj  is  the  phase  associated  with  the  scatterer.  The  coefficient  cj 
is  proportional  to  the  square  root  of  the  radar  cross  section  of  the  j** 
scatterers.  Note  that  independence  of  the  scatterers  (single  scattering)  is 
assumed.  For  most  clutter  situations  this  assumption  should  be  reasonable. 

is  the  envelope  of  the  received  signal;  in  Section  4  the  probability  density 
function  associated  with  £, ,  and  with  the  real  and  imaginary  parts  of  the 
summation,  were  computed  and  discussed.  We  are  now  interested  in  comparing 
the  signal  at  with  the  signal  at  r  .  Correlation  functions  relating 
and  are  therefore  of  interest  here.  However,  as  will  be  seen  shortly,  such 
correlation  functions  cannot  be  obtained  in  any  convenient  form,  in  general, 
and  a  slightly  different  approach  is  required. 

First,  consider  the  signal  received  at  .  From  Equation  B-l  we 

have: 

=  h,  +ZZ  £  cmcn  (B-2) 

(  m>n  /fa/  J 
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where; 


*  £  ci 

m*  f 

Phase  differences  affect  the  signal  envelope  f, .  Equation  B-2  can  be  written 
as: 

#t  •  +  Z  cmcn  cos\^Y(rm  *r„)l  I1  (B-3) 

<■  /»-/  L  A  -I  J 

where; 

X  m  radar  wavelength 

rm  •  range  of  m **  scatterer  (from  any  arbitrary  reference  plane) 

At  some  later  time  t-v  ,  a  new  envelope,  /fg  ,  will  exist.  Suppose  the  scat- 
terers  do  not  change  their  scattering  coefficients  cm  and  that  each  of  them 
has  maintained  a  velocity  component  um  along  the  radar  line  of  sight  during 
time  z  .  Then  we  have: 

Z  cmcn  (B_4) 

1  m>n  nm /  *■  *  *  J) 

if  it  is  assumed  that  none  of  the  scatterers  has  moved  out  of  the  space 
observed  by  the  radar. 

To  obtain  the  correlation  function  relating  /?,  and  l2  »  it  is 
necessary  to  find  the  ensemble  average  (over  scatterer  locations)  of  £<  and  £?. 
Because  of  the  irrational  (square  root)  nature  of  R,  and  R2  ,  such  an  average 
cannot  be  found  very  easily  for  general  ;  for  large  Nf  and  a  Rayleigh- 
distributed  envelope,  the  averages  can  be  found  using  the  properties  of  multi¬ 
variate  Gaussian  distributions.  An  alternative  approach  can  also  be  used: 
investigate  the  correlation  of  intensities  (power  levels)  rather  than  the 
correlations  of  envelopes  (voltage  levels).  Thus  we  consider  Tf  ^  ,  which 
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is  a  rational  function  rather  than  Ff  itself.  The  intensity  is  a  reasonable 
quantity  to  consider,  because  it  can  be  thought  of  as  the  voltage  obtained 
from  a  square -law  detector. 


The  probability  of  density  functions  obtained  in  Section  3  were  for 
p(A) ,  not  p(I)  .  Because  of  the  simple  functional  relationship  between  I 
and  P  it  is  relatively  straightforward  to  find  p(j)  given  p(R)  ,  although 
graphical  methods  (or  digital  computation)  might  be  required.  For  the  limit¬ 
ing  case  of  large,  and  a  resultant  Rayleigh  distribution,  given  by: 


p<*>  . 

the  probability  density,  p(I)  where  1m  rz  t 


F*0 

is  obtained  as: 


(B-5) 


p(I) 


1 ZO 


(B-6) 


and  the  Rayleigh  distribution  is  transformed  into  an  exponential  distribution. 


Before  finding  correlation  functions  of  1,  the  squared  envelope  of 
the  received  signal,  the  expressions  for  Ft  and  Fz  (Equations  B-3  and  B-4 
are  generalized)  and  the  notation  modified  slightly.  Defining: 

Tr»-  a'(£TL)i“'*‘*  (B-7) 

so  that  Tmt  the  equivalent  two-way  time  delay,  can  now  be  used  to  specify  the 
location  of  the  m**  scatterer.  It  is  also  assumed  that  the  radar  operates  at 
(radian)  frequency  when  I0  is  measured  and  at  frequency  u)z  »  ouf+Au>  when 
_T*  is  measured.  It  is  assumed  that  jus  is  small  enough  so  that  cm  can  be 
considered  constant.  Also,  it  is  possible  that  the  radar  observes  a  (partially) 
different  set  of  scatterers  when  It  is  measured.  This  effect,  which  arises 
from  the  scanning  of  the  radar  in  range,  requires  care  in  the  use  of  indices 
in  the  summations.  The  situation  as  shown  in  Figure  B-l  is  assumed.  The 
summation  for  Tt  now  runs  from  f  to  /V,  *  *  rip  ,  and  the  summantion  for  12 
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Figure  B-1  OBSERVED  REGIONS 


runs  from  t  to  -  /Va  ¥■  y.  ;  thus  and  A/z  are  the  total 

numbers  of  scatterers  observed  at  times  0  and  Z  .  We  thus  have: 


I,  *  A,  *  *Z  Zf  cm  c„  cos  \cot(Zm  -a&)] 


m>n  n*t 


(B-8) 


Xz--4z^zZ  Z  C'CB6os\a>t(T'-r$)+-^r(u.r-u.s^ 


r>3  Sm#.*t 


(B-9) 


where; 


A,*Z 


u*  1 


+*'c 


m  - 

From  these  expressions  we  can  find  correlation  functions  describing  effects  of 
scatterer  motion,  frequency  shifting,  and  change  in  position  of  range  bin.  It 
is  assumed  that  scatterers  remain  in  regions  A,  B,  C  (of  Figure  B-1)  during  the 
times  involved;  dimensions  of  these  regions  thus  are  large  relative  to  umZ  , 
in  general. 


First,  it  is  necessary  to  find  the  covariance  of  T4  and  Jz 


-  <I,Iz>  -<ItXh> 
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where  <>  indicates  ensemble  averages  over  Tm  and  um.  Since  it  is  assumed 
that  scatterer  locations  are  independent  of  one  another,  many  terms  in  the 
expanded  form  of  Equation  B-10  drop  out,  leaving: 


C”J\<COS  (?m- 5jJ> 

-<c°3[aJ,(r„-r„  )]><cc5[ii±I(u.m  +  cut(Tm-  r„)j  >J 


CB— 11) 


Note  that  this  summation  is  only  over  the  scatterers  in  the  common  region  B. 

If  there  are  no  scatterers  in  this  common  region,  or  if  the  range  bins  do  not 
overlap,  the  covariance  (and,  hence,  the  correlation)  is  zero.  This  result 
is  simply  a  consequence  of  the  assumed  independence  of  scatterer  locations; 
if  the  scatterers  had  a  high  degree  of  spatial  correlation,  nonzero  values  of 
co/(l,,Ig)  would  occur,  but  this  situation  is  not  of  interest  here.  To  obtain 
the  correlation  function,  the  covariance  must  be  normalized;  the  normalization 
factor  is  the  square  root  of  the  product  of  the  variances  of  Tf  and  Iz  ,  viz: 


D(l,.It)-2\Z  y  Z  c’i.cj l4  f <cos*[u>t(Tm-Tm)\> 


a  ..14 


-<cos  j  +u>t(Tr-  ts) J  > 


CB— 12) 


-<cos[rr^  (“,-«*)*■  *t(zr -**)]>*}*' 
so  the  correlation  function  is  : 


G(r,Aoo,N)=  Iz ^  £  G(N)G(rt&co) 

0  (lt,  I  z)  ' 


(B— 13) 
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The  correlation  function  indicated  here  gives  the  effect  of  receiving  pulses 
a  time  r  apart,  of  changing  frequency  (even  for  the  case  of  t  negligibly 
small),  and  of  changes  in  observed  sets  of  scatterers. 

Note  that  G(N)t  given  by: 


&(*)* 


2  2 

n  ■  +>' 


(B-14) 


expresses  the  effect  of  radar  scanning  which  changes  the  observed  set  of 
scatterers.  If  Nk*0  f  that  is,  if  there  are  no  common  scatterers,  G(N)  «  o  . 
If  A/^mA/cTO  and  *4  *  N  (same  set  of  scatterers  contributing  toJ,  and  lt ) , 
G(A/)xf  .  Further  discussion  of  O(N)  will  be  given  later  in  this  appendix. 

The  expressions  given  in  Equations  B-ll  and  B-12  are  quite  general 
and  can  be  applied  for  any  distributions  of  Tm  and  To  obtain  manageable 
expressions,  and,  from  them,  a  useful  form  for  G(r,  let  us  make  some 

reasonable  assumptions.  First,  assume  that  scatterers  are  uniformly  distri¬ 
buted  within  their  regions  (A,  B,  and  C)  and  that  velocities  and  positions 

★ 

and  Tm )  are  independent.  .  We  thus  have  a  probability  density: 


6  r 

-  —  < 
2 


AT 


-  O  , 


(B- 15) 


* 

But  note  that  T„  and  u„  can  be  statistically  dependent  in  Equations  B-ll  and 
B-12:  it  is  only  in  the  following  analysis  that  the  assumption  of  independence 
is  made. 
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It  must  be  remembered,  however,  that  the  appropriate  tv  may  be  different  in 
regions  A,B,  and  C,  because  we  are  not  assuming  a  homogeneous  distribution 
of  scatterers.  We  therefore  will  describe  thtp(rM)in  Equation  B-15  with  the 
AT  parameter  distinguished  by  a  subscript.  Each  of  the  ensemble  averages 
over  tm  -  2^  must  use  the  AV  parameter  appropriate  for  the  region  involved. 
For  the  averages  required  for  the  covariance,  Atg  will  be  used,  this  being 
appropriate  for  region  B  in  Figure  B-l.  For  the  averages  in  the  terms  of 
Equation  B-12  involving  rm  -  Tn  ,  the  IT  parameter  must  be  A  Vit  corresponding 
to  regions  B  and  C.  For  identically  distributed  and  independent  tm  and  Tn 
and  for  p(Tm)  given  by  Equation  B-15  we  have: 


r„ 


AT 


(AT)* 


,  -Jrs  -r„  <o 

, elsewhere 


(B- 16) 


where  in  each  case  the  appropriate  (subscripted)  value  of  AT  must  be  used. 
Since  only  ?m  -  7„  is  involved,  it  is  obviously  unimportant  what  V  value  is 
used  as  a  reference;  it  is  also  clear  that  the  reference  can  be  chosen  dif¬ 
ferently  for  each  region,  so  that  p(V„)  can  be  made  symmetrical  about  zero  for 
each  region  of  interest. 


We  now  can  use  Equation  B-16  to  simplify  the  averages  in 
Equations  B-ll  and  B-12.  Ihe  ensemble  averages  over  scatterer  locations  will 
thus  be  achieved,  and  it  will  only  be  necessary  to  average  over  the  velocity 
distribution  to  obtain  the  desired  answer.  We  can  show  that: 


<c« [aj,(rm -rn)\cos \f^ (um-«„)  +  (*f+Acv)(rm- )]  >  - 


CB- 17) 


le 


).  ~i 

;  *  1 

•  T“  — 

t 

/(AuJf  +  Aco)  &Z  | 
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as  before,  the  value  of  &r  appropriate  to  the  region  of  interest  must  be  used. 
Equation  B-ll  now  can  be  written: 


(  L(A*oAZ0\ 


seh*j 

/(2u>,  JAw)  A  t0\ 

[  z  1 

1 

( (2uJ,  +Aw)  AT0  \ 

<  2.  ) 

sen 


-2 


*/< _ _  . 

■n  I — = - 1  Stn 


^UJ,A  Z0  y-  ((c 


Ywf*Aw)Ar0\J  (  <  C'5  \rT~  (u,f"  '“»'>]  > 


(B-18) 


Note  that  effects  of  frequency  shift  are  now  separated  from  scatterer-velocity 
effects.  If  region  B  is  at  least  several  wavelengths  long,  or  if  Auj»cj, 
as  will  nearly  always  be  the  case,  the  second  and  third  terms  of  the  factor 
in  curly  brackets  can  be  ignored  relative  to  the  first;  this  approximation 
will  be  made  in  the  following  discussion.  A  similar  approximation  can  be  made 
in  the  evaluation  of  Equation  B-12;  in  this  case,  dependence  on  the  velocity 
difference  drops  out  altogether,  giving: 


0(1 I,  )  ** 
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We  thus  have: 


Gr(Z.Aw) 
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A  CO  A  Z0 


^AtoA  70y 
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If  the  Um  are  assumed  to  be  uniformly  distributed  over  the  range  u.0-=-j-  to 
*  we  obtain: 


i  f  A  oJ At 0\ 

Sen  ( 


Sen‘( 


A  '  ^AeoATgJ1  ' 


(B-21) 
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If  the  <4*  are  assumed  to  be  normally  distributed  with  mean  ue  and  standard 
deviation  05  ,  we  obtain: 


G-(rt&uj) 


a/ A  Ve  j ^ 


r*#r* 


Ira 

e* 


(B-22) 


We  have  thus  shown  that  &(z,Au>)  *•&(?)&(&*>),  very  ndfcrly,  so  long  as 
the  observed  region  is  several  wavelengths  long.  Thus,  the  correlation  func¬ 
tions  for  time  (scatterer  velocity  effect),  frequency  shift,  and  change  in 
observed  region  are  simply  separate  factors  (under  the  conditions  assumed) 
and  can  be  considered  separately.  Of  central  importance  here  is  the  fact  that 
nowhere  in  our  derivation  was  it  necessary  to  assume  a  large  number  of 
scatterers  or  a  set  of  equal -amplitude  scatterers:  the  results  derived  here 
are  valid  for  situations  involving  only  a  few  scatterers  or  strongly  dominant 
scatterers. 
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Appendix  C 

PHASE  RELATIONSHIPS  OF  ISOLATED,  FREQUENCY  INDEPENDENT  SCATTERERS 


The  analysis  performed  shows  that  so  long  as  the  scattering  ampli¬ 
tudes  and  phases  associated  with  individual  scattering  centers  are  frequency 
independent,  only  the  relative  scattering  phases  associated  with  the  carrier 
frequency  are  significant,  regardless  of  the  bandwidth  of  the  pulse.  Only  if 
the  scatterers  have  scattering  coefficients  that  vary  over  the  bandwidth  of 
the  pulse  will  the  straightforward  approach  of  adding  weighted  time-delayed 
images  of  the  radar  pulse  give  incorrect  results .  There  are  no  new  results  in 
this  appendix,  but  it  seems  appropriate  to  give  these  results  in  the  context 
of  the  present  study. 

Assume  the  radar  transmits  a  waveform: 

ee(t)*e,(t)co3t4j.t  (C-l) 

where  uJ0  is  an  RF  carrier  frequency  and  It)  is  an  envelope  function  which  is 
real  and  finite  but  is  in  no  way  restricted  as  to  length.  The  envelope 
function  has  a  Fourier  transform: 

/£» 

e,U)eluJ  dt  (C-2) 

The  Fourier  transform  of  the  transmitted  waveform  is  then: 

**(*)  «  {]”*,(*) \eiuJ'*  +  ciu>9*}eiuJ*dt  (C-3) 


or,  using  Equation  C-2: 

£t(u>)  =  j£~,(ou+u)4)  +  -j  E,  (u> -cOe)  ,  (C-4) 


a  well-known  result. 
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Next,  consider  a  set  of  /V  scatterers,  each  at  a  distance  producing 
a  delay  r4  for  propagation  from  radar  to  target  and  back  and  each  having  a 
scattering  amplitude  a.k  that  is  independent  of  frequency.  Note  that  txk  can 
be  complex  and  thus  can  include  any  phase  shift  inherent  in  the  scattering 
process;  this  phase  shift  is  assumed  here  to  be  frequency  independent. 

Suppose  the  transmitted  signal  were  ,  then  the  received  signal  would  be: 

Z  CC-SJ 

k  >/ 

from  which  we  have,  for  the  transfer  function, 

&(<*)- f  CC-6) 

k*t 

We  can  now  use  this  transfer  function  to  obtain  the  received 
spectrum;  we  have: 

E.(u>)  -  Ak  }  (C-7) 

C  I 


so  that  received  waveforms  can  now  be  found  using  the  inverse  Fourier  trans 
form. 


.(*)=  421  &k  \f  £’/(&+e4)0)eu'rke“utUa/  +  [£t(<i>-u>0)eia/rke'*lot<lu>\  (C-8) 

k-/  Zir  <■-"> 


If  the  phase  shift  were  linearly  dependent  on  frequency  the  following  results 
would  still  hold;  Zk  would,  in  that  case,  be  increased  to  include  the  phase 
shift. 


132 


which  can  also  be  written  as: 


t,«) -  {T  a.t  r-  \  A •r->* '*» » ■  r>  *  W-*) * 

k*f  Zrr  K- 

/_.  .  %  —L  (<AJm  t  "Tfje)  ~  LOJtfo  m  )  i  /  i 

Ef(cu-aJa)&  *  e  t£(c*j -«j0  ) 


(C-9) 


Some  rearrangement  and  an  obvious  change  of  variable  leads  to: 


-icu/t -T^) 


<Lou 


(C- 10) 


' ' *  ■ >  •  fpfj,  (oo)e-iw(i  ’  r‘W  j 


Both  integrals  are  the  same  and  are,  in  fact,  just  the  Fourier  inverses  of 
Equation  C-2.  Urns  we  have: 

-v 

B-k  Ct  C<J3  [^0  A"  t£)J  (C“  11) 

i*/ 


The  conditions  on  o-^  are  simply  those  that  are  required  so  that  no 
distortion  of  the  pulse  occurs  on  scattering  from  the  Ith  scatterer.  Since 
the  phase  shift  produced  by  a  fixed  spacing  is  linearly  dependent  on  fre¬ 
quency,  it  is  not  surprising  that  only  a  delay,  and  no  distortion,  occurs.  If 
amplitude  and  phase  properties  of  a  scattering  center  vary  rapidly  with  fre¬ 
quency  in  the  band  occupied  by  the  pulse  (as,  for  example,  in  the  case  of  a 
resonant  scatterer),  then  further  distortion,  not  accounted  for  by  Equa¬ 
tion  C-ll,  results.  This  type  of  scattering  could  occur  with  some  types  of 
clutter,  so  the  possibility  cannot  be  completely  ignored.  As-yet  unpublished 
work  at  CAL  on  internally  funded  Project  PURE  indicates  that  moderate  varia¬ 
tion  of  amplitude  across  the  radar  bandwidth  produces  essentially  no  effect. 
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Whether  er(t)  consists  of  a  series  of  replicas  of  the  transmitted 
pulse  or  is  a  highly  complex  waveform  depends  upon  whether  the  various  tk 
values  differ  by  more  than  the  pulse  length;  if  the  elemental  pulses  overlap, 
a  complicated  waveform  can  result. 
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Appendix  D 


The  detection  of  low  altitude  aircraft  by  ground  radars  is  significantly 
affected  by  terrain  clutter  and  masking  effects.  In  addition,  the  target 
characteristics  (for  example,  cross  section  versus  aspect  angle)  will  affect 
radar  system  performance.  During  the  performance  of  measurements  to  evaluate 
target  clutter  interactions  or  during  field  tests  to  evaluate  the  performance 
of  specific  radars,  it  is  difficult  to  interpret  resulting  data  unless  informa¬ 
tion  is  available  on  both  target  and  clutter  characteristics  and  their  inter¬ 
actions.  In  this  appendix,  some  of  the  factors  involved  in  interpreting  field 
test  data  are  described  for  a  specific  example  for  which  target,  clutter  and 
masking  effects  are  important. 

The  approach  taken  is  to  present  an  approach  for  gross  evaluation  of  the 
radar  detection  performance  of  an  MTI  radar  which  includes  terrain  clutter  and 
terrain  masking  effects.  Little  or  no  account  of  local  masking  is  in¬ 
cluded,  and  it  will  be  assumed  that  the  radar  can  be  elevated  above  the  local 
mask.  For  the  terrain  that  can  be  seen,  an  estimate  of  the  terrain  slopes 
as  a  function  of  range  provides  information  to  obtain  an  approximate  measure 
of  the  magnitude  of  clutter  competing  with  the  target.  The  latter  estimate 
can  be  obtained  by  utilizing  available  data  on  terrain  reflectivity  as  a 
function  of  incidence  angle.  This  type  of  information  is  presented  in 
Figure  D-l  for  X-band.  Data  on  tree  covered  terrain  at  X-  and  L-band  does 
not  differ  significantly  within  the  boundaries  shown  in  Figures  D-la  and  D-lb. 
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BACK-SCATTERING  COEFFICIENT  O"  (d§) 


where; 

terrain  reflectivity  (dimensionless)  ■ 

2 

07 *  normalized  terrain  cross  section  (m  /unit  area  of  surface) 

5  •  incidence  angle  of  radar  beam  with  local  horizontal  (radians) 

slope  of  terrain  with  local  horizontal  (radians) 

R  ■  range  to  terrain  patch  (m) 

6  ■  3  dB  antenna  azimuth  beamwidth  (radians) 

T  ■  radar  pulsewidth  (s) 

C  *  velocity  of  light  (m/s) 


-  AVERAGED  DATA 

AIRCRAFT 

< 

-  DATA  SPREAD 

APPROX.  70*  OF  DATA 

*x 

TERRAIN 

Figure  0-la  BACK-SCATTER  I NQ  COEFFICIENT  VERSUS  ASPECT  ANGLE;  ARIZONA 
COTTON  SEEDLINGS 
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FREQUENCY  CLUTTER  STUDY 


Figure  D-lb  NORMALIZED  RADAR  CROSS  SECTION  OF  WOODED  TERRAIN 
(LEAVES  ON  THE  TREES) 


r 

A  composite  curve  of  T'vs.  aspect  angle  can  be  produced  from  the 

curves  in  Figure  D-l  as  shown  in  Figure  D-2.  Using  this  curve,  an 

estimate  of  the  average  terrain  cross  section  can  be  obtained  for  a  given 

radar.  If  the  target  aircraft  cross  section  is  ,  then  the  target-to- 

clutter  ratio  (T/C)  can  be  obtained  by  the  relationship  T/C  ■  <TA/(Tr  . 

For  the  MTI  radar  system  to  be  effective,  the  subclutter  visibility  (SCV)  must 

2 

be  equal  to  or  greater  than  C/T,  i.e.,  to  reduce  a  1000  m  clutter  signal  to 
the  receiver  noise  level  requires  that  the  SCV  «  30  dB.  Then  depending  on  the 
other  system  parameters  (peak  power,  bandwidth,  receiver  noise  figure,  losses, 
integration,  etc.),  the  signal-to-noise  ratio  (SNR)  to  the  target  and  result¬ 
ing  probability  of  detection  can  be  estimated  using  the  curves  shown  in 
* 

Figure  D-3  . 

To  evaluate  this  approach,  a  study  was  made  of  actual  detection  per¬ 
formance  data  obtained  by  the  U.S.  Air  Force  using  a  radar  with  the  objective 
of  establishing  a  verification  of  the  performance  predicted  by  the  employed 
evaluation  analysis.  The  specifications  of  the  radar  are  shown  in  Table  D-l. 

The  detection  criteria  suggested  earlier  in  this  appendix  which  is 
used  to  evaluate  the  performance  of  the  search  is  based  on  a  comparison  of 
the  clutter  power  appearing  at  the  radar  MTI  filter  output  wich  that  of  a  par¬ 
ticular  size  target.  When  this  S/C  ratio  fell  below  0  dB  in  areas  of  high 
terrain  reflectivity,  the  system  was  said  to  clutter  limited.  In  practice, 
many  false  PPI  scope  indications  would  continuously  exist  in  these  areas 
because  of  uncancelled  cljtter  residue  power  exceeding  the  thresholds  moni¬ 
toring  the  MTI  filter  outputs.  It  is  apparent,  therefore,  that  an  assumption 
was  necessary  in  the  selection  of  the  MTI  filter  output  threshold  level,  viz., 


* 

These  curves  assume  that  the  clutter  residues  after  filtering  by  the  MTI 
system  have  the  characteristics  of  receiver  noise  (white  and  Gaussian) 
which  is  not  always  true. 

138 


! 


(«p)  *-£> 


139 


rigure 


level  corresponding  to  the  assumed  target  size.  For  this  analysis,  a  minimum 

2 

or  threshold  target  size  of  5  m  is  compatible  with  the  system  sensitivities 
which  can  be  determined  using  the  data  in  Table  D-l. 

The  U.S.  Air  Force  at  Rome  Air  Development  Center  conducted  tests 
of  the  radar  indicated  in  Table  D-l  on  a  test  course  near  Lake  Oneida,  N.Y., 
using  a  helicopter  as  the  target  vehicle.  Test  runs  incorporating  dive  and 
climb  maneuvers  were  conducted  at  altitudes  ranging  from  500  feet  to  10,000 
feet  above  the  radar  site.  A  target  detection  was  recorded  if  a  spot  was  at 
all  detectable  during  a  sweep  of  the  PPI  indicator.  Results  of  these  tests 
can  be  used  to  substantiate,  to  a  degree,  the  analysis  suggested  above  to 
evaluate  detection  capability.  Figure  D-4  shows  a  compilation  of  the  Air  Force 
data,  with  a  MTI  detection  indicated  by  a  short  vertical  mark  and  a  loss  of 
detection  by  a  small  open  circle.  Outbound  runs  are  indicated  by  marks  above 
the  flight  path  line  and  inbound  runs  by  marks  below  the  flight  path  line. 

In  order  to  apply  the  detection  criteria  to  the  above  particular 
terrain  and  target  (helicopter)  combination,  it  was  necessary  to  obtain  terrain 
slope  as  a  function  of  range  along  the  aircraft  ground  track.  In  addition, 
the  radar  cross  section  of  the  helicopter  was  estimated  by  comparing  this 
aircraft  to  a  helicopter  which  is  very  similar  in  size  and  configuration. 

Cross  section  data  on  this  A/C  is  available  for  S-band  (3.5  GHz)  along  with 
translated  values  at  L-band  (1.3  GHz).  Measurements  were  taken  on  a  1/10  scale 
model  using  Ka~band. 

The  terrain  profile  over  which  most  of  the  test  runs  were  made  is 
known  on  Figure  D-5  (heading  of  200  degrees).  It  consists  of  a  relatively 
flat  interval  out  to  a  range  of  approximately  5.0  n.mi.  followed  by  a  800  foot 
hillside  rise  terminating  at  approximately  7.5  n.mi.  A  line-of-sight  check 
from  the  radar  site  to  the  top  of  this  hill  indicates  that  terrain  clutter 
signals  cannot  be  received  beyond  7.5  n.mi.  because  of  masking.  Discussions 
with  AF  project  personnel  who  conducted  these  tests  indicates  that  the  heli¬ 
copter  may  have  strayed  from  the  desired  heading  such  that  complete  clutter 


masking  for  ranges  greater  than  7.5  n.mi.  may  not  have  been  experienced  for  a 
small  percentage  of  the  data.  Map  data  indicates  that  ten  degree  heading 
error  could  allow  some  clutter  (see  Figure  D-5)  returns  out  to  approximately 
14  n.mi.  where  complete  masking  takes  place. 

An  examination  of  the  data  of  Figure  D-4  for  the  outbound  run  at 
1000  feet  indicates  a  complete  detection  loss  during  a  climb  maneuver  from 
1000  feet  to  2000  feet,  but  not  during  the  return  run  over  the  same  terrain 
from  2000  feet  to  1500  feet.  Terrain  returns  are  absent  (masked)  in  this 
interval  making  MTI  detection  of  the  noncoherent  radar  dependent  on  a  "beat" 
signal  formed  by  the  combination  of  returns  from  fixed  and  moving  parts  of 
the  target,  i.e.,  the  fuselage  return  is  the  reference  signal  and  the  rotor 
return  is  the  variable.  Because  of  the  pitch  attitude  of  the  helicopter  under 
normal  flight  conditions,  the  aspect  angle  from  the  radar  site  to  the  aircraft 
on  the  inbound  run  is  approximately  1.5  degrees  above  the  fuselage  centerline, 
and  approximately  3.0  degrees  above  during  the  dive.  The  outbound  run  results 
in  an  approximate  aspiect  angle  of  4.5  degrees  below  tail  during  normal  flight 
and  1.5  degrees  below  tail  aspect  in  the  climb  maneuver.  It  will  be  shown 
below  that  aspect  angle  changes  during  maneuvers  can  explain  some  the 
detection  behavior  that  occurred  during  the  test. 

2 

A  system  sensitivity  check  of  the  radar  indicates  that  a  4  m 
target  can  be  detected  out  to  the  system  maximum  range  of  50  n.mi.  with  a 
detection  probability  of  90  percent  and  a  false  alarm  probability  of  10-6. 

It  is,  therefore,  assumed  that  the  system  threshold  is  approximately  set  for 
this  size  target  with  effective  S/N  improvement  being  accomplished  by  the 
integrating  action  of  the  PPI  indicator.  Clearly,  target  or  clutter  returns 
not  exceeding  this  value  at  the  MTI  output  will  not  be  detected  at  all. 

An  examination  of  the  helicopter  fuselage  cross  section  reveals  a 

2 

generally  constant  nose-on  cross  section  of  approximately  8  m  over  the  aspect 

angles  experienced  during  the  above  mentioned  1500  foot  run.  The  tail -on 

2 

cross  section  is  approximately  6.0  m  for  5  degrees  below  tail  aspect  and 
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2 

about  1.0  m  at  0  degrees  aspect.  The  cross  section  of  the  rotor  blades 

appears  to  be  somewhat  uniform  at  L-band  in  that  a  cross  section  of  approxi- 
2 

mately  5  m  has  a  high  probability  of  occurance  over  most  azimuth  aspect 

angles.  Since  the  maximum  aircraft  cross  section  variation  resulting  from  the 

summation  of  the  fixed  and  variable  (changing  relative  phase)  returns  is  equal 

to  twice  the  amplitude  of  the  smaller  return ,  it  becomes  apparent  why  no  detection 

was  observed  for  the  near  0  degree  aspect  tail-on  case.  Although  the  rotor  return 

remains  relatively  constant  for  the  aspect  angles  under  consideration,  the 

2 

fixed  fuselage  return  of  approximately  1  m  theoretically  results  in  no  more 
2 

than  a  2  m  effective  MTI  target  cross  section  for  detection  -  which  is  below 
2 

the  4  m  threshold  level.  The  actual  quantities  used  are  probably  in  error 
to  some  extent,  but  the  detection  loss  is  believed  to  be  due  to  this  reference 
signal  loss  effect.  Similar  losses  will  be  observed  on  Figure  D-4  during  a 
stair  case  type  run  from  1500  feet  to  5500  feet  while  the  aircraft  is  proceed¬ 
ing  outbound.  Aspect  angles  during  the  climb  maneuvers  are  similar  to  those 
just  examined. 

A  second  test  of  the  evaluation  detection  criteria  can  be  made  for 

the  available  data  as  the  aircraft  appears  over  that  terrain  which  results  in 

the  highest  clutter  returns.  An  examination  of  Figure  D-4  shows  that  reliable 

detection  is  maintained  in  this  region  (5.0  -  7.5  n.mi.)  for  most  of  the  runs. 

Ihe  radar  to  terrain  aspect  angle  is  about  7.0  degrees  for  this  hillside  which 

amounts  to  a  normalized  terrain  coefficient  oz  of  -24  dB  according  to  the  data 

of  Figure  D-3.  For  the  radar  beamdwidth  and  pulsewidth  and  range  a  <r0  of 

-25.0  dB  or  greater  will  cause  the  radar  to  be  cluf.er  limited  for  the  hillside. 

Although  the  criteria  predicts  marginal  performance  because  of  large  clutter 
2 

returns  for  a  5  m  target  threshold  setting,  the  fact  that  good  detection  was 

indeed  observed  can  be  explained  by  the  fact  that  the  helicopter  cross  section 

2 

is  generally  greater  than  5  m  (thereby  allowing  a  higher  threshold  setting) 
and  that  both  the  fuselage  and  rotor  returns  combine  with  the  clutter  refer¬ 
ence  signals  producing  greater  energy  in  the  nonzero  (detectable)  doppler 
spectral  regions. 
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Two  inbound  runs  (10,000  feet  and  3,000  feet)  result  in  a  detection 
loss  in  the  region  from  approximately  8.5  to  10  n.mi.  The  clutter  reference 
is  completely  absent  in  this  region  due  to  terrain  masking  but  detection  is 
evident  on  other  inbound  and  outbound  runs  over  the  same  terrain.  Consequently, 
detection  must  be  accompanied  by  a  sufficiently  strong  fuselage  return.  A 
plausible  explanation  appears  to  be  the  presence  of  nulls  in  the  fuselage 
cross  section  pattern  as  a  function  of  elevation  aspect  angle.  There  is  some 
evidence  which  indicates  that  helicopter  average  cross  section  tends  to  drop 
off  for  viewing  angles  10  to  15  degrees  below  nose-on. 

The  detection  loss  evident  at  a  range  of  approximately  5  n.mi. 

during  the  10,000  feet  inbound  run  is  attributed  to  a  decrease  in  sensitivity 

2 

(gain)  of  the  elevation  plane  CSC  antenna  pattern.  The  radar  used  in  this 

study  is  advertised  as  being  able  to  detect  targets  at  an  altitude  of 

10,000  feet  in  as  close  as  5  n.mi.  This  corresponds  to  a  maximum  elevation 

2 

angle  of  18.4  degrees.  Targets  of  nominal  cross  section  (5  m)  evidently  do 
not  result  in  sufficient  signal  amplitude  to  exceed  the  system  threshold  at 
higher  elevation  angles. 

Although  some  of  the  flight  test  data  is  difficult  to  explain  by 
means  of  the  supplied  and  derived  information,  most  of  the  detection  perform¬ 
ance  experienced  was  predicted  using  the  detection  evaluation  criteria  used 
in  this  appendix.  The  data  which  apparently  does  not  conform  to  the  criteria 
is  believed  to  be  the  result  of  the  uncertainty  associated  with  helicopter 
position  (true  clutter  amplitude)  and  the  differences  between  the  actual  radar 
cross-section  of  the  helicopter  used  in  tests  and  the  data  available  for  this 
analysis. 


147 


Appendix  E 
CLUTTER  SIMULATION 


E.l  DIGITAL-COMPUTER  OPEN-LOOP  SIMULATOR 

Open-loop  simulation  relies  upon  knowledge  of  received  signal  char¬ 
acteristics  to  permit  the  generation  of  a  realistic  clutter  signal  to  add  to 
a  similarly  generated  target-return  signal.  If  one  were  willing  to  simulate 
the  clutter  return  by  random  noise,  without  concern  for  time-,  frequency-,  and 
spatial -correlation  effects,  then  digital  simulation  of  the  clutter  signal 
would  be  easy.  Time  correlation  could  be  approximated  by  appropriate  (numerical) 
filtering  of  the  sequence  of  random  numbers  used  to  simulate  the  noise.  Fre¬ 
quency-  and  spatial-correlation  effects  would  be  much  harder  to  include  in  an 
open-loop  simulation.  Artificial  modification  of  the  noise-simulated  clutter 
signal  to  introduce  these  effects  would  probably  involve  as  much  complexity  as 
is  involved  in  the  closed-loop  method  but  with  none  of  the  advantages  (primar¬ 
ily  greatly  increased  flexibility)  of  that  method.  If  dominant  clutter 
scatterers  existed,  so  that  the  probability  density  of  the  signal  amplitude 
distribution  was  non-Rayleigh,  it  would  be  considerably  harder  to  produce  a 
simulated  clutter  signal  using  open-loop  simulation. 

Inclusion  of  polarization  effects  would  require  still  further  com¬ 
plexity  if  the  directly  generated  signals  were  to  have  appropriate  character¬ 
istics  in  the  channels  corresponding  to  the  two  orthogonally  polarized  radar 
returns.  Again  it  appears  likely  that,  except  for  a  very  crude  modeling  of 
the  clutter  characteristics,  open-loop  simulation  will  be  as  complicated  as 
closed-loop  simulation.  Similar  difficulties  arise  if  masking,  target-clutter 
interaction,  and  bistatic-radar  effects  are  to  be  included  in  the  simulation. 

Open-loop  simulation  is  recommended  only  for  a  very  simple,  first- 
order  simulation  that  is  not  to  be  used  for  detailed  investigation  of  sophisti¬ 
cated  radars. 


E.2 


ANALOG  OPEN-LOOP  SIMULATOR 


Because  analog  computers  frequently  offer  advantages  for  some  types 
of  simulation,  it  is  reasonable  to  consider  briefly  the  possibility  of  using 
analog  rather  than  digital  simulation.  Open-loop  simulation  is  considered  in 
this  section,  and  closed-loop  simulation  in  the  following  section. 

First,  consider  the  simulation  of  the  return  from  the  target. 

Suppose  the  radar  is  a  conventional  pulse  radar.  Tbe  received  signal  is  a 
sum  of  weighted,  possibly  slightly  distorted,  time-delayed  replicas  of  the 
transmitted  pulse,  one  such  pulse  being  received  from  each  scattering  center 
on  the  target.  For  a  moving  nonfluctuating  target  a  similar  waveform  would 
result,  but  the  frequency  would  be  offset  by  the  Doppler  shift. 

Most  radar  targets  produce  returns  that  are  rapidly  fluctuating, 
because  of  changes  in  relative  phase  of  the  scattering  centers  of  the  target. 

To  a  first  order,  one  might  represent  the  fluctuations  by  a  random  process 
and  impose  noise  modulation  on  the  waveform  representing  the  signal  received 
from  a  nonfluctuating  target.  This  model  basically  assumes  that  there  are 
many,  randomly  phased  scatterers  on  the  target  and  will  not  permit  realistic 
simulation  of  correlation-function  properties  of  the  target  signal.  Thus, 
from  pulse  to  pulse  the  signal  amplitude  (pulse  level)  varies  randomly  at  a 
mean  rate  depending  on  target  size,  relative  motion  of  scattering  centers, 
and  radar  wavelength.  Depending  upon  the  instantaneous  (line-of-sight)  differ¬ 
ences  in  scattering-center  velocities,  different  Doppler  shifts  might  be  re¬ 
quired  for  the  individual  pulses.  Note  that  this  form  of  target  signal  is 
much  like  a  clutter  signal,  since  it  has  an  amplitude  that  is  at  least  par¬ 
tially,  and  may  be  entirely,  random;  the  rate  of  fluctuation  will,  in  general, 
be  quite  different,  however. 

Tbe  many-scatterer  and  random-phasing  approximations  for  scattering 
from  such  targets  as  aircraft  has  been  criticized,  largely  because  there  are 
usually  relatively  few  scatterers  on  an  aircraft,  for  example,  that  dominate 


its  radar  return  at  a  particular  time  (or  aspect  angle).  It  is  possible,  for 
some  types  of  target,  to  compute  the  radar  return  as  a  function  of  aspect  angle 
and  radar  parameters.  In  such  cases,  a  more  accurate  waveform  can  be  generated 
to  simulate  the  return  from  the  target.  However,  this  approach  probably  fits 
in  better  with  closed-loop  simulation  than  with  open-loop  simulation. 

Another  possibility  for  open-loop  simulation  of  the  target  return 
lies  in  the  recognition  that  for  many  aircraft,  as  stated  in  the  preceding 
paragraph,  significant  scattering  arises  from  a  relatively  few  scattering 
centers  that  have  random  phases.  A  random  process  having  appropriate  statistics 
(other  than  the  Rayleigh  envelope  obtained  from  many  scatterers)  can  then  be 
used  to  establish  the  output  of  the  target-signal  generator;  statistical  char¬ 
acteristics  of  such  random  processes  have  been  discussed  previously  (see 
Section  4  and  Appendix  A). 

Finally,  one  could  use  analog  recordings  of  radar  returns  from  real 
targets  or,  possibly,  from  scale  models.  In  this  way,  excellent  simulation  of 
returns  from  targets  is  possible.  The  price  paid  for  this  accuracy  is  a  loss 
in  flexibility  arising  from  the  need  for  new  recordings  for  each  change  in 
radar  parameters  and  for  various  types  of  target  motion. 

In  the  above  discussion,  analog  techniques  for  open-loop  simulation 
of  radar  returns  from  a  target  were  considered.  Consider  now  the  simulation 
of  returns  from  clutter.  If  the  clutter  arises  from  a  large  number  of  inde¬ 
pendent  scattering  elements,  none  of  which  dominates  the  rest  and  which  move 
in  such  a  way  as  to  produce  random  phasing,  the  clutter  signal  has  a  Rayleigh- 
distributed  envelope.  Clutter  within  a  particular  radar  resolution  cell  could, 
in  this  case,  be  simulated  by  white  noise  passed  through  a  filter  to  impose 
appropriate  time-correlation  properties.  If  there  were,  in  addition,  one,  or 
more,  larger  scatterers  that  caused  the  clutter-signal  envelope  to  be  non- 
Rayleigh,  these  scatterers  could  be  treated  as  additional  targets  and  the 
signals  from  them  added  to  the  noise  signal  representing  background  clutter. 

The  total  simulated  clutter  signal  then  would  be  non-Rayleigh.  In  the  limit, 
one  could  synthesize  the  complete  clutter  signal  through  the  summation  of 
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many  randomly  phased  signals  whose  amplitudes  corresponded  to  typical  clutter- 
element  returns.  Since  more  than  5  or  6  equal-amplitude,  randomly  phased 
signals  would  result  in  approximately  a  Rayleigh  envelope  distribution 
(Section  4),  a  very  large  number  of  signals  would  not  be  required,  in  practice. 
This  technique  is  analogous  to  that  used  earlier  for  digital-computer  simu¬ 
lation  wherein  and  defined  reflection  properties  of  the  scatterer. 

At  first,  the  analog  clutter-simulation  technique  appears  straight¬ 
forward.  Several  serious  difficulties  exist,  however.  First,  consider  the 
most  common  radar,  an  incoherent  pulse  radar.  Returned  pulses  from  clutter 
elements  arrive  at  different  times  because  of  the  different  distances  to  the 
elements.  The  resultant  waveform  is  generally  quite  complicated.  Consider  a 
configuration  in  which  a  signal  generator  or  oscillator  is  used  to  form  the 
clutter  signal  as  shown  in  Figure  E-l.  Randomly  phased  signals  of  appropriate 
amplitudes  (and  frequencies,  if  Doppler-shift  effects  are  to  be  included)  are 
not  alone  sufficient  to  simulate  the  received  waveform:  each  signal  must  be 
pulsed  or  sampled  so  that  its  output  has  the  pulsed-sinusoid  form  of  the 
signal  received  from  a  single  scattering  element  as  shown  in  Figure  E-l.  These 
keying  pulses  must  be  randomly  set  so  that  clutter  returns  occur  at  random 
times,  at  least  initially;  at  later  times,  it  might  be  desirable  to  let  the 
(initially  random)  pulse  times  vary  continuously,  so  as  to  pj  'vide  realistic 
spatial  and  time  correlations. 

The  technique  discussed  in  the  preceding  paragraph  would  be  suitable 
to  represent  the  return  from  stationary  clutter  in  a  single,  range -resolution 
cell  for  a  nonscanning  radar  operating  at  a  fixed  frequency.  But  suppose  the 
radar  uses  some  form  of  signal  processing  (e.g.,  integration)  that  depends 
upon  the  signals  received  from  several  transmitted  pulses.  It  is  necessary  in 
this  case  to  account  for  the  correlation  properties  of  the  signal.  At  least 
three  types  of  correlation  can  be  significant,  depending  on  the  situation,  as 
discussed  in  Section  5.  As  stated  above,  appropriately  filtered  noise  could 
be  used  in  simulation  of  single-frequency  returns  from  randomly  moving  scatter- 
ers;  here  the  filtering  operation  imposes  time  correlation  on  the  noise  volt¬ 
age  to  produce  the  correlation  that  usually  arises  from  the  finite  time 


151 


t 

PHASE 

SHIFTERS 


t 

ATTENUATORS 


MIXERS 

(S I NOLE—S I  DE-BAND 
SUPPRESSED-CARRIER) 


Figure  E-l  ANALOG  CLUTTER  SIMULATION-OPEN  LOOP 
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required  for  the  relative  phasing  of  the  scatterers  to  change  significantly. 

If  the  scatterers  and  radar  were  stationary,  however,  this  form  of  simulation 
would  fail  to  produce  the  same  clutter  signal  on  each  pulse  and  thus  would 
not  simulate  the  actual  (fixed)  clutter.  Next,  suppose  the  scatterers  were 
stationary  and  the  radar  frequency  changed  from  pulse  to  pulse,  or,  if  not 
pulse  to  pulse,  at  some  rate  producing  significant  changes  in  relative 
phases  of  scatterers  during  the  integration  time.  If  clutter  were  being 
simulated  by  band-limited  noise,  the  effect  of  a  shift  in  radar  frequency  on 
clutter  return  could  not  be  simulated,  because  the  noise  signal  cannot  include 
actual  relative  phases  of  the  scatterers.  Before  such  clutter  simulation 
could  be  used,  it  would  be  necessary  to  establish  the  correlation  properties 
of  the  clutter  signal  for  the  actual  (variable -frequency)  radar  and  to  gener¬ 
ate  a  signal  having  similar  statistical  properties. 

If  clutter  were  instead  being  simulated  by  the  output  of  a  set  of 
(separately  pulsed)  oscillators,  effects  of  scatterer  motion  could  be  included 
by  setting  the  oscillators  to  slightly  different  frequencies,  corresponding 
to  the  Doppler-shifted  frequencies  of  the  radar  returns  from  the  separate 
scatterers.  Effects  of  changes  in  radar  frequency  could  be  introduced  by 
changing  the  frequencies  of  all  of  the  oscillators  by  the  same  amount. 

A  third  form  of  correlation  that  must  be  accounted  for  if  a  signal¬ 
processing  (e.g.,  integrating)  radar  is  used  is  the  spatial  correlation  of 
clutter  returns.  Assume,  for  the  moment,  that  the  radar  antenna  is  stationary 
but  that  the  radar  returns  from  successively  farther  away  are  being  received 
(for  a  single  transmitted  pulse).  Since  this  scan  is  occurring  at  half  the 
velocity  of  propagation,  scatterer  motion  is  relatively  unimportant  for  the 
single  scan  except  for  Doppler  shifts  that  are  important  for  a  coherent  radar. 
Spatial  (in  range  direction)  correlation  now  can  arise  from  two  sources.  First, 
suppose  the  clutter-producing  scatterers  are  uniformly  and  homogeneously  dis¬ 
tributed.  Correlation  in  the  clutter  signal  now  arises  from  the  fact  that 
some  of  the  same  scatterers  produce  the  return  for  a  finite  time;  parts  of 
the  clutter  signal  separated  by  more  than  the  radar  pulse  width  are  uncorrelated. 
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This  form  of  correlation  was  discussed  in  Section  5.  Correlation  can  also 
result  from  actual  spatial  correlations  of  the  scattering  elements  (ground 
contour,  vegetation,  buildings,  etc.)*  This  form  of  correlation  can  produce 
a  clutter  signal  that  is  correlated  for  times  (ranges)  much  greater  than  a 
pulse  width.  In  practice,  of  course,  the  first  (pulse-width)  form  would  also 
occur  along  with  the  second. 

Next,  consider  range-scan  to  range-scan  correlation  (still  assuming 
a  stationary  radar).  At  a  particular  range,  from  scan  to  scan,  the  clutter 
return  will  change  only  if  the  scatterers  move  or  if  the  radar  frequency  is 
changed;  in  either  of  these  cases  the  appropriate  time  and  frequency  correla¬ 
tion  functions  should  be  satisfied,  as  discussed  above.  Thus  some  means  must 
be  used  to  permit  range-scan  to  range-scan  correlation  to  be  maintained  while 
clutter  returns  as  a  function  of  range  (during  any  single  range  scan)  are  much 
less  strongly  correlated.  If  the  clutter  signal  were  being  simulated  by  a 
noise  source  with  a  bandpass  filter  adjusted  to  produce  appropriate  scan-to- 
scan  correlation,  the  time  correlation  during  one  scan  would  then  be  much  too 
high.  If,  on  the  other  hand,  the  noise  were  filtered  to  produce  appropriate 
time  correlation  during  a  single  range  scan,  the  scan-to-scan  correlation  would 
be  essentially  zero.  To  maintain  correct  correlation  both  in  time  during  one 
range  scan  and  from  scan  to  scan  thus  requires  a  considerably  more  complicated 
implementation  than  a  simple  bandpass-filtered  noise  source.  One  way  in  which 
such  correlation  could  be  introduced  is  through  the  use  of  a  noise  source 
filtered  to  produce  appropriate  time  correlation  during  a  single  scan  feeding 
a  tapped  delay  line  or  other  delay  device,  as  shown  in  Figure  E-2;  the  (possibly 
weighted)  sum  of  the  outputs  of  the  taps  of  the  delay  line  then  is  a  simulated 
clutter  signal  with  both  the  appropriate  time  correlation  and  scan-to-scan 
correlation.  It  is  necessary,  of  course,  that  the  time  delay  between  taps  on 
the  delay  line  be  equal  to  the  time  between  scans.  The  principal  disadvantage 
of  this  technique  is  its  lack  of  flexibility;  for  example,  to  change  scan  rates 
the  delay-line  taps  must  be  shifted  correspondingly.  It  is  also  difficult  to 
allow  for  offsets  of  scatterer  motion  or  radar  frequency  shifts  from  scan  to 
scan  when  this  technique  is  used.  One  might  also  use  a  large  number  of  noise 
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/TIME  DELAY  =  SCAN-TO-SCAN  TIME 


SIMULATED  CLUTTER 


Figure  E-2  NOISE  SOURCE  CLUTTER  GENERATOR  WITH  RANGE  AND  SCAN-SCAN  CORRELATION 
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sources,  each  of  which  represents  the  clutter  signal  from  a  patch  of  ground 
that  is  smaller  than  the  resolution  cell  of  the  radar.  Each  noise  source 
would  be  filtered  to  produce  a  time  correlation  function  equivalent  to  that 
produced  by  scatterer  motion.  Effects  of  frequency  correlation  could  not 
easily  be  included,  however.  The  signal  received  during  a  range  scan  would 
now  be  found  by  adding  the  outputs  of  those  noise  generators  representing  the 
patches  of  ground  within  the  range  resolution  cell  of  the  radar;  the  noise 
generators  being  sampled  at  any  instant  would  be  changed  to  simulate  the 
motion  of  the  range  gate  along  the  ground.  This  method  of  clutter  simulation 
offers  much  greater  flexibility  than  the  one  using  delay  lines,  but  it  requires 
many  oscillators  and  filters. 

If  the  clutter  signal  is  to  be  simulated  by  adding  the  outputs  of  a 
large  number  of  randomly  phased  oscillators,  as  discussed  above  for  the  single¬ 
range  case,  range  scanning  can  be  included  directly  by  having  a  (possibly  large) 
set  of  oscillators  whose  outputs  are  pulsed  at  random  times  (corresponding  to 
scatterer  locations)  over  a  range  equivalent  to  a  range  scan.  A  range  gate 
can  be  used  to  select  the  signal  coming  from  any  desired  range.  Scan-to-scan 
correlation  is  unity  if  no  parameters  are  changed  between  simulated  scans. 

If  radar  frequency  is  changed,  relative  phasings  of  returns  from  scatterers 
are  taken  into  account  and  the  correlation  correspondingly  reduced.  Scatterer 
motion  can  be  introduced  by  slowly  changing  the  times  of  the  output  pulses  of 
appropriate  oscillators  and  changing  oscillator  frequency  to  correspond  to  the 
Doppler  shift.  Hie  multiple-oscillator  simulation  has  the  added  advantage  of 
permitting  effects  of  dominant  scatterers  to  be  included  in  a  realistic  manner, 
so  that  one  is  not  limited  to  clutter  producing  Rayleigh-distributed  signals. 

Hie  principal  disadvantage  is,  of  course,  the  large  number  of  oscillators  that 
are  required. 

In  the  discussion  above,  a  nonscanning  antenna  was  assumed.  If  the 
antenna  scans  in  azimuth,  some  new  clutter  elements  are  observed  on  each  range 
scan  (or  all  new  clutter  elements  if  the  azimuth  scan  rate  is  high  relative  to 
the  PRF) .  To  include  these  effects,  as  is  necessary  if  appropriate  spatial 
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correlation  is  to  be  retained  in  the  simulation,  requires  still  more  elabora¬ 
tion.  If  filtered  noise  sources  are  to  be  used  in  the  simulation,  ground 
patches  can  be  set  up  over  a  two-dimensional  array  and  sums  of  signals  from 
sources  corresponding  to  appropriate  elements  of  the  array  used  to  represent 
the  clutter  signal.  The  effect  of  the  antenna  pattern  can  also  be  included  by 
appropriately  weighting  the  sums.  Similarly,  if  a  set  of  randomly  phased 
oscillators  is  used  to  produce  the  clutter  signal,  as  described  in  the  pre¬ 
ceding  paragraph,  it  is  possible  to  use  a  larger  set  than  was  needed  there 
and  to  add  the  outputs  of  appropriate  members  of  the  set.  Thus,  in  essence, 
one  gives  each  of  the  randomly  phased  oscillators  in  the  azimuth  angle 
and,  for  any  particular  range  scan,  takes  a  sum  of  the  outputs  of  those  oscil¬ 
lators  corresponding  to  azimuth  angles  observed  during  that  range  scan.  Again, 
outputs  of  various  oscillators  can  be  weighted  by  a  coefficient  depending  upon 
the  antenna  patterns. 

In  the  above  discussion,  no  consideration  was  given  to  the  frequency 
at  which  the  simulator  will  operate.  Conventional  analog  computation  at  radar 
RF  frequencies  is  not  feasible.  What  considerations  are  necessary  for  scaled- 
frequency  operation?  First,  since  interference  between  pulses  of  various 
phases  and  delays  is  important,  it  seems  important  to  use  pulsed  oscillator)' 
waveforms:  the  use  of  envelopes  alone  is  not  feasible.  If  the  IF  waveform, 
rather  than  the  RF  waveform,  is  used,  frequencies  of  the  order  of  30  MHz  are 
required.  Scaling  of  time  and  frequency  to  permit  operation  at  an  a  log -computer 
frequencies  is  possible,  of  course,  but  very  slow  computation  results.  To 
permit  operation  at  30  Hz,  time  is  increased  by  a  factor  of  106,  so  a  10-sec 
(real  time)  simulation  would  require  about  4  months.  If  the  computer  could 

4 

operate  at  3  kHz,  time  would  only  be  increased  by  10  ,  but  the  10-sec  simulation 
would  still  require  28  hours.  Scaled-frequency  operation  is,  for  numbers  of 
this  size,  obviously  suitable  only  for  an  investigation  in  which  a  few  pulse 
shapes  are  to  be  found.  It  is  possible,  of  course,  that  in  some  situations 
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scaled- frequency  operation  may  prove  quite  useful:  the  numbers  given  here  are 
intended  only  to  serve  as  a  warning  that  exact  frequency  scaling  is  not  neces¬ 
sarily  an  answer  to  all  of  the  problems  of  analog-computer  simulation. 

Alternatively,  one  might  use  components  (filtered  noise  sources  or 
oscillators)  actually  operating  at  the  radar  RF  or  IF  frequency.  Conventional 
analog-computer  elements  might  be  used  to  compute  (real-time)  target  location, 
aspect  angle,  and  other  significant  parameters  which  could  be  used  to  adjust 
parameters  in  the  RF  or  IF  circuitry.  This  technique  appears  to  be  the  only 
practical  form  of  analog  computation  for  general  use,  although  in  special 
cases  it  might  be  possible  to  use  scaled  frequencies  on  a  conventional  analog 
computer  to  obtain  some  information. 

One  significant  parameter  that  has  been  omitted  in  this  discussion  is 
the  polarization  properties  of  the  received  radiation.  Radar  scattering  pro¬ 
perties  of  targets  and  clutter  sources  depend  upon  the  polarization  of  the 
transmitted  wave.  The  received  signal  also  depends  upon  the  polarization  of 
the  receiving  antenna  and  its  relationship  to  the  polarization  of  the  scattered 
wave.  One  way  to  include  polarization  diveristy  in  an  open-loop  analog  simula¬ 
tion  is  to  provide  two  separate  channels  corresponding  to  horizontal  and  verti¬ 
cal  polarization;  amplitudes  of  separate  clutter-signal  and  target -signal 
generators  would  be  fixed  at  appropriate  levels  for  their  respective  polariza¬ 
tions.  The  output  signal  from  the  antenna  (i.e.,  the  input  to  the  processor) 
would  then  be  synthesized  from  an  (appropriately  phased)  sum  of  signals  from 
the  two  channels.  The  principal  disadvantage  of  this  technique  is  the  need 
for  twice  as  many  signal-generating  components  plus  additional  complications 
arising  from  the  need  for  an  appropriate  combination  of  signals  for  the  pro¬ 
cessor  input. 

Inclusion  of  masking  effects  or  of  target-clutter  interactions  and 
bistatic-radar  effects  are  probably  not  practical  in  open-loop  simulation, 
although  they  might  be  approximately  simulated  in  some  cases.  Open-loop  simu¬ 
lation,  whether  digital  or  analog,  does  not  appear  practical  for  these  more 
complicated  situations. 
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E.3 


ANALOG-COMPUTATION  CLOSED-LOOP  SIMULATOR 


Analog  generation  of  a  signal  representing  the  radar  transmission 
poses  no  problem,  aside  from  the  question  of  what  frequency  is  to  be  used  in 
the  simulation;  as  noted  in  the  preceding  section,  frequency  scaling  must  be 
handled  with  care  because  of  the  time  requirement  involved  if  large  scaling 
factors  are  used.  For  the  moment,  assume  the  generation  of  a  pulse  at  the 
actual  radar  RF. 

Clutter  returns  can  now  be  formed  by  means  of  delay  lines  of  different 
lengths  having  variable  transmission  coefficients  or  adjustable  attenuators. 

The  variable  attenuation  permits  adjustment  of  a ^  ,  and  the  variable  delay 
adjustment  of  t \  ,  where  a ^  and  V ^  are  as  defined  in  Section  5.  Range  scan¬ 
ning  is  automatically  accounted  for  by  time  variation;  spatial  correlation  in 
the  range  direction  is  likewise  automatically  accounted  for.  Spatial  correla¬ 
tion  for  the  scanning-antenna  case  is  not  so  easily  included.  It  appears 
necessary  to  provide  sets  of  delay  lines  corresponding  to  targets  at  different 
angles  from  the  radar  and  to  switch  from  one  set  to  another,  gradually,  as  the 
(simulated)  antenna  scans  in  azimuth. 

Simulation  of  the  returns  from  the  target  can  be  effected  in  the  same 
way  as  simulation  of  the  clutter  returns,  provided  the  scattering  from  the 
target  can  be  expressed  in  terms  of  amplitudes  and  phases  of  signals  returned 
from  the  scattering  centers  on  the  target.  To  simulate  target  motion,  the 
delays  of  the  delay  lines  associated  with  the  target  scattering  centers  must 
be  changed  continuously.  Obviously  the  main  problem  in  simulating  target  and 
clutter  signals  using  this  method  is  the  need  for  many  delay  lines  producing 
low  distortion,  variable  delay,  and  variable  attenuation. 

Inclusion  of  polarization  effects  would  seem  to  require  two  sets  of 
delay  lines  and  their  controls  with  provision  for  coupling  (artificially  intro¬ 
duced  crosstalk)  of  outputs  of  these  sets  to  correspond  to  polarization  trans¬ 
formation  caused  by  the  scattering  process.  Inclusion  of  masking  effects  is 
possible  through  the  use  of  an  auxiliary  computation  of  shadowed  regions  to 
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blank  out  the  appropriate  delay-line  outputs.  Target-clutter  interaction  might 
be  treated  using  combinations  of  delay-line  outputs,  but  the  simulation  would 
be  exceedingly  complex. 

Another  technique  would  involve  the  recording  of  actual  terrain- 
clutter  radar  signals  for  a  variety  of  situations.  Quite  realistic  clutter 
signals  can  be  obtained  in  this  way,  of  course.  This  method,  while  effective 
for  the  situation  involved  in  the  recording  (i.e.,  for  the  particular  radar 
parameters  and  terrain  characteristics  used) ,  lacks  the  flexibility  that  should 
be  available  in  a  general-purpose  simulator. 
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Appendix  F 

CLUTTER  MEASUREMENTS  DATA  AND  ANALYSES 


In  this  Appendix  is  given  a  list  of  references  on  clutter  measurements 
and  analyses  which  were  used  during  the  program.  Table  F-l  shows  the  types 
of  data  available. 


Admiralty  Signal  and  Radar  Establishment  Tech.  Note  R4/50/15,  "Sea-Clutter 
Investigation,"  7  September  1950  (SECRET). 

Boeing  Airplane  Co.,  Doc.  No.  D-11738,  "An  Experimental  Investigation  of  Radar 
Ground  Clutter,"  15  January  1952  (CONFIDENTIAL). 

K.  Bullington,  "Reflection  Coefficients  of  Irregular  Terrain"  Proc.  IRE,  42, 
1258-1262  (1954). 

University  of  Illinois,  Control  Systems  Lab.  "Sea  Clutter  Spectrum  Studies 
Using  Airborne  Coherent  Radar  III"  by  B.L.  Hicks  (May  1958),  Report  No.  R-105. 

R.E.  McGarwin  and  L.J.  Maloney,  "A  Study  at  1046  MC  of  the  Reflection  Coefficient 
of  Irregular  Terrain  at  Small  Grazing  Angles"  NBS  Report  No.  5099,  25  November 
1968. 


R.L.  Cosgriff,  W.H.  Peaks,  '“.C.  Taylor,  "Electromagnetic  Reflection  Properties 
of  Natural  Surfaces  with  Applications  to  Design  of  Radars  and  Other  Sensors 
(Terrain  Handbook)"  Ohio  State  Univ.  Report  No.  694-9,  1  February  1959. 

1959  Symposium  on  Radar  Return,  Part  I  (Unc.  papers)  Univ.  of  New  Mexico, 

11-12  May  1959. 

J.  Atkins,  H.V.  Bikel,  M.  Weiss,  "Realistic  Simulation  of  Radar  Clutter" 
Electronics,  25  September  1959,  pp.  78-81. 

Goodyear  Aircraft  Corp.  "Radar  Terrain  Study  Final  Report:  Measurements  of 
Terrai  Back-Scattering  Coefficients  with  Airborne  X-band  Radar"  GERA-463, 

30  September  1959  (Unclassified). 

M.P.  Bachynski,  "Microwave  Propagation  over  Rough  Surfaces,"  RCA  Rev.  20, 
pp.  308-335  (1959). 

1959  Symposium  on  Radar  Return,  Part  2,  Confidential  papers,  11-12  May  1959, 
(CONFIDENTIAL). 
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Hughes  Aircraft  Co.  Tech.  Memo  656,  "A  Summary  of  Measurement  and  Theories 
of  Radar  Ground  Return,"  September  1960  (CONFIDENTIAL). 

J.W.  McGiven,  Jr.  and  E.W.  Pike,  "A  Study  of  Clutter  Spectra,"  in  Statistical 
Methods  in  Radio  Wave  Propagation,  Pergamon,  1960,  pp.  49-92. 

S.O.  Coleman,  G.R.  Hetrich,  "Ground  Clutter  and  the  Calculation  for  Airborne 
Pulse  Doppler  Radar"  Proceedings  of  5th  MIL-E-CON  National  Convention  on 
Military  Electronics,  26-28  June  1961,  pp.  409-415. 

Hughes  -  "Overland  AEW  Radar  Study(U),"  Final  Report  -  Report  No.  OR-F 
1  October  1961  (CONFIDENTIAL). 

"Final  Report  -  Generic  Ranging  Study,"  Vol.  I,  G.E.  Report  No.  ASER-37-61, 

15  February  1962. 

"Limited  War  Missions  Study  Program  -  Part  I  -  Problem  Definition  Requirements 
and  System  Techniques  -  Vol.  2  of  Final  Report  on  Radar  Techniques  Study, 
Development  and  Evaluation  Study,"  17  March  1962  (CONFIDENTIAL). 

J.L.  Farrell,  R.L.  Taylor,  "Doppler  Radar  Clutter,"  Trans.  IEEE,  ANE-11, 
September  1964,  pp.  162-172. 

G.R.  Corry,  "Measurements  of  UHF  and  L-Band  Radar  Clutter  in  Central  Pacific 
Ocean,"  Trans.  IEEE  MIL-9,  January  1965,  pp.  39-44. 

RRE  Memorandum  No.  2158,  "The  X-Band  Reflectivity  of  Terrain  Targets  at 
Small  Angles  of  Incidence"  by  G.  Cross  ley,  February  1965. 

"Mulitfrequency  Clutter  Study,"  Raytheon  Report  No.  BR-3289,  5  March  1965. 

L.F.  Helgostam,  B.  Ronnerstam,  "Grovnd  Clutter  Calibration  for  Airborne 
Doppler  Radars,"  Trans.  IEEE,  MIL-9,  July-October  1965,  pp.  294-297. 

Army  Missile  Command  Report,  RETR66-13,  "Analysis  of  Clutter  Data  (U)," 

9  May  1966  (SECRET) 

Tech.  Memorandum,  "Report  of  Radar  Clutter  Signal  Processing  Committee, 

Part  I,  Radar  Clutter  Effects  (U),"  JHU  Report  TG  842-1,  September  1966, 
(CONFIDENTIAL) . 

J.M.  Hunter,  TBA  Senior  "Experimental  Studies  of  Sea-Surface  Effects  on  Low 
Angle  Radars,"  Proc.  IEEE,  113,  November  1966,  pp.  1731-1740. 

"Clutter  Model  for  AEW  Radar  Design,"  U.S.  Naval  Air  Development  Center 
Report  No.  NADC-AE-6638,  29  November  1966. 

J.R.  Bamum,  "High  Frequency  Backscatter  from  Terrain  with  Buildings," 

Stanford  U.  Report  SEL-67-002,  January  1967. 
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JHU-APL  Tech.  Memo  No.  TG-899,  "Radar  Precipitation  Echoes:  Experiments  on 
Temporal,  Spatial  and  Frequency  Correlation,"  by  Nathanson  and  Reilly, 

April  1967. 

V.W.  Pidgeon,  "Radar  Land  Clutter  for  Small  Grazing  Angles  at  X-Band  and 
L-Band,"  JHU-APL  Tech.  Note,  May  1967. 

V.W.  Pidgeon,  "Time,  Frequency,  and  Spatial  Correlation  of  Radar  Sea  Return," 
JHU-APL  Tech.  Note,  May  1967. 

N.W.  Guinard,  J.T.  Ransome,  Jr.,  M.B.  Laing  and  L.W.  Hearten,"  NRL  Terrain 
Clutter  Study  Phase  I,  NRL  Report  6487,  May  10,  1967. 

T.A.  Croft,  "Computation  of  HF  Ground  Backscatter  Amplitude,"  Radio  Science, 
VOL.  2  (New  Series)  July  1967,  pp.  739-746. 

T.R.  Benedict  and  T.T.  Soong,  "The  Joint  Estimation  of  Signal  and  Noise  from 
the  Sum  Envelope,"  Trans,  IEEE  IT-13,  pp.  447-454,  July  1967. 

D.F.  DeLone,  Jr.  and  E.M.  Hofstatter,  "On  the  Design  of  Optimum  Radar 
Waveforms  for  Clutter  Rejection,"  Trans.  IEEE,  IT-13,  pp.  454-463,  July  1967. 

J.  Kroszezyinski,  "Efficienty  of  Attenuation  of  Moving  Clutter,"  Radio  and 
Electronics  Engineer  34,  pp.  157-159,  September  1967. 

Reports  on  Masking  and  Clutter,  Including  Information  on  HAWK  Sites: 

1.  Allied  Air  Forces  Central  Europe  Rept.  0RA/72/AMB  (1964) 

Radar  Screening  Data  for  Hawk  Missile  Sites,  Wurzburg  (U) , 

NATO  SECRET  Report. 

2.  Allied  Air  Forces  Central  Europe  Rept.  0RA/74/AMB  (1964) 

Radar  Screening  Data  for  Hawk  Missile  Sites,  Kitzingen  (U), 

NATO  SECRET  Report. 

3.  AIRCENT  (Allied  Air  Forces  Central  Europe)  Rept.  0RA/39/AIL, 

19  Oct.  1962,  NATO  CONFIDENTIAL  Report. 

4.  AIRCENT  Rept.  0RA/66/AMB,  6  June  1966,  NATO  SECRET  Report. 

AIRCENT  Rept.  ORA/80/AMB,  25  June  1965,  NATO  SECRET  Report. 

6.  AIRCENT  Rept.  ORA/82/AMB,  17  May  1966,  NATO  SECRET  Report. 

7.  AIRCENT  Rept.  AEV-1730/6,  29  July  1965,  NATO  SECRET  Report. 

8.  Royal  Armament  Res.  6  Dev.  Establishment  Memo  2/64,  no  date 
(rec'd  CAL  Library  15  July  1964),  CONFIDENTIAL  Report. 
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9.  RRE  (Royal  Radar  Establishment)  Memorandum  No.  1592,  26  May  1959, 
NATO  SECRET  Report. 

10.  RRE  Memorandum  No.  1791,  December  1960,  NATO  SECRET  Report. 

11.  RRE  Memorandum  No.  TIL/BR/7365,  7  June  1964,  SECRET  Report. 
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Appendix  G 

MEASUREMENT  OF  FORWARD-SCATTERING  USING 
HIGHLY  r TABLE  FREQUENCY  SOURCES 


Because  of  the  long  baselines  involved,  coherency  measurements  of  a 
forward  scattered  microwave  transmission  using  conventional  single  source 
interferometry  techniques  may  not  be  feasible.  An  alternative  approach  involves 
the  use  of  separate  sources  which  have  sufficient  frequency  stability  (or 
spectral  purity)  so  that  phase  variations  between  them  during  the  period  of 
observation  are  acceptably  small.  This  technique  has  been  successfully  applied 
to  radio -astronomical  interferometry  where  rubidium  atomic  frequency  sources 
allowed  observation  times  of  2-1/2  minutes  without  loss  of  coherence 
(Reference  1).  For  the  purpose  of  microwave  forward  scattering  measurements, 
however,  observation  periods  over  which  coherent  measurements  are  desired  will 
probably  not  exceed  a  few  seconds.  In  this  case  high-stability  quartz  crystal 
oscillators  would  perform  satisfactorily.  Typical  rms  phase  variation  during 
a  1  second  observation  period  for  a  filtered  5  MHz  quartz  crystal  oscillator 
multiplied  to  10  GHz  is  about  0.2  radians,  from  short-term  frequency  stability 
data  presented  in  References  2  and  3.  It  may  also  be  noted  from  these  data 
that,  although  the  phase  stability  deteriorates  for  observation  periods  greater 
than  1  second  due  to  1/f  noise,  it  generally  improves  or  remains  constant  for 
shorter  periods. 

From  the  spectral  viewpoint,  these  phase  instabilities  may  be  char¬ 
acterized  as  a  wide  band  noise  spectrum  surrounding  the  carrier  frequency. 

This  noise  spectrum  will  limit  the  usable  dynamic  range  of  forward  scatter 
spectral  measurements.  At  10  GHz,  this  limit  is  approximately  45  dB  down  from 
the  carrier  level  for  a  10  Hz  analyzer  bandwidth. 

Since  both  phase  stability  and  S/N  ratio  are  inversely  proportioned 
to  the  frequency  multiplication  factor  (References  2  and  3) ,  tests  conducted 
at  frequencies  lower  than  10  GHz  will  yield  higher  phase  accuracy  and  improved 
spectral  dynamic  range.  Thus,  coherent  forward  scattering  measurements  are 
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quite  feasible  using  independent,  high-stability  quartz  crystal  oscillators 
multiplied  to  microwave  frequencies  up  to  10  GHz. 
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13  ABSTRACT 

The  objective  of  this  program  was  to  prepare  a  clutter  research  plan  for  Army 
air  defense  weapons  that  encompasses  the  anticipated  methods  of  employment  of  these 
weapons,  as  well  as  the  problems  associated  with  the  introduction  of  new  radar 
technology  in  these  systems.  Specifically,  the  investigations  included  studies  of 
the  physics  of  electromagnetic  scattering  from  rough  surfaces,  a  definition  of  clutter 
problems  as  related  to  innovations  in  radar  technology  such  as  electronic  beam 
steering,  signal  processing  and  design,  computer  technology,  study  of  clutter  models 
and  simulation  techniques,  and  a  study  of  clutter  effects  in  a  radar  defense  complex. 

Probability  density  functions,  correlation  functions  and  doppler  spectra  of 
clutter  signals  are  derived  for  the  configuration  in  which  only  a  few  major  scatterei; 
are  present.  Results  are  presented  for  both  coherent  and  noncoherent  detection. 

Methods  of  simulating  clutter,  using  both  analog  and  digital  approaches,  are 
described.  Emphasis  is  placed  on  advanced  forms  of  simulators  such  as  would  be 
required  for  evaluation  of  future  multi-dimensional  radar  systems  rather  than  on 
additive  filtered  noise  for  the  simulation  of  clutter.  The  simulation  techniques  are 
based  on  two  suggested  generic  forms  of  clutter  models,  called  open-loop  and  closed 
loop.  Because  of  the  greater  flexibility  of  closed-loop  simulation,  this  is  presented 
in  greater  detail . 

Techniques  for  performing  radar  clutter  measurements  are  described  and  some  of 
the  significant  problems,  related  primarily  to  terrain  and  cultural  observables,  are 
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Abstract  (Cont.) 

discussed.  Clutter  problems  as  related  to  radar 
systems  and  innovations  in  radar  technology,  such 
as  electronically  steered  phase  arrays,  are  also 
described. 


